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PREF ACE 


The purpose of this study is to analyze and evaluate the factors 
involved in the determination of modern rates of postglacial crustal 
movement in the Great Lakes region in order to ascertain the validity 
of previously determined rates of uplift. 

Because the calculation of uplift rates are based on findings from 
many fields, i.e., geology, geophysics, geodesy, meteorology, oceanog- 
raphy, and engineering, it is necessary to integrate information from 
all of these fields if a solution to the problem is to, be tounge 
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ABS TRACT 


The purpose of this study is to analyze and evaluate the factors 
which are involved in the determination of modern rates of postglacial 
crustal movement in the Great Lakes region in order to ascertain the 
validity of previously determined rates of uplift. 


Because the calculations of rates of uplift are based on findings 
from many fields, i.e., geology, geophysics, geodesy, meteorology, ocea- 
nography, and engineering, it is necessary to integrate information from 
all of these fields if a solution. to the problem is to. be found, 


The problem was studied by the following methods: 

(a) by analyzing the data, methods, and results of previous investiga- 
tors of modern uplift; (b) by comparing modern areas and rates of up- 
lift with areas and rates of uplift based on differential warping of 
former glacial lake shoreline features; (c) by comparing modern uplift 
in the Great Lakes region with modern uplift in Fennoscandia; (d) by 
examining principles of operation, instrument construction, external 
influences, and errors which are inherent in water-level gaging (the 
results of gaging are used to calculate modern rates of uplift); (e) 

by computing daily, monthly and summer season vector winds for Lakes 
Erie, Ontario and Superior (1950-59) in order to test the assumption 
(underlying water-leveling and the calculation of rates of uplift) that 
the summer mean Great Lakes water surfaces are level; and (f) by mak- 
ing a correlation study of effective wind velocities and lake-level 
gage differences to determine if gage differences represent land uplift, 
or wind slope of the lake surface. 

The results indicate that: (a) all previously determined modern 
rates .of uplift on Lake Erie are not valid because Lake Erie is in (the 
postglacial area of horizontality; (b) all rates of modern uplift based 
on pairs of gages in which one gage is south of the Nipissing zero iso- 
base.: ate also not valid; and (c) rates of uplift on Lake Superior and 
on the other Great Lakes which are based on records of gage pairs located 
north of the Nipissing zero isobase are probably erroneous owing to the 
inclusion in the gage differences of errors caused by meteorological 
effects, ‘gage location efieets, instrument error and operator error. 


Six of the more important conclusions are: 
1. lLake-level gages used to detect and measure modern uplift must be 
north of the known area of horizontality (i.e., north of the Nipissing 
zero isobase). Gage records must be corrected for meteorological, in- 


Vase 


strument and operator errors. 

2. The Lake Erie correlation study showed that gage differences repre- 
sent, almost wholly, meteorological effects and not uplift. 

4. The summer season mean water surfaces of the Great Lakes are not 
level. 

4, Ekman's theory of ocean currents should be re-examined and modified 
on the basis of empirical observations. 

5. The Lake Erie correlation study indicated that Ekman's concept of 
water-surface slope direction being in the direction of the wind is in- 
correct; on Lake Erie the water-surface slope direction is about 23° to 
the right of the wind. 

6. The angle of deviation of surface currents from wind directions on 
Lake Erie (c 23° to the right) should also be representative of the angle 
of deviation for the other Great Lakes. ) 
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\ I. HISTORICAL INTRODUCTION 


\ 
\ 


\/ Fennoscandia 


The phenomenon of land rising from a body of water has been a sub- 
Ject of study and speculation since Classical times. Although Greek, 
Roman and Medieval writers recorded the uplift of land in volcanic areas, 
it was not until 1625 that slow imperceptible uplift of land in a non- 
volcanic region was mentioned in a Finnish book of discourses. 

The gradual seaward retreat of the shorelines of the shallow Baltic 
Sea and the Gulf of Bothnia, as well as the shoaling of harbors and water- 
ways, caused Scandinavian naturalists of the late 17th and early 18th cen- 
turies to speculate as to its cause. The question soon arose as to 
whether the land in the Baltic area was rising, or the quantity of water 
in the sea diminishing. Such noted scientists as Linnaeus, Celsius and 
Swedenborg held that the waters of the Baltic Sea were decreasing. 

In an effort to actually measure the change between land and sea, 
Anders Celsius made what is probably the first systematic observations 
of this relative change when, in 1/31, he had marks corresponding to sea 
level and the date chiseled into exposed bedrock. By comparing later 
marks with the 1731 mark, Celsius estimated the lowering of sea level to 


be 4.1 feet per century. (Kiuridinen, 1953; p. 7). 


The theory of diminishing water was shown to be false by the Bishop 
of Abo, Johan Browallius, in 1/55 when he pointed out that the sea main- 
tains an equilibrium and that churches built in Sweden in the 14th cen- 
tury were still near the shore and ancient buildings in southern Swede 
and Denmark were still on the shore; therefore @ general decrease cies 
level could not have occurred. (Ka#ridinen, 1953, p. 8; Thorarinsson, 
1940, p. 1423 Bergsten, 150.7120). 

HK. O. Runeberg of Finland was the first proponent of crustal elée- 
vation as the cause of the shifting of the coastline, of the Gulf of 
Bothnia. His observations and studies led him to declare in 1765 that 
the earth's crust rises while sea level remains the same. (Ka#ridinen, 
1953, p- &; Bergsten, 1940, p. 21,) This view was also supported by 
Playfair (1802, pp. 445-447) in his Illustrations of the Huttonian 
Theory and by Von Buch (1813, p. 386) in his Travels through Norway and 

Charles Lyell (1837, pp: 437-449) discussed the uplift of land in 
the Baltic region in "Chapter XVII, Elevation and Subsidence of Land 
without Earthquakes" in his Principles of Geology. This chapter, which 
summarized the investigations and observations which had been made be- 
fore 1837, concluded with remarks as to the significance of a slow sec- 
ular uplift and speculated as to its cause. Lyell's last paragraph of 
Chapter XVII supplies the views of the possible causes of the Baltic up- 
lift which were current in the first half of the 19th century. , He staves 


(1837, p. 449); 


The foundations of the country, thus gradually uplifted 
in Sweden, must be undergoing important modifications. Whether 
we ascribe these to an expansion of solid matter by continually 
increasing heat, Or to the diquefaction of rock, or to the crys- 
tallization of a dense fluid, or the accumulation of pent-up 
gases, in whatever conjecture we indulge, we can never doubt for 
amoment, that at some unknown depth the structure of the globe 
is in our own times becoming changed from day to day, through- 
out a space probably more than a thousand miles in length and 
several hundred in breadth. 


The modern theory of the cause of crustal movement in Fennoscandia 


and North America, i.e., ice unloading and subsequent isostatic adjust- 


ment, 


was proposed by a Scot, Thomas F. Jamieson, in 1865. His theory, 


quoted in the following paragraph, was later developed in detail in the 


paper "On the Cause of the Depression and Re-elevation of the Land during 


the Glacial Period" published in The Geological Magazine in 1882. 


that: 


tial 


Jamieson (1865, p. 178) stated in his first expression of the theory 


It is worthy of remark that in Scandinavia and North 
America, as well as in Scotland, we have evidence of the 
great ice-covering; and singular to say, the height to which 
merine fossils have been found in all three countries is very 
nearly the same. It has occurred to me that the enormous 
weight of the ice thrown upon the land may have had something 
to do with this depression. Agassiz considers the ice to have 
been a mile thick in some parts of America; and everything 
points to a great thickness in Scandinavia and North Britain. 
We don't know what is the state of the matter on which the 
solid crust of the earth reposes. If it is in a state of 
fusion, a depression might take place from a cause of this 
kind, and then the melting of the ice would account for the 
rising of the land, which seems to have followed upon the 
decrease of the glaciers. 


Fennoscandia continued to provide material for research on differen- 


uplift of ‘the land until, according to E. Ka&ridinen (1953, p- 9), 


fifty or more scientists had investigated uplift in this region before 

the beginning of the twentieth century. Since 1900 more than forty Euro- 
pean and British pone eee geodesists and engineers have studied and re- 
ported upon post-glacial and contemporary differential land uplift in the 
Baltic region. 

The importance of Fennoscandian investigations of this secular up- 
lift to similar studies in the Great Lakes area may be briefly stated 
the Fennoscandian studies pioneered in the use of water gaging and pre- 
cise leveling to reveal uplift, in applying corrections (chiefly meteor- 
ological) to water gaging in order to reduce systematic errors and in 
supplying the necessary mechanism for the cause of the uplift. 

These Kuropean studies, published over a period of more than 200 
years, were probably the source of American ideas on this subjects” ius 
thermore, the current Scandinavian literature dealing with modern crustal 
deformation reveals many concepts and techniques, particularly in the 
application of corrections for systematic errors, which have not been 


used in American studies of the same topic. 


Great Lakes Region 


Detection of crustal movements in the Great Lakes region has been 
hampered by the large area involved, by the sparse population (conse- 
quently by the limited number of observers), and by the relatively short 


time that the region has been settled. The limited distribution of obser- 


vation points coupled with the short period of recorded observation may 
help to explain why the first systematic study of "modern" land uplift 
was not undertaken until 1896 by Dr. G. K: Gilbert. 

Three members of the early geological surveys of New York, Ohio and 
Michigan reported ancient uplifted beaches. Ebenezer Emmons (1837, p. 
123), geologist for the New York Natural History Survey, described con- 
current land uplift along the St. Lawrence River; he compared it to the 
uplift in Norway and said that the only way of discovering its magnitude 
was to establish fixed land marks which would be compared over a number 
of years. Emmons (1838, p. 239) also discussed uplift in the Lake 
Champlain area. Charles Whittlesey (1838, p. 55), topographer of the 
Ohio Geological Survey, reported measuring ancient tilted beaches south 
of the Lake Erie shore in 1838, and Bela Hubbard (1840, pp. 105-106) of 
the Michigan Geological Survey described the formation of raised beaches 
near Lake Erie as follows: 

»»-ln Other words, the land has been subsequently sub- 
jected to an upheaving force, which at last has elevated the 
whole far above the influence of the sea. 

Whether the upheaving of the land was general at this 

era, throughout the continent or was mainly operative in the 

region of the lakes, probably cannot be satisfactorily deter- 

mined. It may be competent, however, to suppose that these 
apparent "lake ridges" were the boundaries of the ancient sea 


formed during intervals of rest in the upward tendency of the 
land. 

There also exist strong reasons for supposing that the 
relative levels of the land did not everywhere remain the 
same, or that disproportionate elevations took place.... 


Despite the early recognition of tilted beaches in the lower Great 
Lakes region it was not until the latter part of the 19th century that 
the studies of deformed postglacial shoreline features by Bell,” Gilverr 
Goldthwait, Lawson, Leverett, Spencer, Taylor and Upham provided the 
proper "atmosphere" for the consideration of modern uplift. 

GOR Stuntz, a Wisconsin surveyor, submitted a paper to the Pee 
can Association for the Advancement of Science in 1854 in which he sug- 
gested that the waters of Lake Superior appeared to be rising at the 
west end of the lake and falling at the eastern end. This papery ware 
was published in 1869, was one of the first to infer that a change was 
taking place in modern times in the relative positions of the lake water 
and the land (1869, pp. 205-210). 

In 1868 (p. 129), N. S. Shaler, while discussing changes of level 


! 


of shorelines. said: .--LOOking still further, we perceive some vere 
peculiar features in the distribution of the changes of level which are 
still going forward, or which have taken place since the close of the 
glacial period." 

The next suggestion that modern land uplift could be occurring was 
made by Robert Bell in his description of his explorations of the Hudsom 
Bay area. Despite the fact that Bell later (1897) spoke of the "Rising 


! 


of the Land around Hudson Bay," in his earlier reports he was not as def- 
inite in asserting that uplift was occurring. Bell's report published in 


1880 spoke of: 


..-the comparatively rapid elevation of the land, or retiring 

of the sea, around James Bay and at York Factory was referred 

to in my reports for 1877 and 1878. ... ‘This recession of 

the sea may be due to a general lowering of its level rela- 

tively to the land, and partly to the silting up of portions 

of Hudson's Bay, interrupting the free flow of the tides (p. 

2 Veu,. 

The first definite statement regarding modern crustal movement in 
the Great Lakes area was made by J. W. Spencer (1894) in a paper "The 
Duration of Niagara Falls." He (p. 472) concluded his paper by declar- 
ing: "...Lastly, if the rate of terrestrial deformation continues as 
it appears to have done, then in about 5000 years the life of Niagara 
Falls will cease by the turning of the waters into the Mississippi." 
Spencer (1907, 1913) later recanted this hypothesis (the same idea has 
been expressed by other writers several times since his original utter- 
ance) after reworking Gilbert's data in conjunction with additional in- 
formation. 

Grove Karl Gilbert (1896-97) made the first determination of land 
uplift in the Great Lakes region using "modern" observations. The under- 
lying principles of his procedures and techniques have been followed by 
all subsequent investigators. Gilbert's original determination of the 
amount of earth tilting was prompted by the speculation that the forces 
which had tilted former beaches were still active and could be detected. 
He made use of the U. S. Lake Survey assumption that a lake surface is 
level if measured over a "protracted" period to provide the necessary 


level for comparison between pairs of gages located in the general direc- 


LOMO ale 


Four pairs of water-level gages were used and the rate of tilting 
was found by comparing the gage differences for each pair of gages for 
1874 with the gage differences of the same pairs for 1896. Because the 
gage differences were compared in relation to the same datum (the level 
lake surface) they would have been the same in 1896 as in 1874 if no 
tilting had occurred. However, the gage differences for 1896 avereeeele 
from those of 1874 which caused Gilbert to conclude that tilting had 
occurred during the intervening’ years. 

Gilbert's awareness of the pitfalls which were present in his as- 
sumptions, raw data, and method of calculation induced him to include a 
discussion of the sources and importance of errors inherent in the data 
and method. His report also included a section on "Plans for Precise 
Measurement," which, if followed, would have greatly reduced or elimi- 
nated the systematic errors whitch are still intrinsic in waver—-Lever 
records. 

The next two studies of modern crustal movement were made by 
J. W. Spencer in 190/ and 1914. Spencer applied Gilbert's techniques 
to the study of gage records taken from 1855 to 1905 in the first paper 
and from 1855 to 1912 in the second paper. He calculated the gage dif- 
ferences for periods of five years (using the mean daily lake levels), 
and concluded from his study that the means of all the five year periods 
from 1855 to 1912 were within the limits of the probable error—therefore 


the earth's crust was stable. 


These investigations were followed by a paper in 1922 by Sherman 
Moore, an engineer of the U. S. Lake Survey. Again using the procedure 
originated by Gilbert, Moore examined 18 pairs of gages whose periods 
of record were from 1870-80 to 1919. The gages were located on all of 
the Great Lakes, including Lake Superior. Moore's determination of the 
rates of uplift for Lake Superior were the first calculated for that lake. 
His conclusion was that tilting of the land occurs on all of the lakes 
but that the rates vary for different lake basins. Moore's estimate of 
the general rate of tilting was about six inches in one hundred miles in 
one hundred years. 

The controversy arising from the diversion of water from Lake Mich- 
igan by the City of Chicago led to two studies of the hydrology of the 
Great Lakes. ‘These studies, one by John R. Freeman in 1926 and the 
other by R. E. Horton and C. E. Grunsky in 1927, included determinations 
of land uplift. Both papers reviewed the literature dealing with this 
topic and both recalculated the rate of movement. The results were 
essentially the same as had been found by previous investigators, al- 
though Horton and Grunsky concluded that no uplift was occurring around 
Lake Erie. 

/ Beno Gutenberg (1944) conducted an investigation of "Tilting Due 
to Glacial Melting," including not only the Great Lakes region but also 
Hudson Bay and Scandinavia. His paper included a compilation of previous 


work on crustal deformation, the calculation of rates of movement on 


LO 


the Great Lakes from 28 pairs of gages, and a discussion of the deter- 
mination of changes of level by means of the examination of ocean tide 
gage records. His’ results” concurred with those of Gilbert, Moore rand 
Freeman. 

Bight years after his first paper on postglacial tilting, B. Guten= 
berg (1941) published a much more comprehensive work "Changes in Sea 
Level, Postglacial Uplift, and Mobility of the Earth's Interior." Using 
additional lake-level gages and tide gages, Gutenberg made new calculations 
of the rates of uplift in both North American and Fennoscandia, which he 
demonstrated could be used to help determine the viscosity and strength 
Of ime eartias 

~/ Sherman Moore's (1948) second paper on "Crustal Movement in the 
Great Lakes Area" was the result of data which had been accumulated in 
order to establish a new datum (1935 Datum) on the Great Lakes. 

The re-leveling, both water and instrumental, which was necessary 
to establish the new datum permitted Moore to determine the rates of up- 
lift in relation to sea level. Moore's analysis of the data led to a 
number” of Interesting econelusions, several of which do not concum witee 
current geologic thought. Moore (p. 697) inferred that, "The entire area, 
except for the extreme northerly part of Lake Superior is subsiding with 
respect to sea level." He (p. 708) also stated, "Whatever the cause of 
the postglacial warping, there seems to be no connection between the pres- 


ent movement and isostatic recovery from the weight of the ice." 


ay 


In 1954 Charles A. Price of the Canadian Hydrographic Service re- 
ported on the "Crustal Movement in the Lake Ontario—Upper St. Lawrence 


River Basin." 


The procedures which he used were those employed by pre- 
vious investigators. Price (1954, Pl. 0-6) reported the over-all change 
in gage relations as 0.5% feet per 100 miles per 100 years with tilting 
occurring in a direction N 40° E. 

The most recent determination of rates of crustal movement in the 
Great Lakes region are those of the Vertical Control Subcommittee of the 
Coordinating Committee on Great Lakes Basic Hydraulic and Hydrologic 
Data. The Subcommittee, which is composed of Canadian and American mem- 
bers, was appointed to establish a new level datum for the Great Lakes. 
Part of its duties included the study of crustal movements. The results 
of the Subcommittee findings have not been published although interim 
reports have béen issued to participating government agencies. 

The techniques which are used by the Subcommittee to calculate rates 
of movement are those which were originated by Gilbert in 1896. Although 
certain minor specifications were made as to the data to be used, the 
Subcommittee has not analyzed the concepts or assumptions which are in- 
herent in the calculation of rates of movement by the lake-level gage 
procedure; therefore it may be expected that its findings will agree for 


the most part with prior determinations. | 


Lite SGAND UEEiIET 


"Land uplift" which is defined as, "Elevation of any extensive part 
of the earth's surface relatively to some other part...." (A.G.1.,7i¢5m 
p. 310) will be restricted in this study to the differential vertical | 
movement of land in the former glaciated areas of Fennoscandia and North 
America. 

The existence of this land uplift has been revealed by the measumes 
ment of warped former shoreline features of late and postglacial water 
bodies; by examination of water-level gage records of modern lakes and 
seas; and by precise léveling-: 

As the ice sheets of Late Wisconsin time retreated in the Great 
Lakes region toward the center of glaciation in, the area. of Hudson Ba 
large lakes were formed by glacier ice blocking the normal drainage of 
the land. The waters of these glacier-margin lakes rose until they were 
able to escape to the south over low areas of the Great Lakes watershed. 
As the glacier front retreated, or occasionally advanced slightly, var- 
ious. cutlets were covered ‘or uncovered, which allowed the lakes to vary 
in size and elevation. ‘The lake levels remained constant long enough 
for shoreline features, i.e., beaches, wave-cut' cliffs, deltas, to be 
developed on the emerging shores. Shoreline features which were level 


when first built are now found to be warped upward toward the northeast 


2 


13 


in the northern part of the lake basins. The greatest tilt is measured 
in the highest (oldest) shorelines. 

This progressive northward warping of the shorelines is due to a 
differential uplift of the earth's crust which most geologists believe 
to be caused by recovery of the crust after depression by the load of 
the glacial ice. 

The tendency of portions of the earth's crust to approach a con- 
dition of balance leads to the establishment of a state of equilibrium in 
the crust known as "isostasy." The upwarping of strandlines in former 
glaciated areas has been cited by many geologists and geophysicists as 


one of the most convincing proofs of the principle of isostasy. 


Mechanics of Warping 


The determination of the underlying cause of warped shoreline fea- 
tures in Fennoscandia and North America has unfortunately been compli- 
cated by the fact that the glaciated areas correspond very closely with 
the Baltic Shield and the Canadian Shield, regions where ene movements 
have been upward for a long period of time (see Plate I). 

Scientists who explain the upwarping by isostatic rebound are opposed 
- by those who believe that modern uplift is a continuation of movements 
which have characterized shield areas since Precambrian times. One of 
the adherents of the tectonic, endogenetic theory of land uplift in these 


areas summarized this concept by stating: 
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That both the Canadian and the Baltic shields were being ele- 
vated, clearly by tectonic forees, prior to the siseral, period, 
and there are no grounds for maintaining that the very same 
forces did not play an important role in the recent movement 
of the’ shields ((Lyustikh, 1960;.p. 107): 


The arguments of the proponents for a tectonic cause of uplift were 
answered, for the most part, by R. A. Daly (1940, p. 318) who declared: 


This hypothesis, that there is no, connection between the 
upwarping and the deglaciation shows its full weakness when _ 
confronted with field statistics. We might conceive that the 
observed systematic warping in one or two tracts might be ex- 
plained by independent epeirogenic movements, but it seems 
incredible that in a dozen regions the same type of warping 
should appear as mere accidental products of a stress system 
that has no vital connection with ice-loads. Yet basining 
and recoil have been demonstrated in as many widely separated 
tracts, each having been covered with heavy masses of ice re- 
cently melted away. In some, if not ell, (of whe Gases the 
melting and warping began less than 50,000 years ago. With 
few exceptions there are no signs that the lithosphere out- 
Side these tracts was simultaneously disturbed by anything 
like the same amount. 


Isostasy 


The first clear explanation of the theory of isostasy came out of an 
attempt to explain certain systematic errors in the calculations of the 
Trigonometrical Survey of India in the first half of the nineteenth cen- 
tury. When latitudes and longitudes derived from astronomical observa- 
tions were compared with latitudes and longitudes of the same stations 
computed from the triangulation net, it was found that the relative posi- 
tion determined by triangulation failed to agree with the astronomical 


positions by an amount which could not be accounted for by surveying error. 
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The determinations of astronomic positions are made with instruments 
containing leveling bubbles; therefore the vertical axis of the instru- 
Menveds perpendicular to the geoid; i.e., the vertical axis, like a plumb 
Sige, 1s parallel to the direction of the force of gravity. On the other 
hand, triangulation positions are calculated on the basis of an assumed 
ellipsoid and yield geodetic latitudes and longitudes; perpendiculars to 
these coordinates are normal to the ellipsoid. The angle between the 
normal to the ellipsoid and the normal to the geoid is called "the de- 
flection of the vertical." 

It has long been known that a plumb line (normal to the geoid) is 
attracted by the mass of a mountain (Newton, 1728, p. 41) and that the 
magnitude of attraction is less than would be expected if the attractive 
force depended solely upon the mountain's mass as computed by its dimen- 
sions and the average density of rock (Bouguer, 1749, pp. 364-385; 
este: 17755:pp. 500-508; DeZach, 1814, pp. V-XVI). However, it was 
not until 1855 when two papers dealing with this subject were published 
in the Philosophical Transactions of the Royal Society of London that an 
explanation for this phenomenon was advanced. 

John H. Pratt, an archdeacon of Calcutta, discussed the systematic 
error which had been revealed in a comparison of the astronomic and geo- 
detic determinations of latitude of the terminal points of the northern 
division of the Meridional Arc of India. After proving that the geodetic 


calculations were not at fault, Pratt pointed out that the plumb line of 
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the astronomical observations must have been affected by the mass of the 
Himalayas and the Tibetan plateau. He computed the magnitude of the at- 
traction by, means of a laborious. mechanical adintegration of. the 
topographic effect of the mass of the Himalayas and the mountain region 
beyond. His computations showed that the difference between the calcu- 
lated deflections of the plumb line at Kalianpur on the central plateau 
and Kaliana at the foot of the Himalayas was three times larger than the 
measured deflection. 

Although Pratt (1855, p. 100) declared, "The conclusion, then, to 
which [ come is, that there is no way of reconciling the ditferences as 
tween the error. in latitude deduced in Colonel Everest's work and the 
amount I have assigned to deflection of the plumb line arising from 


attraction. a 


the next paper in the same volume of the Transaction 
provided the answer to the problem. 

George B. Airy, Astronomer Royal of Great Britain, heard the orig- 
inal presentation of Rev. Pratt's paper on December 7, 1854, and, in 1855, 
submitted &@ short paper to the Royal Society which expressed for the times 
time the principles of that condition.of equilibrium in the earth's crust 
which was later called "isostasy." In his succinct paper Airy explained 
that the fluidity of the earth's interior was imperfect; that it was 
probably extremely viscous; that the material below the crust was of 


greater density than the crust; and that the strength of the crust was 


insufficient to support the weight of a table-land or to maintain a moun- 


tain range, although it would support local topographic features such as 


a) 


a mountain. He (1855, p. 163) expressed his views of the underlying sup- 


port for elevated regions of the earth's crust as follows: 


I conceive that there can be no other support than that aris- 
ing from the downward projection of a portion of the earth's 
Pigniecrussenvo- Loe dense lavas... the depth of its: projec- 
tion downwards being such that the increased power of floatation 
thus gained is roughly equal to the increase of weight above 
from the prominence of the table-land. 


In regard to the effect on gravity of the elevated portion of the 


crust and its downward projection, Airy (1855, pp. 103-104) said: 


Let us consider what will be its effect in disturbing the 
direction of gravity at different points in its proximity. It 
will be remarked that the disturbance depends on two actions; 
the positive attraction produced by the elevated table-land; 
and the diminution of attraction, or negative attraction, pro- 
duced by the substitution of a certain amount of light crust 
(in the lower projection) for heavy lava. 

The diminution of attractive matter below, produced by 
the substitution of light crust for heavy lava, will be sen- 
sibly equal to the increase of attractive matter above. The 
difference of the negative attraction of one and the positive 
attraction of the other, as estimated in the direction of a 
line perpendicular to that joining the centers of attraction 
of the two masses (or as estimated in a horizontal line), will 
be proportional to the difference of the inverse cubes of the 
distances of the attracted point from the two masses. 

The general conclusion then if this, In all cases, the 
real disturbance will be less than that found by computing 
the effects of the mountains, on the law of gravitation. Near 
to the elevated country, the part which is to be subtracted from 
the computed effect is a small proportion of the whole. Ata 
distance from the elevated country, the part which is to be 
subtracted is so nearly equal to the whole, that the remainder 
may be neglected as insignificant, even in cases where the 
attraction of the elevated country itself would be considerable. 
But in our ignorance of the depth at which the downward immer- 
sion of the projecting crust into the lava takes place, we 
cannot give greater precision to the statement. 
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Airy recognized that isostasy may be incomplete—this is revealed 
in his statement, (1655, p.. 104) that: 


In all the latter inferences, it is supposed that the 
Crustiiss floating saimeerstate: of-equilibrium. BUC ainRoOuanen= 
tire ignorance of the modus operandi of the forces which have 
raised submarine strata to the tops of high mountains, we can- 
‘not insist on this as absolutely true. We know (from the 
reasoning above) that it will be so to the limits of breakage 
[strength of the crust] of the table-lands; but within those 
limits there may be some range of the conditions either way. 
It is quite possible that the immersion of the lower projec- 
tion in the lava may be too great, as that the elevation may 
be too great; and in the former of these cases, the attrac- 
tion would be negative. 


Four years after Airy made public his theory of the equilibrium of 
the earth's crust, J. H. Pratt (1859) published a second paper on the 
deflection of the plumb line in India which supplies the rudiments for 
the hypothesis known as "Pratt's Theory of Isostasy." 

Although Pratt admitted in the beginning of his paper that Airy's 
concept of a deficiency of mass beneath the mountain mass was correct, 
he rejected Airy's explanation of the phenomenon and proposed one of his 
own. Pratt (1859, pp. 746-747) states his objection to Airy's theory 
as follows: 

This hypothesis appears, however, to be untenable for 

three reasons: (1) It supposes the thickness of the earth's 

solid erust to be considerably smaller than that assigned by 

the only satisfactory physical calculation made on the sub- 

ject—those by Mr. Hopkins of Cambridge. He [Mr. Hopkins ] 

considers the thickness to be about 800 or 1000 miles at 

least. (2) It assumes that this thin crust is lighter than 

the fluid on which it is suppose to rest. But we should ex- 

pect that in becoming solid from the fluid state, it would 


contract from loss of heat and become heavier. (3) ‘The same 
reasoning by which Mr. Airy makes it appear that every pro- 


a 


tuberance outside this thin crust must be accompanied by a 
protuberance inside, down into the fluid mass, would equally 
prove that wherever there was a hollow, as in deep seas, in 
the outer surface, there must be one also in the inner sur- 
face of the crust corresponding to it; thus leading to a law 
of varying thickness which no process of cooling would have 
produced. 


Pratt (1859, p. 751) explained the deficiency of matter beneath the 
mountain mass in the following way: 


At the time when the earth had just ceased to be 
wholly fluid, the form must have been a perfect spheroid, 
with no mountains and valleys nor mountain hollows. As the 
crust formed, and grew continually thicker, contractions and 
expansions may have taken place in any of its parts, so as 
to depress and elevate the corresponding portions of the sur- 
face. If these changes took place chiefly in a vertical di- 
rection, then at any epoch a vertical line drawn down to a 
sufficient depth from any place in the surface will pass 
through a mass of matter which has remained the same in amount 
all through the changes. By the process of [thermal] expansion 
the mountains have been forced up, and the mass thus raised 
above the level has produced a corresponding attenuation of 
matter below. This attenuation is most likely very trifling, 
as it probably exists through a great depth. 


The deflection of the plumb line caused by the mass of the Himalayas 
and the mountain region beyond was calculated by Pratt as 27.978 see at 
Kaliana (at the foot of the Himalayas) and the difference of the deflec- 
tion of the plumb lines at Kaliana and Kalianpur as 15.941 sec. He then 
computed the deflection as modified by the supposed attenuation below the 
mountains assuming that the level of attenuation extended down to a depth 
of 100, 300, 500 and 1000 miles. His figures showed that the deflection 
at Kaliana, if the attenuation were extended down to the 100 mile level 


would be 1.538 sec and the difference between Kaliana and Kalianpur would 


22 


be 1.602 sec. However, the measured difference of deflection at Kaliana 
and Kalianpur was 5.236 sec which would require a depth of attenuation 
of almost 4300 miles. The problems posed by these results led Pratt 
(1859, p. 762) to declare, "... No hypothesis of deficiency of matter 
below, which we a habe ive will remove the anomaly. ‘The disturbing 
cause must be elsewhere; 4 

Despite the fact that all of Pratt's concepts of the nature of cae 
earth's crust (except that rock densities vary horizontally) have been 
disproved, his theory of isostasy (in greatly modified form) is stake 
being used to reduce gravity anomalies and, indirectly, to compute the 
Size and shape of the earth. 

The concepts contained in the hypotheses of Airy and Pratt are the 
basis for several isostatic systems jin use at the present time; the prine 
cipal ones being the Pratt-Hayford system, the Airy-Heiskanen system and 
the Vening Meinesz system. 

Figure 1 presents a schematic representation of the structure ofp mie 
earth's crust based on the concepts of: (a) Airy, (b) Pratt (modified), 
and (c) a composite view—probably the closest approach to actual con- 


ditions in the erust. 
ISOS TATIG. RECOVERY 


Glacial Loading 
The Hudson Bay-Great Lakes area was covered by a continental glacier 


of Wisconsin age which was probably more than 3,000 meters (10,000 ft) 


ee 


Ga). Airy CRUST 
hypothesis 

(b) Pratt 
(modified ) 
hypothesis 


ye: CRUST ome USiOl ieee: 
(c) Combined 
hypothesis 


Fig. 1. Schematic diagrams of the earth's crust as implied by the 
theories of (a) Airy, (b) Pratt (modified), and (c) combined hypoth- 
esis. Rock density increases with increased density of stippling. 
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thick if an analogy may be made with the neeteur ice sheets of Antarctica 
and Greenland. The maximum thickness of the Antarctic ice sheet is 4270 
meters (14,000 ft), and that of Greenland is 4,400 meters (10,800 ft) 
(Bentley and Ostenso, 1961, p. 886; Holtzscherer and Robin, 1954, p. 196). 
The shape of the continental glacier was that of a flat-topped shield 
with relatively steep gradients along the margins. 

The weight of this tremendous mass of ice piled on the earth's 
crust created a stress which exceeded the strength of the lithosphere, 
and, according, to the concept of isostasy, the lithosphere sank intowea. 
weak layer below (the asthenosphere) until a condition of approximate 
equilibrium was again achieved. As the ice accumulated the increasing 
Weight first. caused: the dithosphere to be derormed clastical ain. 
basin-like shape; then, as the stress exceeded the strength of the crugs 
and the asthenosphere, plastic flow occurred in the weak subcrustal layer. 

The movement of plastic material from beneath the glacial area to- 
ward the ocean basins, the floors of which were rising due to the de- 
creasing load caused by the progressive loss of water, was accompanied, 
according to many geologists (e.g., Jamieson, 1882, p. 461; Barrell, 
1915, p. 13; Nansen, 1927, pp. 40, 42-44; Burger and Collette, 1958, pp. 
227-240), by an upward bulging of the crust Saeedicle the area of depres- 
Slon.-) The tcormatronsoLeuic Roles eae by Nansen (1927, pp. 39- 
40) as follows: 


It is, however, obvious that when, within a certain area, 
the crust is depressed by deposition of sediments ‘ice] on its 


surface, this will cause an: outward flow of mobile matter in 


aD 


the plastic substratum; but owing to the resistance caused by 

the great internal friction, this flow cannot extend over in- 

definite areas, but will at first be limited to the surround- 

ing regions, where the surface will be somewhat elevated in a 

belt around the depressed area. The elevation will be highest 

near the latter, and will gradually decrease outward toward 

Zerg. 

The existence of a noticeable peripheral bulge has been denied by 
other workers who pointed out that the high bending strength of the lith- 
osphere would preserve a bulge to the present time, and that a permanent 
tilt should be recorded in the postglacial beaches south of their "hinge 
lines"—neither condition has been found. However, the apparent lack of 
a bulge may be explained if the peripheral bulge is so low and broad that 
it cannot be detected by present methods of field measurement. 

An estimate of the maximum extent and magnitude of the depression of 
the land due to glacial loading may be gained from an examination of the 


limits of the Wisconsin glaciation as shown in Plate I and from a con- 


Sideration of the simple formula 


dy 
Doe! 1 
a ey, 


where D = depth of depression, t = thickness of ice, d; = density of ice, 
and dg = density of the asthenosphere. 

The Wisconsin continental glacier extended from its center in the 
latitude of Hudson Bay to the southern boundary of the Great Lakes Lobe 
which terminated in the southern half of the States of Illinois, Indiana 


and Ohio. It must be realized, however, that the limits of the crustal 
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depression did not coincide with the southern limits of glaciation] sige 
depression was tare Hy the thickness of ice necessary to overcome the 
strength of''the underlying crust. Near the edge of the ice sheet theme 
was probably not thick enough to provide the necessary weight for depres- 
sion; whereas, farther in from the edge, the ice was thick enough formes 
weight to cause elastic deformation, and still farther north the weight 
of the ice exceeded the strength of the lithosphere and asthenosphere so 
that plastic as well as elastic deformation occurred. With suificiers: 
time the lithosphere under the glacier would be depressed until hydro-= 
static equilibrium was established. 


At the center of the ice sheet, an approximate depth of depression 


d- 
may be calculated from the formula D = tr, if the following assumptions 
a 


are made: (a) the ice thickness was ec. 10,000 ft, (b) the ice density 
was c 0.92, (c) the density of the athenosphere is approximately 4.3, 

and (d) the isostatic adjustment was complete. With these qualifications, 
the depth of depression at the glacial center would be on the order of: 
2,800 ft, decreasing in magnitude southward as the ice thinned. The 
profile of depression was. probably, the inverse protile or since olagier 


surface. 


Glacial Unloading 
The process of glacial unloading began after the regimen of the 
Wisconsin glacier had changed from a positive to a negative balance. 


The retreat of the glacier'’s terminus, .as well as the thinning of the 
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ice sheet by melting of its upper surface, brought about the reversal of 
those processes described in "Glacial Loading." The first reaction of 
the depressed crust to the decreased weight of ice was an elastic recoil 
of the crust; this, in turn, was followed by uplift due to plastic flow 
in the substratum. Plastic deformation occurred when the melt waters of 
the glacier returned to the ocean, and as the load on the land decreased, 
that on the ocean bottom became greater and greater—once again disturb- 
ing the condition of hydrostatic equilibrium—which caused a reversal of 
flow of material in the asthenosphere. 

The ice front retreated about 140 miles north of its southernmost 
extension (crossing the southern boundary of the Great Lakes watershed) 
before glacial margin lakes began to form between the southern moraines 
and the ice front. The glacier continued to retreat in a series of par- 
tial retreats and advances to the time of the Cary-Port Huron interval 
when, according to J. L. Hough (1958, p. 287), the ice front in the 
Great Lakes region probably extended from the vicinity of Oswego, New 
York, to the neighborhood of Beaver Bay, Wisconsin, on Lake Superior; 
passing near Gravenhurst, Ontario, the northern tip of the Saugeen Pen- 
insula, the Straits of Mackinac, Escanaba, Michigan, and the southern 
end of the Keweenaw Peninsula. 

During the time that the glacier was receding c 440 miles in a per- 
iod of 2000-3000 years (Hough, 1958, p. 278) elastic and plastic defor- 


mation were acting to restore the southern glaciated areas to their pre- 
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glacial elevations. The amount of land uplift which occurred during this 
period may answer the point brought up by Sherman Moore (1948, p. 708) 
when he declared: "... The ice of the Wisconsin glaciation extended al- 
most to the Ohio River, yet the old beaches show no differential move- 
ment below the hinge lines which cut across the Great Lakes." 

The old beaches failed to show differential movement because the 
land south of the “hinge lines" had almost completely resumed its pre- 
glacial elevation before beaches were formed. This serves to emphasize 
the fact that the amount of land uplift revealed by differential warp 
ing of glacial lake strandlines since ice retreat is a minimum value 
and represents an unknown proportion of the total recovery from glacial 
loading. 

The process of glacial lake beach formation has been used to explain 
two different concepts of the continuity of land uplift. Proponents ten 
both views agree that the strong beaches were formed either during per- 
iods of relatively constant relationships between the water level and 
the land, or during periods of slightly rising water levels. However, 
opinions differ as to the way in which the constant relationships were 
maintained. 

The more widely-held opinion is that the relationships were constant 
and shoreline features were formed during: times when land uplift had 
stopped so that stability prevailed. When land uplift resumed, the shore- 
line features were warped in those areas where isostatic adjustment was 


not complete. 
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The second viewpoint, that shoreline features were formed during 
periods: of continuous uplift when the rise of water level kept pace with 
the uplift of the land, has been summarized by the Norwegian scientist 
and explorer, F. Nansen (1927, pp. 47-48), as follows: 


The now raised post-Glacial strand-lines, beaches, and 
terraces of the formerly glaciated regions, have obviously 
been formed during periods of nearly constant relation be- 
tween sea-level and land, or when the sea was rising slightly 
relatively to the land. In an earlier paper [1922] the author 
has tried to prove that such conditions were not created by 
any break in the upheaval, still less by any new depression 
of the land in late-Glacial or post-Glacial time, as assumed 
by most previous authors. Owing to the lag in the isostatic 
adjustment there must obviously have been so great an excess 
of "buoyancy" in the crust during the period of its upheaval 
that any temporary pause in the melting of the ice-cap, or 
even a temporary increase of it, cannot have stopped the up- 
lift of the land, and still less have produced a temporary 
subsidence. 

The conditions for the incision of strand-lines, or the 
formation of marked beaches or terraces, have existed during 
periods when the sea-level rose at the same rate as the coast 
or slightly faster, and the coast-line remained more or less 
stable sufficiently long. 


The progressive return of the earth's crust to approximately its 
pre-Wisconsin elevation is manifested by the northeastward shift of the 
areas of horizontality of the strong beaches (Whittlesey, Warren, Iroquois, 
Algonquin and Nipissing) in the Great Lakes region. Plate I shows the 
location of the zero isobase of each of the five strong beaches, the 
Whittlesey beach is the oldest (c 12,700 yrs B.P.). Each beach was 
horizontal when it was formed; Ppeedient upwarping northeast of the zero 
isobase indicates that uplift continued here long after it had ceased 


southwest of the zero isobase. 
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According to one concept, after upwarping of the older shoreline 
northeast of the zero isobase occurred, crustal movement must have 
stopped in order that the shoreline features of the succeeding lake be - 
formed; if Nansen's coneept is followed, thevrates of uplift and visemes 
the lake's water level must have been synchronous in order for shorelime 
features to form. After the formation of the shoreline fesvures oF the 
younger lake, upwarping either began again or proceeded faster than the 
rise in water level, depending on which concept is followed, so thavuene 
shoreline features of the younger lake were warped. This process oc= 
curred for each of the *elacial take shorelines. 

The last distinct postglacial strand line to show upwarping is the 
Nipissing, which ended about 3,200 years before the present time (Farrand, 
1960, p. 125, Table III). ‘The area of horizontality of. the Nipissing 
beach outlines that portion of the Great Lakes area where pre-glacial 
elevations have been restored to approximately their former altitudes: 
The degree of restoration is probably not complete due to the thick over= 
burden of glacial drift which covers the glaciated region south of the 
Canadian Shield. 

The drift, composed of material carried to the area by the four 
continental glaciers, has an average thickness of about 400 feet in 
southern Michigan and northern Indiana (Leverett and Taylor, 1915, p. 61) 
and a maximum thickness of more than 700 feet in Michigan (Aker, 1938). 


The thickness of the drift deposited in the Great Lakes basins should 
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also be hundreds of feet. Although the portion of the total drift thick- 
ness contributed by the Wisconsin glacier is unknown, it was certainly 
appreciable and its weight would probably keep the bedrock from attaining 
its full pre-Wisconsin height. 

The above situation in the Great Lakes region may be contrasted 
with the one which exists in the Canadian Shield—the source area of the 
drift. Material from the Shield was removed by the ice and deposited 
elsewhere, making lighter the load on the underlying bedrock. Because 
the weight which formerly existed in the central area of glaciation is 
now distributed either in the outer areas of glaciation or carried away 
to the oceans, isostatic adjustment, when complete, should allow the 
Canadian Shield bedrock surfaces to become higher in elevation than they 


were before glaciation. 


III. DETECTION AND MAGNITUDE OF UPLIFT 


Land uplift in Fennoscandia and the Great Lakes region has been 
measured quantitatively by three methods: (a) by measuring the warping 
of former shorelines, (b) by comparing water-level records over a long | 
period of time, and (c) by precise leveling. In addition, gravity meas- 
urements and earthquake intensities have been used in a Qualitative oe 
to both affirm and deny the existence of residual depression and conse- 


quent modern uplift in these areas. 


Warped Glacial Lake Shoreline Features 


The most positive method of demonstrating the existence of postgla- 
clal uplift is by measuring the warping of the depositional and erosional 
land forms produced by the action of the former lake surface. Due to the 
length of time in which the upwarping action was able to act, the upwarp 
of the shoreline features amounts to tens or hundreds of feet as they are 
traced from south to north. 

The shoreline. features used to determine the amount of uplift in-= 
clude: beaches and bars (deposits of sand or pebbles accumulated by on- 
shore and along-shore waves and currents), deltas, wave-cut and wave- 
built benches, wave-cut cliffs and outlet. channels. These glacial lake 


features are now found at a higher elevation and farther inland than 


present lake features. 
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In order to cover the greatest area possible, the majority of the 
investigators who mapped the shoreline features of former lakes used al- 
timeters or hand levels to determine their elevations. The inaccuracies 
inherent in the instruments and methods coupled with the difficulties of 
determining the elevation of the former water surface—i.e., of correla- 
ting wave-cut features with wave-built features of lakes whose surfaces 
varied several feet in elevation from year to year—result in elevations 
which are probably accurate from as Sy G0 My 10 feet. Although some eleva- 
tions were determined by spirit level or by hand level from near-by known 
elevations, the uncertainities of measuring features formed by a fluctu- 
ating water surface remain; thus making an accuracy of better than + 5 
feet rather unlikely (Robinson, 1908, pp. 348-458). 

Gilbert's (1890, p. 368) method of representing a warped plane was 
followed by G. De Geer (1892, p. 457), the Swedish glacial geologist, 
who connected points of equal deformation with a line which he called an 
"isobase." The isobase of zero deformation which marks the boundary be- 
tween the area of horizontality of the glacial lake features and its area 
of warping was called the "hinge line" by J. W. Goldthwait (1908, p. 473) 
and subsequent workers. However, the term "hinge line" connotes a def- 
inite demarcation between the horizontal and warped areas which does not 
actually exist; as Farrand (1960, p. 9) has recently pointed out, "the 
apparent hinging effect shown in most profiles of former water planes 


(... ) is ‘produced by the exaggeration of vertical scale which is neces- 
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sary in such diagrams." Farrand (1960, p. 90) later called attention to 
the fact that the "isobase of zero uplift for the Nipissing beach ~ = eu 
difficult to determine within 10 or 20 miles because of the extremely low 
slope of the water plane as it approaches horizonality /1:135,900]." 

The difficulty of locating the position of the intersection ofeGa. 
warped and horizontal portions of glacial lake shorelines having gradi- 
ents ranging from c 1:1,400 to 1:14,000 may justify substituting the 
term "transition zone" for "hinge line." 

This problem of illustrating on a vertical profile the change from 
horizontal to upwarped lake shoreline features is one facet of the over= 
all problem of representing in two dimensions a phenomenon (occurring on 
an ellipsoidal earth) which has extremely large horizontal dimensions 
when these dimensions are compared to its vertical extent. 

The extreme shallowness of the depression caused by the continental 
glacier at its maximum may be illustrated by comparing the southward ex- 
tent of the ice sheet of 1,300 miles with the possible depression of the 
earth's crust of 2,800 feet at the center of ‘the glacier. ‘The average 
gradient in this case would be c 1:2,450 or the slope angle would be c 
1 mtin.24 sec. 

Distortions induced by the grossly exaggerated vertical profile neces- 
sary to show the relationships of horizontal and uplifted areas of shore- 
line features have led to the use of the term "doming" or "updoming" to 
describe the action of the land recovering from depression. An examina- 


tion of Fig. 2 reveals that the "doming" is illusory, caused by the up- 
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(a) Land sky line profile before glaciation. 
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Fig. 2. Schematic diagrams showing recovery of the depressed regions 
and upwarping of shoreline features. 
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warping of sections of former shorelines. The crust which was depressed 
under the weight of the ice rises only until it assumes its preglacial ele- 
vation; the elevation of the land at the center of the ice sheet may be 
higher after isostatic recovery than it was before glaciation (see p. 31). 

The glacial lakes were formed in topographical low areas of the 
regions depressed by the weight of the ice sheet; usually (as shown by 
the attitudes of former beaches) the lakes extended from non-depressed 
regions to the depressed peers in front of the glacier. The shore. 
lines reached from an area of horizontality, across an area free of mee 
but still depressed due to the lag in isostatic recovery, to the glacier 
front. After ia period of time in which the glacier retreated to 7am 
position, another lake formed at a lower elevation than the first; iso- 
static recovery occurred and the depressed area rose to approximately 
its former elevation; carrying with it “the former shoreline.” ius soa. 
former depressed areas are now essentially horizontal, although part of 
the former shorelines are now upwarped. 

The positions of the isobases of zero deformation shown on Plate I 
demonstrates the fact that the zero isobases are almost parallel with each 
other and that the locations of the zero isobases show successively 
younger isobases to be northeast of the older ones. 

Isobases found on maps depicting postglacial uplift of the land 
(Hough, 19598, Migs.) 34) 42, 51; Flint, 19>7jj0tics. eh Oe eee 


6, 14-9; Daly, 1954, Fig. 64; Sauramo, 1939, Pigs 1, ete.) are resionel 


if 


isobases, and, by smoothing out irregularities, show the broad gently 
curving trends of deformation,’ However, detailed studies of the local 
trends of the isobases have revealed that the amount of uplift, as well 
as the direction of the trend, was influenced by geological and topo- 
graphical features (DeGeer, 1892, p. 458; Leverett and MacLachlan, 1944, 
p. 550; Sauramo, 1939, pp. 15, 21-23; MacLachlan, 1939, pp. 63, 80-82). 

MacLachlan (1939, pp. 80-81) studied the Glacial Lake Warren shore- 
line and found that: 

The water plane of Lake Warren was not only deformed 

by northeastward tilting as indicated by straight line iso- 

bases connecting widely separated points, but there was also 

local deformation of the area along these isobases of regional 

tilt and this local deformation reflects the major details of 

buried structures. 

In addition, the spacing of the regional isobases indicates that the 
former shorelines were not tilted as flat planes, but were warped upward 
into gentle curves. The progressively closer spacing of the isobases to 
the northeast shows that the concave profile of the shorelines became 
steeper to the northeast, and, consequently, that the rates of uplift 
increased in that direction. 

Average rates of uplift for six shorelines in the Great Lakes region 
(Whittlesey, Warren, Grassmere, Algonquin, Iroquois and Nipissing) are 
given in Table 1. The rates were found by measuring from the zero iso- 
base to the highest definite isobase of the shoreline being measured. 


The average rates do not illustrate the progressive steepening of the 


shoreline profiles to the northeast, e.g., the average rate of uplift for 
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the Nipissing shoreline in the Superior basin is 0.49 ft per mi; whereas, 
Farrand (1960, p. 56) finds that: "It is horizontal south of Knife River, 
Minnesota, but slopes 0.38 ft per mi in the Tofte-Lutsen area, 0.42 ft per 
MenieareGneand, Frorrage Minnesota, and O.50 ft per mi in» the Nipigon, 
Ontario, area." 

Rates of uplift computed from the elevations of former shoreline 
features are based on quantities which could be the result of up to 
13,000 years of land uplift; furthermore, the differences in elevation 
between the zero isobases and the isobases of maximum deformation are 
measured in tens or hundredsof feet. The magnitude of these quantities 
is great enough for elevations determined by ordinary methods of spirit 
leveling (corrected for meteorological effects and for the fluctuation 
of the modern lake datum) to be used in the calculation of accurate 
rates of past uplift. 

This situation contrasts greatly with the determination of modern 
rates of upliit by means of precise leveling; or by means of water-level 
gage records; where the period of record is measured in tens of years and 
the magnitude of uplift is in hundredths of a foot, in tenths of a foot, 
and, in some cases, in feet. Under these circumstances, errors in the 
determination of elevations which would be insignificant in the calcu- 
lation of rates of past uplift constitute such a large proportion of the 
total amount of modern rates of uplift that the modern figures could be 


xf 
valueless. \ 


LO 
Water-level Gage Comparisons 


Modern land uplift values have been calculated from comparisons of 
water-level gage records (tide gage records and lake-level gage records ) 
taken over a period of 50-100 years. ‘The principle forming the basis for 
this type of computation is that of "water leveling;" which, as defined 
by the’ U. 8. Coast and Geodetic Surveys ls: 

A method of obtaining relative elevations by observing 
heights with respect to the surface of a body of still water. 
The surface of 4a body ofsstiss) waver, a6 of lake, eis 

a level surface (equipotential surface), and the relative 

elevations of objects along its shores may be obtained by 

taking the differences of their heights with respect to the 

surface of the water. ... (Mitchell, 1948, p. 46). 

The comparison of heights taken from tide gage records furnished 
the necessary data for the computation of rates of uplift around the 
Gulf of Bothnia and the Gulf of Finland in Fennoscandia, and differences 
in elevation of lake-level gage stations have supplied information nec- 
essary for the calculation of rates of uplift in the Great Lakes region 
and the lake district of Finland. 

Tide: gage records. are used not only to provide a direcuameanc: ou 
calculating rates of uplift, but are also used indirectly when uplitc 
rates are determined by precise leveling. Mean sea level (MSL), the 
standard datum for elevations, is established by means of tide observa- 
tions taken over a number of years; in addition to providing the standard 


datum, primary tide gage stations are used as starting and "tie-in" 


points 
in..the precise leveling net (Marmer, 1951, p. 67;.Hayford; 1922, pp. 131- 


152). 


Ka 


Great Lakes rates of land uplift are also influenced by this indi- 
rect function of tide gages determining mean sea level if absolute rates 
of movement are desired—lake-level gage records alone allow only relative 
rates of uplift to be found. Tide gage stations in the New York City area 
(Sandy Hook, New Jersey; Governors Island, New York; Fort Hamilton, New 
York; and the Battery, New York City) have been the starting points for 
level lines used to establish elevations in the Great Lakes region; there- 
fore tidal records of the New York City area are important in the deter- 


mination of Great Lakes absolute rates of uplift. 


WATER-LEVEL GAGES 

The component parts of the water-level gage station illustrated in 
Fig. 3 are common to almost all permanent installations. A recording 
device is mounted on a stable platform situated over a stilling well 
which is connected to open water by a pipe of relatively small diameter. 
The stilling well with its restricted inlet eliminates short period fluc- 
tuations without affecting the recording of the water level. ‘The re- 
corder consists of two basic parts, one is a clock mechanism to move a 
roll of paper under a stylus at a uniform speed, and the other is a 
height registering mechanism with a float which, through a linkage of 
float-line and gear train, moves the stylus a distance which is propor- 
tional to the change in water level. 

An index (zero) mark on the stable platform is connected 


by a double line of spirit leveling to three or more permanent bench 


he 
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marks. At regular intervals the leveling from bench marks to gage is 
repeated in order to detect any changes in the zero of the gage. Peri- 
odically, an observer, using a staff gage or electric contact tape, 
measures the distance from the index mark to the surface of the water; 
he records this figure and the time so that during tabulation of the 
record the automatic gage record may be correlated with the elevation of 
the index mark. 

After the automatic gage roll has been removed from the recorder, 
hourly elevations are tabulated from the continuous neeoe the recorder 
roll in order to provide a basis for computing daily, monthly and yearly 
mean water-level elevations. 

Water leveling, the determination of mean sea level, and the cal- 
culation of rates of land uplift from water-level records are all based 
on the tacit assumption that the mean elevation of the water surface as 
determined by the record from the gage station represents the actual 
mean elevation of the body of water. This assumption would be very nearly 
true if the surface of the water body were a level surface (equipotential 
surface), and if the tabulated values of the water elevation were repre- 
sentative samples of the level surface elevations. However, the water 
SUriace ie rarciy, Jf eyer, 2 level surface and its elevation, as re- 
vealed by water-level gage records, is affected by a number of factors 
whose influence must be removed orrendered insignificant if the gage 
readings are to an accurate representation of the actual elevation of 


the water surface. 
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Some of the factors which disturb the level of the surface of “a ™poq. 
of water or distort its true elevation, i.é., meteorological effects; sep. 
erator and instrument errors, benchmark and/or gage movements, are common 
to both sea and lake surfaces; whereas, other influences, i.é., the ‘eusrae 
tic rise of sea level, density changes and tides, are important only in 


computing elevations of the sea surface. 


Tide Gage Comparisons 

The principle underlying the calculation of rates of land uplifeee, 
means of tide gage comparisons is a simple one—-the elevations of certain 
fixed points along a coastline (tide gage index marks) are compared with 
a plane of zero elevation (mean sea level) over a long: period of time: 
If the differences between the zero plane (MSL) and the elevations of 
the tide gage marks remain constant, no uplift ae taken place; on the 
other hand, if the differences increase with respect to the zero plane, 
uplift has occurred; and the rate of uplift is found by means of the 
equation for a straight line 

5 A ee eI Co) Sle! (2) 

where y = the differences in the elevation of the tide gage; x = the 
time, in wears; )b== uplift ab thestars of the time series; and m = the 
angle coefficient, the value of the rate of yearly land uplift. 

The principle is’ a simple one, ‘but difficulty arises in practice 
due to the problem of establishing the zero plane; that is, the true 


mean sea level. Because the factors, e.g., meteorological effects, 


45 


tides, operator and instrument errors, influencing the determination of 
mean sea level vary their effects from one time to another, as well as 
from one geographical location to another, they must be discovered, 
their magnitudes computed, and their effects removed from the recorded 
data before an accurate mean can be found. 

Requirements as to the degree of accuracy which is necessary for 
establishing mean sea level varies for different purposes—due to the 
very small magnitude of modern land uplift, about 0.008 ft/100 mi/yr 
(see page 91), the accuracy requirements which are necessary when deal- 
ing with this phenomenon are stringent; thus the effects of the disturb- 
ing elements must be removed until their residual influence is very 


small. 


Tides 

One of the factors producing the greatest deviation of the sea's 
surface from that of a theoretical equipotential surface is the tide, 
which is defined as: "The periodic rising and falling of the water that 
results from the gravitational attraction of the moon and sun acting upon 
the rotating earth." (Schureman, 1949, p. 36). The equipotential sur- 
face necessary for use as a datum plane is approached when the effects 
of the tide and meteorological factors are eliminated by using tide gage 
records tc determine the mean sea level. 

Because tide-producing forces vary on a daily, monthly and yearly 


basis, it is necessary that tidal records be taken over a number of 


LO 


years before mean sea level can be established. H. A. Marmer COSA s PRs 
64-64) of the U. S. Coast and Geodetic Survey states that: 


A period of 19 years is generally considered as consti- 
tubing a full tidal ‘cycle, for during) this period of (Lime tae 
more important of the tidal variations will have gone through 
complete cycles. It is therefore customary to regard results 
derived from 19 years of tide observations as constituting mean 
values. 

If the mean level of the sea remained constant over long 
periods of time and if the coast were absolutely stable, we 
might expect sea level at any place determined from one 19- 
year series to be the same as that derived from another such 
series even if separated by, aonumber ofiyears. wAppaveuiily, 
however, this is. net. the casey ‘and form precise= pur posecmin 
isvtherefore necessary toespecil y the particular epociwusced 
invthe determination of mean sea level, 

For New York Harbor there are available 56 years of obser- 
vations, from 1893 through 1948. This permits three 19-year 
series, 1894-1911, 1912-1930 and 1930-1948, the last two hav- 
ine the year 1940 in ceomions | For the series: 192-7 Omeoce 
evel referred to a number of bench marks In thes vicinity. 
the tide station was 0.09 foot higher than for the series of 
1693-1911; for 1930-1948 it was 0.29 foot higher than for 
1893-1911, and 0.20 foot higher than for 1912-1930. 


It would appear from the above quotation that mean sea level derived 
from periods of 19 or more years would eliminate almost all of the dis- 
turbing influences produced by the tides and meteorological effects; 
however, the gradual rise of mean sea level at tide stations outside 
areas of known crustal movements points to an eustatic change in sea 
level; and, considering local areas, to residual meteorological effects 
of long period climatic changes: 

A consideration of long period variations in mean sea level is im- 
portant in the calculation of rates of land uplift by precise leveling 


where re-leveling is carried out 40-80 or more years after the initial 
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precise leveling. The long period rise in sea level may have two effects— 
both causing a relative subsidence of the land—if compensation for the 
rise is not carried out. In the case of tide gages and accompanying 
bench marks which are located in an area of land uplift, such as Fenno- 
scandia, the failure to correct for long period changes in mean sea level 
will cause the rate of land uplift to appear to be too small because the 
zero reference plane is rising at the same time that the land is being 
uplifted. The rate of rise of mean sea level must be added to the rate 
of land uplift if the true rate of crustal movement is to be determined. 

The second effect of apparent subsidence occurs when a series of 
bench marks are re-leveled after a number of years, and the calculation 
of elevations is made on the assumption that mean sea level remains con- 
stant over long periods of time. In this case, the zero datum is assumed 
to have the same relation to the tide gage bench mark for the second lev- 
eling as it had for the first leveling; whereas, it is actually closer to 
the bench mark. Because the zero datum plane is held constant, it ap- 
pears that the bench marks are subsiding; although, no actual crustal 
movement has occurred. This effect of apparent subsidence may be illus- 
trated by an example based on information contained in the quotation 
taken from H. A. Marmer (see page 46). 

Mean sea level in New York Harbor rose 0.29 foot from 19@ (mid-point 
1893-1911 series) to 1939 (mid-point 1930-1948 series)—if a line of pre- 


cise levels were run in 19@ and again in 1939, and if it were assumed 
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that no change had taken place in the elevation of mean sea level; then 
the bench marks would appear to have subsided 0.29 foot in 4/ years, or 


at a rate of about 0.76 foot per 100 yeans. 


Meteorological Effects 

It has. long been known that.the: level of the sea is affected not 
only by astronomical tides, but also by the force of the wind which blows 
across the surface, causing a piling up of water on the lee shore andea 
lowering of water level on the windward shore. Although this effect of 
the wind has been known since ancient times, it was not until 1804 that 
the influence of barometric pressure on the level of the sea was demon= 
strated. In 1804, a Swedish physician named Schulten correlated changes 
in barometric pressure with variations in the height of water level in 
the Baltic Sea and discovered that at a given locality the greatest in- 
crease in height of the water corresponded to the greatest depression of 
the mercury column of the barometer. He found that the maximum varia- 
GON of 2.5 inches of the mercury column was the equivalent of a change 
invsea level of 34 inches, which corresponds very closely to the theo- 
retical ratio of 1:13.21. Schulten's measurements together with his ob- 
servation that the rise in sea level always preceded the drop in the mer- 
cury column, led him to conclude that the change in level was due "to the 
unequal pressure of the atmosphere upon different parts of the surface; 


i O(Mi Cul hoehyerolisen par eOo)). 
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Despite the fact that the term "meteorological effects" could be 
applied to changes in water level brought about by local additions of 
water to the sea by rain and river discharge, as well as by density dif- 
ferences caused by dilution of the salinity of the sea and temperature 
changes, it usually includes only the effects of wind and barometric 


pressure (see Fig. 4a and 4b). 


Wind Slope 

As wind blows over a water surface an interfacial stress is pro- 
duced by the viscous drag of the moving air on the water and by the form 
drag of the waves and wavelets (Montgomery, 1952, p. 132). This stress 
causes the water to move in a general downwind direction until it reaches 
shallow water or the shore where a piling up of the water creates a slope 
OF the water surface. The surface slope,in turn, creates a gradient or 
gravity current which flows beneath, and in opposite direction to, the 
Wwind-drift current. 

Various authors have labeled this wind-induced surface slope (after 


corrections have been made for the barometric effect) a "wind slope," a 


! 1 


"wind denivellation," a "wind tide," a "wind effect," a "windstau," and 


wind "set-up." The term "wind slope" best describes this phenomenon; as 


"wind denivellation," 


while the most precise term is probably too long 
to be commonly used; "wind tide" has the connotation of periodic rising 


and falling due to astronomical influences; "wind effect" is not definite 


enough; "windstau" when translated as an "accumulation" or "banking up” 
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by the wind has the proper meaning but is not as precise as "wind slope;" 
and the engineering term "set-up" is best used as 
either the displacement of the water surface at the lee- 
ward end of the channel with respect to the undisturbed level 
of the water or the difference of the displacement of the 


water levels at the windward and leeward ends of the channel, 
that at the windward end being negative (Keulegan, 1951, p. 


365). 

Two equations taken from G. H. Keulegan (1951, pp. 360, 576) illus- 
trate the relationship which exists between the wind and a body of water 
when a slope is created through wind action. The theoretical develop- 
ment of the equations for determining the interfacial stress and wind 
slope assumes that: 

(a) the air and water are homogeneous, 

(b) the eddy viscosity is independent of depth, 

(c) the stability of the air is nearly "indifferent," and 

(d) an "equilibrium" state of the sea exists. 

(Sverdrup, et al., 1942, pp. 473, 476, 491, 492, 495, 498; Montgomery, 
Popeye pe 155; Huteningon. 1957 ,—0pp. 251; 205, 266, 26/.) 

The set-up, the difference in elevation of the water surface between 

the. windward and leeward ends of a body of water, may be determined by 


the equation 


where 
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and where 
S = total set-up 


A Rar To 
n = numerical coefficient based on cae a 1, where To = bottom stress 
S 


and T, = surface stress 
Tg = tangential surface stress 
L = distance between the maximum and minimum water elevations 
o = density of the liquid 
go =secceleration due to graviuy 
H = undisturbed depth of water 
of density of moist air 
K, = drag or resistance coefficient, also written Cq and xe 


yes component of wind causing the maximum set-up 


An accurate determination of set-up from the quantities involved in 
equations (3) and (4) is difficult due to the wide variances existing in 
the values for -n, K.and the exponent or Ve. 

The quantity n (also labeled x, H and Cal)s the ratio of the bottom 
stress to the surface stress, eS from 1 to. 1.5, depending uUponstue 
assumptions made by the investigator, and whether the flow is laminar 
(n= 1.5) Mormmpurbulent. (1 <n 15.) (Helis trom: ho Mis uepe ne mao ls 
Keulégan, 1951," pp. 36, 61; Van Dorn, 1955, pe 255.) Huschinson, 1957p. 
215). 

The extreme range of values for the drag coefficient K. is aptly 


illustrated by B. W. Wilson's (1960) review of the work of some 47 
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investigators dealing with wind-stress measurements. After "... adapt- 
ing the laboratory data for wind speeds at 10-cm height (usually) to 


" 


prototype conditions of wind speeds at 10-m height ... and adjusting 
the results of various workers to the "U- Law," Wilson (1960, p. 3381) 
found that: 
The over-all average value (for strong winds) is 2.37 x 10°° 
with a standard deviation of 0.56 x 10°, 
The average value of Cg for light winds is found to be 

1.49 x 10°73 with a standard deviation of 0.83 x 10°. ... The 

over-all range of values in this case extends from 0.4 x 10% 

Loeice: Sel Ome sandeeven O.2. Xx aloes though the latter value has 

been discounted in the averaging. 

It may be seen from the conclusions reached in Wilson's paper that 
the determination of the value of the drag (or resistance) coefficient 
has improved only slightly since 1952 when R. B. Montgomery (1952, p. 
134) declared "... the present evidence is so conflicting that at no 
wind speed is the resistance coefficient known confidently within half 
its value.” 


The third key quantity, the exponent of the wind velocity, as ex- 


pressed in the general formula of equation (4), that is 
Tg. = Kgpg(V-vo)™ (5) 


(where Vo = the velocity of the surface water, which is small as compared 
to V and is dropped) is usually taken as the square power; although 
values from 1.5 to 5 have been used (Wilson, 1960, p. 3378; Hutchinson, 


1957, pp. 273-275; Sverdrup, et al., 1942, p. 490). 
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The determination of the wind velocity term—usually considered as 
ve—is further complicated by the fact that the set-up (as measured at 
any point) is the result of two component winds, a regional wind and a 
local wind. Therefore, it is necessary to determine the velocity of two 
winds whose influence varies in relative importance in response to chang- 
ing velocities, durations and areal extent (Miller, 1958, p. 1; Gallé, 
1925, pp. 917-918). 

Theoretical determinations of the types of wind-driven currents 

in the ocean and their direction of motion with respect to the wind have 
been derived by V. W. Ekman and eee from the Euler-Navier theory of 
jaminar flow with the coefficient of viscosity in the equations being 
replaced by an "eddy" viscosity, a mechanical viscosity which depends 
upon the nature and state of motion of the fluid. The development of 
these theories also requires that the assumptions listed in a previous 
paragraph be made; i.e., that: (a) the air and water be homogeneous, 
(b) the eddy viscosity be independent of depth, (c) the stability of air 
and water be nearly indifferent, (d) equilibrium state of the sea exists, 
and in addition that, (e) the water be so deep that. no see stress: 
exists, and (f) that the horizontal stresses and their frictional forces 
not be considered (Sverdrup, et al., 1942, pp. 469-500; Hutchinson, 1957, 
PP. 259-2743 Hellstrom, 1941. pp, .28-51;. fekert, L060; ppd) O>— an 
95-97). 

Because conditions in nature seldom conform to the assumptions under- 


lying the theoretical basis for determining the types of currents and their 
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TeLiebous Loe NeaWwind,. becalice, laminar ftlowerarely, if ever, occurs. in 
large natural bodies of water, and because doubt has been cast on the 
validity of using the Euler-Navier theory with the replacement of u 
(dynamic coefficient of molecular viscosity) by € (coefficient of eddy 
viscosity) and y (kinematic coefficient of molecular viscosity) by K 
(kinematic coefficient of eddy viscosity) for turbulent flow (Dryden, 
et al., 1956, pp. 389-390), the numerical values derived by Ekman for 
the angles between winds and wind currents, for the depths of frictional 
Pesiouance, LOr the sizes of the bottom Ekman spiral, etc:; at various 
water depths do not have the universal application that other investi- 
gators have given them. 

The improbability of constructing an exact and complete mathematical 
model of the varying, interacting factors which exist in nature in the 
oceans and atmosphere demands a constant analysis and adjustment of 
theory when theoretical results conflict with observations of natural 
conditions. Very often empirical observations of the angles between 
wind direction and wind-drift currents in the oceans and in shallow 
bodies of water have been dismissed as being in error when they dis- 
agreed with Ekman's angle of 45°; in addition, the results of current 
measurements and directions of water slope have been questioned because 
EKkman's theory predicted very small or no angles between wind direction 
and current direction and water slope. Although Ekman's theory was the 


first and most important workable explanation of ocean currents, as with 
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any theory it must be evaluated in the light of new knowledge. 


Long Period Meteorological Changes 

Changes in wind velocities and barometric pressures throughout time 
are important as the existence of land uplift revealed by tide gage rec- 
ords taken over a number of decades may be obscured by a progressive 
change in mean sea level caused by long period changes in the atmospheric 
circulation. This type of change in mean sea level can be brought apege 
by a residual meteorological effect, as well as by eustatic change in 
sea level—neither of which is eliminated by averaging the tidal records 
over a-full tidal cycle of 19 years (Marmer, 1935, pp. 30, 33, 48; 19465 
pp. 201-02; 1951, pp. 63-64). At least two investigators (Kd#ridinen, 
LOD, et; Bergsten, 1950, sp. 52) have pointed out that presumed vari- 
ations, in the amount of land uplift, as revealed by tide gage records 
eorrected.for known meteorological effects and eustatice change, are 
actually mean sea level variations resulting from long term climatic 
fluctuations. 

The inerease in mean sea level. of 0.29: Toot from 1902 to 1959 12 
New York Harbor (or 0.78 foot per 100 years) may be also partially or 
wholly explained by long range temperature and wind field changes. 

Recent studies (Mather, 1954, pp. 287, 297; Baum and Haven, 1956, 
pp. 441-42, 447-48; Bjerknes, 1959, pp. 65-69; Landsberg, 1960, p. 1519) 
have shown that changes in the over-all circulation of the atmosphere, 


revealed by increased temperature and pressure changes, have occurred 


eal 


in the last fifty years. During this period, which corresponds approx- 
imately with the time that continuous tide records of New York Harbor 
have been kept, the temperature over the Northern Hemisphere has in- 
creased and the pressure differences in the North Atlantic off the east 
coast have also increased. 

The water supplied by the melting of the glaciers and the increased 
heating of ocean waters shown by the increased annual temperature may be 
the chief causes of a eustatic change of sea level amounting to about 11 
em/100 years (Gutenberg, 1941, pp. 729-30; Kuenen, 1950, pp. 532-35). A 
eustatic rise of sea level of this magnitude would increase the elevation 
of mean sea level in New York Harbor by 0.13 foot from 1902-1939. 

The height of mean sea level in New York Harbor is affected pri- 
marily by astronomic tides and by water piled up by barometric pressure 
Rr rrerenced and the wind. Another factor which may influence the height 
of mean sea level at this location is variation in the discharge of the 
Hudson River—the discharge averages 26,000 cubic feet of fresh water 
per second at the Narrows, according to H. A. Marmer (1935, p. 2). 

The barometric pressure effect and wind slope are composed of local 
and deep water components, both of which depend upon large scale atmos- 
pheric pressure differences. If mean sea level is derived from 19 years 
of tidal observations, the meteorological effects as well as the astro- 
nomical tides should be eliminated. However, if a progressive increase 


in pressure and geostrophic wind has taken place since the 1890's as 
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suggested by Bjerknes (1959, p.69),the residual effect of the increased 
piling up of water would be a rise in mean sea level. This increase 
could account for the 0.16 foot of mean sea level increase at New York 
remaining after the 0.14 foot change due to the eustatic rise of sea 


level had been accounted for. 


Bifeco of Vide Gage Location 

Primary tide stations which provide the necessary data for the com- 
putation of land movement, as well as for furnishing end points of pre- 
cise leveling lines, are usually located near coastal cities and towns 
which are almost always built along estuaries and harbors. ‘The factors 
that influence the elevation of mean sea level (tides, wind slopes and 
barometric pressure differences), have their greatest effect at those 
places where the water is shallow. In addition, the phenomenon of reso- 
nance which Proudman (1954, p. 200) calls: "the primary cause pie 
pronounced amplification (of the tide) observed in some gulfs and 


" 
° 


estuaries occurs in bays, harbors and estuaries. In other words, 

in many cases primary tide gages are situated in those areas where changes 
in the terrestrial and submarine topography will have the greatest effect 
(Harrie, (1907, pp. 455, 453, 455s Hayrord, eo pp. 124, 1243) Marner. 
19555. PP> a5, 45,0525 019515: pp. 25, 455 Sverdrup, satya, C2 sop eae 


554, 5623 Harris) 1954, p..5; Proudman,; 11954, psidi99,62 00,9202 se panauen 


1956, spp 60; iene 
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Land uplift, due to its small magnitude (e.g., 0.75 foot/100 
miles/100 years in the Baltic region), cannot be accurately measured 
unless index marks are compared with mean sea level over a period of 
several decades. If mean sea level (as determined by tide gages) changes 
during this period, the variation must be compensated for before land up- 
lift can be calculated. During the decades in which the records are 
being taken, natural and man-made changes occur in the coast line and 
in the submarine topography which affect the range of the tide (Marmer, 
1951, pp. 132-133). For example: the dredging of channels and dumping 
of spoil in shallow water; the filling in of shallow areas to extend 
the shore; the building of breakwaters, groins, piers, etc.; the growth 
of spits, bars and deltas; all change the configuration of the bottom 
and the shore, which in turn influencesthe height of the tide and of 
wind setup thus effecting the determination of mean sea level. 

Two additional factors which may influence the determination of 
mean sea level at tide gages which are located in estuaries and straits 
are: (a) variations in the discharge of the stream, and (b) the dif- 
ference in range of tide on opposite shores of a strait due to the de- 
flective force of the earth's rotation (Marmer, 1935, pp. 19, 21, 57, 
67-68; 1951, p. 45; Jensen and Sinding, 1945, pp. 16-19). Adjustments 
must be made for these factors when the records of two different tide 
gages are compared; or when records covering a number of decades from 


one tide gage are compared. 
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Instrument and Operator Error 

The use of modern recording gages, modern operating procedures, and 
modern methods of record reduction have reduced the limits of error in 
the determination of mean sea level to negligible proportions for ordi- 
nary purposes; however, the computation of the rates of land movement 
from water-level records must use still more accurate methods because of 
the very small magnitudes involved. An additional problem stems from the 
fact that calculation of rates of crustal movement requires that mean 
sea level figures derived from relatively accurate modern instruments 
be compared with mean sea level determinations of forty, fifty, or more 
years ago—a period when neither the instruments nor the procedures were 
as precise as those of today. In order to evaluate the accuracy of 
water-level records, it is necessary to examine the errors which may 
‘have their source in the recording instrument and its operators. 

Errors may enter in the determination of mean sea level in one or 
more of the following steps: (a) in the recording of relative sea level 
by the automatic gage, (b) in the recording of "absolute" sea level by 
the operator using a staff, hook, or electric tape gage, and (c) in the 
reduction of data from the continuous curve of the automatic water-level 
gage to daily, monthly, and yearly means. ‘The discussion of instrument 
and operator error will not include the effects of obvious gross errors 
resulting from equipment malfunctions (e.g., leaky floats, improperly 


calibrated staff gages, freezing of float in the stilling well, etc.), 
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or personnel blunders (e.g., failure to read and record numbers correctly, 
failure to keep stilling well inlet clear, failure to take staff gage 
readings, errors in computations, etc.). 

Water-level recorders, being float operated instruments, are sub- 
ject to.certain:- errors inherent in their construction.. The errors are 
usually very small in instruments of modern design with large floats 
and light lines and counterweights; they assume much greater propor- 
tions, however, in the more crude instruments of 50-100 years ago. 

As may be seen in Fig. 3 (p. 42), the float, through the medium of 
the float line and counterweight, performs work in operating the recorder 
pen. Owing to the friction in the instrument there is a lag between the 
force being applied by the water and the recording of the water level on 
the recorder chart. If the index is set while the float is rising, the 
rising stages will be recorded correctly, but the falling stages will be 
above the true level by the amount of the "float lag;" if the index is 
set for a falling stage, the rising stage will be low due to the lag. 

When the float rises, part of the float line passing over the pulley 
adds its weight to the counterweight; thus lifting the float and causing 
the index pen to register a height greater than the elevation of the 
water level. When the water level goes down, the weight of the increased 
length of line is added to the float, causing the pen to record an ele- 
vation lower than the actual water level. If the counterweight and the 


line enter the water, their weight will be reduced by buoyancy; the float 
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will sink deeper into the water and the recorded water level will be too 
low (Stevens, n.d., pp. 20-27; 1919-1920, pp. 394-395; Corbett, et al., 
1943, pp. 183-189). 

The magnitude of the above errors may be computed by the following 
empirical formulas from J. C. Stevens's Hydrographic Data Book (no date, 
pp. 27-28). Stevens's formulas, which give the errors in feet with ia 


struments using lead counterweights and steel float lines, are quoted 


below: 
; F 
Maximum error due to float lag = 0.37 peccrctee (7) (6) 
Error from submergence of 2 
counterpoise tS ORO LN mec (a) ra) 
Error from line shift with 
counterpoise in air Se ae 5B Ar a Cie) (8) 
Brror from line shift with 7 
counterpoise submerged = Oa 7 AG sees | (9) 
in which 
F = force (pull on float line) in ounces to move index 
D = diameter of float in inches 
u = weight of float ine in ounces’ per) foot 
@e = weight of counterpoise in ounces 


AH = change in stage in feet from previous setting. 

Float records may also be affected by changes in humidity which can 
cause an expansion or contraction of the recorder paper, or, in those 
instruments using hemp float line, which cause a change in the length of 
the line. If the index pen records on the gage paper in a humid environ- 
ment, where the paper is expanded, and measurements of the water stage are 


taken from the record in a normal environment, where the paper has resumed 
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its normal.dimensions, an error will be incorporated in the readings. 
The error can be compensated for only if two base line pens (a fixed 
distance apart) record in opposite margins of the chart. In this case 
the actual distance between the margin lines can be compared with the 
known fixed distance and the corrections computed (Stevens, n.d., pp. 
2-50; Chrystal, 1908, pp.. 361-370). 

The role of the operator in the recording of :relative sea level by 
the automatic gage has been summarized by J. C. Stevens (n.d., pp. 39-40) 
who said: 

The Human Equation constitutes the greatest possible 

SOurCe Of .erTror., | Incorrect. gage reading,. inaccurate setting 

of pens and pencils, failure to wind the clock, failure to 

start clock after being wound, failure to put stylus on the 

paper, wrongly attaching float so that pen moves up when it 

should go down, failure to release set screws or to tighten 

them as the case may be, failure to note gage reading and 

time on record sheet; failure to see that inlet to well is 

open; failure to ofl bearings occasionally, or the -use.of 

gummy oil, are among the few things that all too frequently 

result in erroneous or incomplete records. 

The records of relative heights of sea level are changed to actual 
sea level elevations by comparing elevations on the gage record with 
separate readings made simultaneously by the operator with a staff, hook 
or electric tape gage at periodic intervals. The separate height deter- 
mination measures the distance between an index mark on the water-level 
recorder and the surface of the sea. Because the index mark could, and 


often does, change due to accident, replacement of the equipment, etc., 


it is essential that the elevation of the index mark be tied by precise 
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spirit leveling to three or more permanent bench marks. It is obvious 
that for a study extending over a number of years, it is of utmost im- 
portance that regular connections be made between tide gage and bench 
marks in order to detect any change in the elevation of the index mark; 
probably one of the greatest errors in the determination of mean sea 
level over a long period of time is the failure to check the herghe ce 
the tide gage index with its bench marks, 
Other errors which may occur in the determination of "absolute" 

sea level by means of a staff gage are those due to the "error of the 


" 


interval" of the staff gage and the "personal error" of the observer. 


The “error of the interval" is, according to R. Gibbs (1929, p. 71), 
he the same as the probable error of the last figure; namely one-dquarcer 
of the interval." H. A. Marmer (1951, p. 30) in speaking of U. S. Coase 
and Geodetic Survey procedure, states that: "The observer reads the 
staff to the nearest half tenth of a foot if “the water is free from 
waves, or to the nearest tenth giving the highest and lowest readings." 

In this case the. error shold be £:0-01e.fdot orvtvO.025 foot. 

The "personal error" of the observer results from the particular 
way in which an observer: takes his observations; e.g., the observer may 
always assume a certain position, so that his eye is a little above the 
index mark; thus introducing a degree of parallax which would differ from 


that in the reading of another observer; or one observer may tend to 


mentally subdivide the interval between the graduations so as to favor 
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recording even tenths of a foot, whereas another observer may favor 
recording half tenths. Thus, if "... the maximum error of any function 
due to errors in the variables is the sum of the errors due to each 
separately" as stated by R. Gibbs (1929, p. 77); then in the course of 
50-100 years the sum of the observers' "personal errors" and of the 
"errors of the interval" could be significant. 

Errors in the determination of mean sea level which originate in 
the tabulation of mean values from ie continuous gage curve (excluding 
errors in calculation) may be classed under three headings: (a) samp- 
ling error, (b) errors of interval, and (c) errors of interpretation. 
These errors are probably of insignificant size for short term calcu- 
lations, but their effects taken over a period of fifty or more years 
may be appreciable—apparently no studies as to their significance have 
been published. 

The gage record of an automatic tide gage is a continuous sample of 
the fluctuating surface of the sea; if the effects of the disturbing fac- 
tors were removed from the trace of the gage, the remaining curve would 
be a representative sample of mean sea level. If the effects of the dis- 
turbing elements and instrument errors are not removed from the gage rec- 
ord, a primary sampling error occurs. Furthermore, monthly and yearly 
averages are not compiled from the irregular gage trace, but are compiled 
from hourly tabulations taken from the curve. The discrepancy between 


the hourly values of sea level and the values obtained from the integrated 


area under the curve may be called the secondary sampling error. 
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As monthly means are computed from the hourly tabulation, yearly 
means from the monthly, and 19-year means from the yearly; the samples 
become less and less representative of the actual elevations of the sea 
(this assumes that the effects of the disturbing elements have been re- 
moved, and that the discrepancies are due solely to using averages based 
on hourly tabulations and other averages rather than being based on the 
continuous curve). The error arising from these procedures is a sample 
averaging error. 

The errors of. interval in the tabulation of mean sea level occumia 
two places: (a) in determining the factor used to reduce the gage trace 
to a fixed datum, and (b) in making the hourly tabulation from the trace. 
The values for the height of the tide at the time of the staff gage read-= 
ings are measured (to tenths and half tenths of a foot) from the prelim- 
inary datum line to the curve with a reading scale graduated to the same 
scale as is used on the tide gage. This height (the relative height) is 
subtracted from the staff gage reading (the absolute height); the dif- 
ferences are summed and divided by the number of readings to obtain the 
mean difference. The mean difference, together with a econstanv [or tame 
datum, is added to the preliminary datum setting to obtain the correct 
height of the base line of the gage curve. The hourly readings are then 
measured to tenths and half tenths from the base line to the curve. 

The gage trace (if the stilling well inlet is not small enough) is 
not a smooth curve, but, instead, is an irregular curve with small sharp 


fluctuations ("saw teeth") due to the action of waves and swells and 
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larger irregularities due to seiches (periodic oscillations of the water 
meee ete irregularity of the trace makes tabulation difficult. There- 
fore, "For use in the determination of tidal datum planes it is prefer- 
able to consider a smooth curve through such irregularities and tabulate 
the hourly heights directly from this smooth curve" (Marmer, 1951, p. 41). 

The error of interpretation is the algebraic sum of the differences 
between the actual tidal curve and the smooth curve. Another source of 
interpretative error takes place when an interpolation is made for breaks 
in the tidal record. As in the previous case, the error is the algebraic 
sum of the differences between the interpolated curve and the actual sea 
level. 

The influence of instrument and operator error is relatively small 
in modern records; nevertheless it is still too large if accurate uplift 
rates are to be computed. Rates of postglacial crustal movement in 
Fennoscandia illustrates the need for very accurate measurements. ‘The 
rate of uplift 87 miles from the former ice center is 2.68 feet per cen- 
tury (Kad#ridinen, 1953, p. 59, Fig. 14) or 0.0% foot/year (0.00 foot/ 
month), and only 450 miles from the ice center uplift ceases; the average 
rate is 0.75 foot/100 miles/100 years. ‘The small magnitudes of these up- 
lift rates emphasizes the fact that the cumulative effects of the various 
factors which distort the true’ value of mean sea level must be considered 
and eliminated, or at least compensated for, if accurate rates of uplift 


are to be determined. 
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LAKE LEVEL GAGES 


Records of pairs of lake-level gages have been the basis for crustal 
movement studies in the lake plateau of Finland and the Great Lakes re- 
gion of North America. The Finnish lake study of A. Siren in 1951 Gta 
lized the principle that the means of lake-level elevations, corrected 
for meteorological effects and the effect of slope due to discharge, 
would provide a level surface for the comparison, over a number of years. 
of two gages which were remote from each other. When the average values 
of the gage differences are plotted against time and the points connected 
by a straight line, the slope of the line gives the average land upliig. 

The calculation of rates of uplift in the Great. Lakes area by various 
investigators have followed the method of G. K. Gilbert (1896-97)—using 
the procedures described on pages 7-8, which briefly, is that elevations 
of pairs of lake gage stations are measured from a level surface (the 
lake surface), and the differences between the elevations of the gages 
plotted for @ period of time. Jf 4 Hine of best-fit through ther pom 
has a slope, land uplift exists. When this method is used, no corrections 
are made for the effects of disturbing agents on the lake surface nor are 
corrections applied to the gage readings; it is assumed that the lake 
surface is level during the summer months (Comstock, 1876, D.5 $0 LOG2.aspe 
595) (see Appendix I). The fact that each rate calculated by this method 
depends upon the records of two gages brings out the importance of detect- 


ing and removing those factors which distort the recording of true lake 
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level at each gage station. 

Lake level gages and the procedures used to convert their recording 
curves to mean lake level elevations differ only in minor detail from 
tide gages and tide gaging procedure; therefore, the discussion of the 
factors which influence tide gage records also applies to lake-level 


records. | 


Tides 

Astronomical tides on those lakes large enough to exhibit detect- 
able tides are measured in inches, whereas ocean tides are measured in 
feet. In the Great Lakes, Defant (1958, Table 4) shows Lake Erie as 
having a spring tide of 8.0 cm (0.26 ft), Lake Michigan of 7.3 cm (0.24 
ft), and Lake Superior of 5.9 em (0.19 ft). 

The effects of tides of this magnitude should be very small, unless 
the lake-level gages are located in shallow water or in sites with con- 
verging shores; under these circumstances the tidal range would be in- 
creased. Rates of crustal movements in lake regions are computed from 
the differences in elevation of two gages; if dredging, construction of 
breakwaters, etc., occurred at one or both of the gage sites, it is pos- 
sible that the change in the range of tide would appear in the gage 


height differences to be incorporated in the rate of uplift. 


Meteorological Effects 
The influence of wind and barometric pressure differences was dis- 


cussed in many of the early papers dealing with fluctuations of the Great 
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Lakes (e.g., Dwight, 1622; p. 96; Clinton, 1@27, poe 292) 1295; Wheres 
1831, p. 214; Hall, 1843, p. 410; Mather, 1848, p. 1617; Whittlesey, 
1859; Lachlan, 1855, pp. 165, 168, 172). ‘These investigators of the late 
18th and early 19th century observed the increase in the depth of water 
which occurred on the eastern shores of the lakes with prevailing west- 
erly winds; were aware that barometric pressure differences piled up 
water in the area of lower pressure, and that this created seiches; ‘and 
described the effects of shallow water and converging shores on wind 
slopes and seiches. Later studies of wind slope and barometric effect 
were concentrated on Lake Erie which produces greater wind slopes than 
any of the other lakes due to its shallowness, its long narrow shape with 
pointed ends, and its alignment with the prevailing winds (Blunt, 1897; 
Henry, 1902; Hayford, 1922; Hayford, 1923; Hellstrom, 1941, pp. 115-128; 
Keulegan, 1954; Hunt, 1958; Gillies, 1959). 

Wind slope in the Great Lakes, as in the ocean, varies directly 
With the square of the wind velocity and inversely with the depth of 


water (see p. 51). In the formula (4) for setup 
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S = eH 


the undisturbed depth of water H is taken as the average depth of water 
along length L (e.g., in Lake Erie, H equals 58 feet); however, the es- 
tablishment of "thermoclines" in the Great Lakes during the summer months 


may require a modification of the usual definition of H. 
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During the summer months in the Great Lakes a thermal stratification 


of the water in which an upper layer of " ... more or less uniformly warn, 


circulating and fairly turbulent water termed epilimmion, "is divided 


Prom. i a deep, cold and relatively undisturbed region termed the 
hypolimmion...'" by the "thermocline ... defined as the plane of maximum 
rate of decrease in temperature ..." (Hutchinson, 1958, pp. 427-28). 


(Eckart, 1960, pp. 69-71). ‘The depth of the thermocline is about 50-100 
feet in Lake Huron and in the deeper portions of Lake Michigan and 20- 

50 feet in the shallow southern basin of Lake Michigan. (Ayers, et al., 
io, pp~s9, 41) 603 1950,pp. 5, 415 54; 80).. “he thermocline in Lake 
Erie is similar to the thermocline in the southern basin of Lake Michigan, 
usually being from 25-45 feet deep in the more shallow western two-thirds 
of the lake and 45-60 feet deep in the eastern third. At times during 
the summer the epilimmion extends to the bottom (up to 60 feet) in the 
central and western parts of Lake Erie, and maintains this depth across 
the deeper eastern third of the lake (Powers, et al., 1959, pp. 150-164). 
Presumably Lake Ontario has about the same depth to the thermocline 
(20-100 feet) as occurs in the deeper portions of Lake Michigan and Lake 
Huron (Hachey, 1952, pp. 326-328). 

The significance of the depth of the thermocline is that "... the 
boundary between the two layers acts as if it were a temporary bottom" 
(Hellstrom, 1941, p. 112). In other words, eine the summer when the 
lakes are stratified, the depth of water H in formula (3) could be con- 


sidered as the average depth to the thermocline rather than the average 
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depth to the actual bottom. This concept is reinforced by observacione 
made on the Scottish Lochs Ness and Gary (Murray, 1908, pp. 416-417; 
Wedderburn and Watson, 1908-09, pp. 629, 635; Wedderburn, 1909-10, pp. 
bia ot, O10) Lake Georce Ncw ork ((hengmudr, 1938, pp. 121-123) and 
the Sakrower and Walchensee, Germany (Hutchison, 1957, pp. 287-288). 

If the thermocline acts as a quasi-bottom, and since the depth ite 
the thermocline is approximately the same in all of the Great Lakes, the 
amount of deflection between wind direction and surface current for Lake 
Erie during the summer months should be representative of the amount of 
deflection for the other lakes. 

H. L. Langhaar (1951, pp. 279-280) has suggested that the total wind 
slope in lakes is made up of two component parts, a "statical tide" and 
a "dynamical tide." He declares that: 

The tide at the leeward end of the lake is the superpo- 

Sition of the> tide due. to the seiches ‘and. the statieal tide 

that. the wind would maintain if it persisted indefinitely. 

The tide due to the seiches will be called the "dynamical 

tide." The total height h of the tide at the leeward end of 

the lake is h =h, + hg, .-. the term hg {dynamical tide] 

varies periodically with time. 

The maximum value of the ratio hg/h, depends upon the 
planform of the lake, and upon the rate at which the wind 

develops. 

Differences in elevation (corrected for barometric pressure effect) 
which are recorded by two gages over a number of years. consist of a com- 
ponent resulting from the differences in the average height of water- 


level fluctuations produced by seiches at one gage being subtracted from 


the average height of water-level fluctuations produced by seiches at the 
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other gage (the "dynamical" portion), and a larger component made up of 
the net differences in water-level caused by the force of the wind (the 
"statical" portion). The magnitude of both components would be influ- 
enced by changes in the gage site environment occurring over a number of 
years (see pp. 58-59, 73-75). 

The barometric pressure effect on the Great Lakes is of less impor- 
tance than it is’on the oceans owing to the smaller pressure gradients 
which exist over the lakes. Over the summer month period the differences 
in barometric pressure, and thus the barometric effect, between two gage 
sites on the Great Lakes has an order of magnitude of 10°* or ror 
whereas the order of magnitude of the net wind slope during the summer 
months is 1071 (see Table 6, Appendix II, p. 179), e.g., if the effects 
of convergence and resonance are not taken into account, the barometric 
effect between Toledo and Buffalo on Lake Erie ranged between 0.002 foot 


and 0.017 foot for the summer months, 1950-59. 


Effect of Lake-Level Gage Locations 

Lake-level gages, like their counterparts in the oceans, are usually 
located near population centers. This means that almost all gages are in 
harbors or bays, within river mouths, or within a river mouth located in 
a harbor or bay; e.g., in Lake Ontario the gages at Kingston and Cape 
Vincent are situated on the north and south channels of the St. Lawrence 
River; in Lake Erie, Buffalo is at the very narrow eastern end of Lake 


Erie and at the entrance of the Niagara River and Toledo lies at the 
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apex of the converging shores of Maumee Bay and at the mouth of the 
Maumee River; in Lake Michigan-Huron, Collingwood, Ontario is in 
Nottawasaga Bay, Thessalon, Ontario, is at the western end of North 
Channel and Mackinaw City is on the Straits of Mackinac; in Lake Supe- 
rior the key gage at Point Iroquois is at the mouth of the St. Marys 
River at the southeast end of Whitefish Bay; Port Arthur is within 
Thunder Bay at the outlet of several rivers and Duluth is within Supe- 
rior, Bay at thevoutlet, of the Sv. Lovis River. 

The discussion on pages 55-59 of the effects of converging shores, 
shallow water, and resonance on the piling up of water by wind slope and 
barometric effect also applies to lakes. Manmade and natural modifica 
tion of the shores and underwater topography during the last 100 years 
have probably been'ias extensive in the Great Lakes regions as in ocean 
ports, which would result in a progressive change in lake level eleva- 
tions being incorporated in the gage records. 

Recent advances in telemetering equipment used in oceanography and 
Space satellites suggests that certain key water-level gages on each of 
the Great Lakes could be located well off-shore in deep water; from 
these sites wind velocities, barometric pressures, temperatures and 
water-surface elevations could be automatically determined, transmitted 
and recorded. The location of water-level gages away from shallow water, 
Days: Harbors; “Cl. woUldmadcereasemthneseclLects Om Soh certeotes and reso- 
nance; thus the errors caused by meteorological effects, tides, river dis- 


charges and man-made changes would be greatly reduced. ‘The recording of 
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temperature, wind and barometric pressure data at each of the key gage 
sites would provide the necessary information for the correction of the 


remaining meteorological effects. 


Instrument and Operator Errors 

Errors in the recording, correlating and reduction of lake-level 
gage records are of the same nature and order as those found in tide 
gage records; therefore the discussion of instrument and operator errors 
on pages 60-67 will also apply to lake-level gaging. Again it must be 
emphasized that rates of land uplift are computed from records which 
include those taken 45 or more years ago when neither the instruments 
nor the procedures were as refined as those of today, and errors which 
are very small in modern records are of much greater proportions in the 
older records. 

The “error of the interval" of lake-level records will be smaller 
than that found in tide gage records because both staff gage readings 
and tabulated values are given in hundredths of a foot rather than tenths 
or half-tenths of a foot. The errors of interpretation caused by the 
replacement of the actual recorded trace by a smooth curve is much re- 
duced or absent in lake-level records owing to a virtual lack of "saw- 
tooth" fluctuations in the trace, apparently the result of better damp- 
ing of short period waves. 

In an effort to obtain some concept of the magnitude of the secondary 


sampling error caused by obtaining a mean lake level from hourly tabu- 
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lations of the gage record rather than from the continuous curve itself, 
a comparison was made between a monthly mean lake level as published by 
the U. S. Lake Survey and a monthly mean level computed from the contin- 
uous curve of the same gage record. 

The lake-level gage record from Toledo, Ohio, for September, 1959, 
was used for the comparison. The area under the continuous curve, which 
was measured with a planimeter, was divided by the length of the base 
dine in order to obtain the mean altitude of the lake surface curveuae 
bove the bottom line of the gage record. The mean altitude of the curve 
was converted to feet (1.618 feet) and added to the mean elevation of 
the bottom line of the gage record of 570.512 feet (obtained from staff 
gage readings). This resulted in a mean elevation of 570.131 feet for 
the lake level—the mean elevation for September, 1959, published by the 
U. S. Lake Survey was 5/0.13 feet. 

Despite the fact that only one comparison of one gage was made for 
one month of one year, the close agreement of the two means suggests that 


the secondary sampling error for monthly means is small. 


MEASUREMENT AND ERRORS 


A review of the literature concerning the results of water-level 
gaging and the determination of rates of land uplift reveals that often 
the distinction between precise measurements and accurate measurements 
is not made, that statistical results are offered as proof of a state- 


ment without other supporting evidence, and that the Method of Least 
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Squares has been used without the prior removal of constant and system- 
atic errors. These points suggest the need of a brief discussion of 
measurement and the influence of errors. 

Measurement has been defined by N. R. Campbell (1957, pp. 267, 
524) as, "... ‘the process of assigning metblind to represent qualities," 
and "the primary object of measurement is to find a way of assigning to 
each of an ordered series of magnitudes a numeral, so that to each mag- 
nitude is assigned one and only one number and so that the order of the 
numerals is the order of the magnitudes...." 

Numbers are assigned "by comparing systems with the standard 
series." Each "system" has a hypothetical true magnitude which is re- 
lated to the true magnitudes of other "systems" by a numerical law 
(equation of condition) which is capable of being represented by an 
analytic function. 

If measurements are carried out with the most accurate instruments 


available, where ". 


- by the most accurate 'instrument' is meant, not 
merely a piece of apparatus which would ordinarily be called an instru- 
ment, but any arrangement whatever which permits measurement" (Campbell, 
eS i eyes Die Behe and with no discernable error of method, the observed 
magnitude may fail to approach the true magnitude by a maximum error E 
which is characteristic of the instrument of measurement and is dependent 
upon the smallest "step" of the instrument. 


When measurements are carried out with the above qualifications it 


will be found that a series of measurements will fail to agree within 
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fixed limits. The range of the inconsistency is determined by the error 
of consistency which N. R. Campbell (1957, p. 475) defines as: 


Nothing but errors of method magnified until they can: be 
directly detected by experiment. The magnification is ef- 
fected through the equation of condition, and results either 
from the addition of several partial errors of method, each 
of which could not be detected separately, or from the trans- 
ference of the error from a magnitude less accurately measur- 
able to one which is more accurately measurable. 


Campbell (1957, p. 509) states further that: 

The fundamental fact in the whole theory of errors of in- 
consistency {consistency] is that the true value of a complete 
collection of inconsistent measurements on a single magnitude 
is the arithmetic mean. If the collection (still consisting 
of measurements on a single magnitude) is incomplete, then 
we cannot determine the ‘true. value, that. is really all thar 
there is to be said about it. It may be convenient for some 
strictly limited purpose to express the results by a single 
numeral, and, if that is so, we shall probably select again 
the arithmetic mean as that numeral; but it cannot be too 
strongly insisted that the selection of that numeral does not 
imply a belief that it is the true value. 

Measurements are usually divided into two general classes; (a) di- 
rect measurements which are made by observing the comparison of a stand- 
ard with the system being measured, and (b) indirect measurements which 
are obtained by computation from direct measurements. It is apparent 
that the errors which affect indirect measurements are functions of the 
component direct measurements. 

Errors are usually classified as: (a) gross errors ("blunders," 
"mistakes"), (b) constant errors which have the same effect (maintain 


the same sign and magnitude) on all the observations of a particular 


series of observations, (c) systematic errors whose algebraic sign and 
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magnitude are determined by a fixed relation to some condition (being 
based on a law, they may be detected if the law and its coefficients are 
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known), and (d) accidental ("erratic," "residual," "experimental") errors 
which are the errors that remain after the gross, constant and systematic 
errors have been removed from the observation. Accidental errors are 
independent errors which are small in magnitude, are as likely to be 
positive as negative, and are much more likely to:'be small than large. 
Accidental errors are the only errors amenable to adjustment by the 
Method of Least Squares (Holman, 1904, pp. 4-10; Wright and Hayford, 
POO DOTI =O, 415,005, 25-2705 Goodwin, 19207 pp. (-123 Blunt,” 1941, 
Pewee yeeoeers, 955, pp. >-/; Anderson, 1955; p: (3; Rainstord, 1957, 

pp. 1-5). 

The customary classification of errors in the preceding paragraphs 
would be encompassed by "errors of method" [if errors of method are con- 
sidered to include the magnified errors of method (errors of consistency) 
in N. R. Campbell's (1957, pp. 267-294, 437-521) theory of measurement 
and error]. If errors of method exist in the observed magnitude the 
numerical law which relates that magnitude to other magnitudes cannot 
be determined; as long as these errors exist, the true magnitude, as 
represented by the arithmetic mean, cannot be found. When the errors 
of method are removed to the point that only errors of consistency re- 
main, the arithmetic mean may be used as the true value and possible laws 


may be formulated which relate the magnitude to other systems. 
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The concepts of accurate measurement and precise measurement are 
also clarified by Campbell's theory. An accurate measurement is one in 
which both the errors of method and errors of consistency are insignifi- 
cant$ thus the observed value will approach the true value and may be 
related by numerical law to other magnitudes. Precise measurement, on 
the other hand, is measurement in which the errors of consistency are 
low but errors of method remain. In this situation several different 
series of determinations would yield similar magnitudes, but the equa- 
tion of condition would be false and the magnitudes would not have the 
proper relation to other magnitudes. Because the goal of measurement 
is the determination of numerical laws, the taking of precise measure- 
ments (in that it reduces the error of consistency) is only a step in 
the determination of accurate values and is mot an endeanitsels. 

f 

The Method of Least Squares which is: 

A mathematical method for determining: (a) the most prob- 
able value of a single quantity from a number of measurements of 

that quantity; (b) the probable error of the mean value of a 

number of observations; (c) the best curve which may be drawn 

for a series of observed values of the ordinate over a range of 

values of the abcessa (AGI, 1957, p. 184). 

rests upon the mathematical demonstration that where each 
of a very large number of observations of any quantity is of 
the same quality as the others, the most probable value of 


the quantity is the one for which the sum of the squares of 
the residual errors (or corrections) is a minimum...(Mitchell, 


1948, p. 4b). 
The method has been a useful statistical tool for many years, and 
when used correctly provides the most probable values for observations. 


However, the fundamental requirement that first constant and systematic 
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errors be removed from the observations is, at times, not followed; thus 
erroneous conclusions are reached. 

The importance of this requirement is emphasized by the following 
quotations: 
T. W. Wright and J. F. Hayford (1906, p. 278) declared that: 


dhe devtectionsoalesystematicsor constant. errors necessarily 
involves least squares as a basis, but this must be supple- 
mented by something else, as the method of least squares 
deals with accidental errors only. 


D. Brunt (1941, p.v:) states that: 


It cannot be too strongly insisted upon that the methods of 
Least Squares cannot in any way improve upon the actual ob- 
servations. ‘The application of these methods to a large num- 
ber of carelessly conducted experiments cannot in general be 
expected to yield results as reliable as could be obtained 
from two or three carefully conducted experiments.... 


W. E. Deming (1943, p. 2) states: 


The principle of Least Squares provides a method for getting 
an adjusted value. It can be applied whether or not the data 
are worth adjusting, but the results are useful only when the 
data are good in the first place; no purely mathematical pro- 
cedure can make a good figure out of any number of bad ones. 
Data not in statistical control—i.e., not random, are not 
usefully adjusted. It is important to know when data are 
worth adjusting. 


and he continues (pp. 11-12) by declaring: 


The method of least squares can be applied to a single set of 
data, but no matter how carefully the least squares adjustment 
is carried out, the curve so fitted, or the observations so 
aijusped, JoO*snot have scientific validity unless there is 
other evidence at hand to show under what conditions the same 
or similar results will be obtained, and how these conditions 
are to be brought about and controlled. 
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The preceding discussion of errors and measurement has been made 
for three principal reasons: 

(a) To point out that measurement is fundamentally the determination 
of numerical laws which relate true magnitudes to each other, and that it 
is necessary to approach the true magnitude as closely as possible (make 
accurate measurements) if these laws are to be discovered. For example, 
when we determine rates of postglacial uplift by means of water-level 
gages, we are trying to find the numerical law which relates the gage 
readings or difference in gage readings to the magnitude of the uplift; 
if ithe true relationship between the two quantities) 1s to be: found. aa: 
is necessary to recognize and remove the errors which exist in the meas= 
urements. The mates of uplift are valid only to the degree) tiatmene 
errors are removed or compensated for. 

(b) To indicate the types of error which affect measurements and 
to bring out the fact that observations can be adjusted only after con- 
stant and systematic errors have been removed so that only accidental 
errors remain. 

(c) To emphasize the fact that the application of statistical 
tools, such as the method of least squares, cannot be used to produce 
valid results unless the assumptions upon which the determination is 


based are true. 


PRECISE LEVELING 


Precise leveling (also called "leveling of high precision," 
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"precision leveling" and "first-order leveling") is the determination of 
elevations of sequential points on the earth's surface (bench marks) with 
respect to each other and to a datum plane by means of a refined leveling 
instrument with a very sensitive spirit level to indicate the horizon. 
Careful methods of taking and processing the observed elevations, to- 
gether with a number of corrections (i.e., index, level, rod length, rod 
temperature, orthometric, curvature, and refraction corrections) reduce 
the magnitude of errors to certain prescribed limits and permit the 
establishment of accurate elevations. 

Three international meetings held in 1867, 1912 and 1936 set the 
limits of error which would be allowed for each class of leveling. ‘The 
1867 meeting defined precise leveling as leveling with an average prob- 
able error not in excess of 43 millimeters per kilometer and a maximum 
probable error of not more than 5 millimeters per kilometer. 

The meeting in 1912 prescribed the following classification (Rappleye, 
1948b, p. 150): 

Therefore the Seventh General Conference of the Inter- 

national Geodetic Association, still preserving unchanged the 

limits of error of 1867 for precise leveling, decides to place 

hereafter in a new class of leveling, to be termed "leveling 

of high precision," every line, set of lines, or net which is 

run twice in opposite directions on different dates as far as 

possible, and whose errors, accidental and systematic, com- 


puted by the formulas hereinafter given, do not exceed— 


= lmmperkmior the probable accidental error, 
or 


I+ 


1.5 mm per km for the mean accidental error; 
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O.2 mm per km for the probable systematic error, 
Or 
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O.3 mm per km for the mean systematic error. 

The Sixth General Assembly of the International Association of 
Geodesy issued new information as to the design of instruments and rods, 
methods of operation, computation and adjustments which were to be used 
for leveling of high precision. The Association also redefined levetimg 
of high precision as a method of leveling with a total probable error 
not exceeding 2 millimeters per kilometer—leveling with a total prob- 
able error exceeding 2 millimeters but not exceeding 6 millimeters per 
kilometer was classified as "precise leveling" (Rappleye, 1948b, p. 154). 

First-order leveling of the U. S. Coast and Geodetic Survey and the 
U. S. Lake Survey includes leveling in which the level lines are divided 
into 1 to 2 kilometer sections and the results of a forward and backwara 
leveling over a section does not differ by more than 4.0 millimeters 
times the square root of the length of the section in kilometers (4.0 mm 
Vk), or its equivalent: 0.017 times the square root of the length of the 
section in miles (Mitchell, 1948, pp. 45-6). 

H. S. Rappleye (1948a, pp. 1-2) states: 

First-order leveling by the United States Coast and 

Geodetic, Survey began with the transcontinental line of levels 

in 1878. Previous to that time the Bureau had done some lev- 

eling, but it was used mostly to control trigonometric leveling 

and, while it served. its purpose, it was not of ayhigh grade or 


accuracy compared with the standards. of today.... 


and 


The first-order leveling done by the United States Coast 
and Geodetic Survey since 1899 falls within the limits prescribed 
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for "leveling of high precision" at the Hamburg meeting [1912]. 

In the 1912 adjustment of the first-order level net the aver- 

age probable accidental error per kilometer was plus or minus 

O./71 millimeter and the average probable systematic error was 

plus or minus 0.08 millimeter. 

Finnish Uplift Rates By Precise Leveling 

Finland presents an almost ideal situation for the determination of 
extensive non-volcanic land uplift by precise leveling. Finland's bor- 
ders are from c 67 miles (108 km) to c 455 miles (574 km) to the south- 
Saou. swe sl Ormer Fennoscandian ice center, thus Finland. is very close 
to the region of maximum uplift. 

The underlying igneous and metamorphic bedrock crops out in many 
areas which allows the placement of numerous bench marks in the stable 
bedrock; thereby reducing the errors due to shifting bench marks. Prob- 
ably the most important factor, other than proper instruments and tech- 
niques, in the determination of accurate absolute rates of uplift is 
the fact that the level net is tied into either 14 primary tide gages 
(first precise leveling) or 12 tide gages (second precise leveling). 

The average distance between the 14 tide gages, measured along the 
leveling lines, was c 78 miles (126 km); the shortest distance was c 36 
miles (58 km); and the longest distance was c 138 miles (222 km). The 
distances between the 12 tide gages of the second precise leveling were: 
(a) average distance c 81 miles (130 km), (b) shortest distance ec 55 
miles (88 km), and (c) the longest distance c 138 miles (245 km). ‘The 


point on the level line polygon which is farthest from a tide gage is c 
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197 miles (317 km) from the gage, and the point on a level line common to 
both precise levelings which is farthest from a tide gage is c 182 miles 
(293 km) from that gage. 

The information in the following paragraphs regarding the Finnish 
precise levelings of 1892-1910 and 1935-1955 is from E. Ké&ridinen's 
"On the Recent Uplift of the Earth's Crust in Finland" (1953, pp. 31-64). 

The leveling net of the First Precise Leveling consisted of 11 closed 
polygons with a ‘total length of Leveling lines in the principal nero. 
2464 miles (3967 km). ‘The smallest polygon was 71.4 miles (115 km) in 
length and the others ranged from 142.2 miles (229 km) to 502 miles 
(808 km) in length. Approximately 3,000 bench marks were placed in bed- 
rock, "immovable" boulders and stone foundations at spacings of 0.94- 
W2y mides (.5=2 km). 

"From the closing errors of the polygons, the greatest of which was 
66.30 mm it was computed that the mean error of the levelings was + 1.24 
mm/km" (1953, p. 31). 

The second Precise Leveling consists of sasmeuwork Of J5,closed 
polygons with circumferences ranging from 146 miles (235 km) to 543 
miles (874 km). After World War II, 2843 miles (4577 km) of leveling 
lines.,and 12 tide gage stations remained within Winnish territory. 

ae The greatest closing error, without taking into account the 
refraction and the land uplift correction is 22.10 mm. ‘The mean error 


of the Second Leveling, as computed from closing errors, is + 0.45 mm/km" 


(1955 DO) 
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The leveling lines common to both the First and Second Leveling 
totaled 2309 miles (3717 km). Although about 70 per cent of the old 
bench marks were found (c 46 per cent of which were on rock), the land 
uplift calculations were based almost entirely upon 900 highly dependable 
rock bench marks; in a few cases bench marks on boulders were used. 

After the necessary computation, corrections and adjustments were 
made to the leveling lines and nets, the mean errors of the yearly land 
uplift values were calculated for 28 tie-points, and a value of "... + 
O.4 mm as the average mean error of the yearly land uplift obtained by 
precise leveling" (1953, p- 59) was computed. The rates of uplift of 


the 12 tide gage stations are listed on page 93 of this paper. 


| Great Lakes Region—Precise Leveling 

Elevations on the Great Lakes have been determined since the first 
complete leveling in 1875 by a combination of spirit-level measurements 
and waver-level measurements. U. S. Lake Survey leveling lines have 
originated at a bench mark in Rensselaer, New York, called Greenbush 
Gristmill (see Appendix I). The elevation of Greenbush Gristmill above 
sea level was determined by the U. S. Coast and Geodetic Survey in 1856- 
57, 1877, 1894, 1902, 1934, and 1955. Instrumental level lines are run 
between Greenbush and Oswego, New York (the first point of elevation on 
the Great Lakes proper, and between each of the other Great Lakes. Water- 
level determinations are used to carry the elevations from the eastern to 


the western ends of the lakes (Comstock, 1882, pp. 595-609). 
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The geographic factors inherent in the location of the’ leveling 
lines used to determine rates of land uplift in the Great Lakes Region 
are in distinct contrast to the favorable,situation which exists in 
Finland. 

Elevations in the Great Lakes Region are based on leveling lines 
run forward and backward from the tide gage at New York City to the Green- 
bush Gristmi I ‘benchmark, ~acdi stance lot ¢ 147 miles (23% km) ; the level 
lines are then, continued by another organization to Oswego, New yorkjgcm 
Lake Ontario, a distance of 173-195 miles (278-314 km) depending upon 
the yroute veaken- 

From Oswego, the elevations are carried by alternating water-leveling 
and spirit leveling a distance of ¢ 1212 miles (1950 km) to Duluth, 
Minnesota. If the rates of uplift are to be calculated for the Lake 
puperior area, elevations for Ports Arthur, Ontario, vareruscd. serone 
Arthur (which is ec 185 miles [298 km] from Duluth) elevations were deter- 
mined by yet another organization, the Canadian Hydrographic Service, 
from a comparison with lake-level gage records at Marquette, Michigan, 
from 1907-1914. 

The total length of the leveling lines, both water-level and spirit 
level, from the tide gage at New York City to the farthest point in the 
level net (Port Arthur) is approximately 1720 miles (2768 km). ‘This 
situation is to be compared with the Finnish leveling lines suitable 
for determining rates of uplift which total 24309 miles (4716 km) and are 


tied into le tide gages. The point which is farthest from a tide gage in 
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the Finnish level net common to the First and Second Levelings is ec 182 
miles (292 km) from the tide gage. 

The tide gage records used in the Finnish level nets are corrected 
for the eustatic rise of sea level, meteorological effects, and instru- 
ment and operator error, as well as for astronomical tides. The tide 
gage record at New York City on the other hand, is of a sufficiently 
long period to eliminate the majority of the astronomical tides and 
meteorological effects, but the other influences remain uncorrected. 
Elevations from lake-level gages records are uncorrected observed ele- 
vations of the lake surface. 

A final comparison between the two areas may be made on the basis 
of the distances of the areas of uplift from their respective Pleistocene 
ice centers. The Finnish tide gage at LeppXluoto, which is c 8% miles 
(134 km) from the ice center, undergoes an uplift of 2.69 feet/100 years 
(8.80 + 0.33 mm/year), and the tide gage at Hamina, 346 miles (557 km) 
from the ice center, is uplifted at a rate of 0.82 foot/100 years (2.70 
Su Ove hy. mm/year ). The distances of these gages from the ice center are 
to be compared with the distances from the Laurentian ice center to points 
in the Great Lakes region. For example, the distance from the Canadian 
ice center to Oswego, New York is c 800 miles (1288 km); from the ice 
center to Duluth, Minnesota, is c 870 miles (1400 km); and from the ice 
center to Milwaukee, Wisconsin, on Lake Michigan the distance is c 1110 
miles (1790 km). 


A concise summary of the caution needed in using precise leveling to 
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7 
detect rates of crustal movement is stated by F. Nemeth (1960, p. 53), 
who: declared: 
Geodesists have long drawn attention to the fact that 
leveling data should be used as evidence of crustal movement 

only with caution. The differences in altitude of leveling 

sections are mostly veliable only to the order Of saqmi la 

meter because of measurement error; however, movements of 

ground and structure can multiply the measurement error. A 

difference in compared altitude values can also stem from 

adjustment of different leveling networks. For investigation 

of changes -of, levelseonily onesand the same) network, spec tfaiky 

stabilized for this purpose and always measured by the same 

methods is suitable. 
FENNOSCANDIAN—GREAT LAKES ANALOGY 

Numerous comparisons have been made of the postglacial land upliga 
of Fennoscandia and the Great Lakes Region (e.g., De Geer, 1892; Guten- 
berg, 1933, 1941, 1954; Flint, 1957) which stressed the over-all eimi- 
larity of uplift in the two areas but did not compare in detail the re- 
lationships between maximum depressions, distances from their respective 
Tce Ceneers,, ana the raves On modern Uplat tc. 

If the two regions are analogous, it would seem that an examination 
of the rate and extent of modern land uplift in Fennoscandia would aid 
in locating, at least approximately, the zero isobase of modern land 
uplift in the region which had been covered by the Laurentide ice sheet. 
Although the Fennoscandian ice sheet was smaller than the Laurentide 
sheet, both areas had approximately the same thickness of ice, 10, 000+ 


feet (2,500+ meters), and consequently about the same depth of depression 


of the crust. Therefore the distance from the ice center to the zero 
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isobase in the Great Lakes-Hudson Bay area should be roughly equal to 
the distance from the Fennoscandian ice center to the zero isobase at 
Leningrad, U.S.S.R. 

The Fennoscandian uplift values in Table 2-A are from Kaaridinen 
(1953, p- 59, Fig. 4); distances were measured from the location of the 
axis of maximum ice thickness in eastern Sweden (Flint, 1957, p. 368, 
Plate 5) to each of the twelve tide gages. The distances from the 
Laurentide ice divide to points in the Great Lakes-Hudson Bay region 
(Table 2-B) were determined from the Glacial Map of Canada (Geological 
Association of Canada, 1958). 

A plot of the uplift of the Finnish tide gage stations, determined 
by precise leveling (in feet per 100 years), against their distances 
(in miles) from the Fennoscandian axis of maximum ice thickness reveals 
that the distribution of points lies almost in a straight line—the slope 
of the line fitted by the method of least squares gives a modern rate of 
Oda te Ot 350 feet /450 miles/100 Veer One Onn) foot/100 miles/100 years. 
Ud Seer BUC, OL 1 Oats foot/100 miles/100 years might then be assumed to be 
of the proper order of magnitude for the modern Laurentide land uplift. 

If the analogy with Fennoscandia holds, the zero isobase of the land 
uplift would follow along a line extending down from the western fifth of 
Hudson Bay to the southwestern shore of James Bay (just southwest of 
Moosoonee, Ontario) then to the vicinity of Kempt Lake, Quebec, crossing 
the St. Lawrence River about 20 miles northeast of Quebec, Quebec. The 


relationship of Churchill, Manitoba, to the zero isobase, i.e., Churchill 
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TABLE 2-A 


FENNOSCANDIAN UPLIFT RATES AND DISTANCES 


Uplift in 


Distance from Axis of | ‘Yearly Uplift in 


age | Maximum Ice Thickness mn and Mean Error Feet/100 

| ___| Kilometers | Miles | By Precise Leveling Years _ 
1. Leppaluoto 134 83 S700> +* 059 2.69 
2. Vaskiluoto 141 87 So (4) 0.00 2.08 
32) Horankel ito 146 91 9.02 + 0.43 eile: 
4. Kaskinen 184 114 71,06 + 0.32 2.34 
Dae Nem 190 LASS 8.524 0.55 260 
6., Toppila 196 ee 8.25 + 0.43 Pee 
7. Mantyluoto 258 160 6.96 + 0.25 2gi2 
8, Rauma 296 184 6.82 + 0.28 2.08 
9. Ruissalo 383 238 52 50 028 1362 
10. Hanko . 4.60 286 5.50 + 0.55 102 
1, -wedeinics 503 pile 2.90 ESOCes 0.88 
12. Hamina 5D i 346 2.70 + 0.24 0.62 
Leningrad, .USSR Hane LLG O:.00 0.00 


TABLE 2=B 
DISTANCES FROM THE LAURENTIDE ICE DIVIDE 
Approximate Distance 


from Tce:Dividem 


| Kilometers | 


1. Father Point, Quebec 645 4.00 
2. James Bay (Moosoonee [co1lis ]), Ont. 10 h9Q 
4. Kempt Lake Quebec 805 500 
4. Ottawa, Ontario 1040 645 
5. Churchill. Manitoba 1060 660 
6. Lake Ontario (Kingston, Ontario) 1210 750 
7. Lake Huron (Manitoulin Island, Ontario) 1255 780 
8. Lake Michigan-Huron (Mackinaw City, Mich.) 1380 G55 
9. Lake Superior (Duluth, Minnesota) 1400 870 
10. Lake Erie (Buffalo, New York) ? 1430 885 
11.’ Lake Michigan (Milwaukee, Wisconsin) 1790 Cae) 
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is c 100 miles farther away from the ice center than the presumed zero 
isobase, may support those investigators who deny the existence of land 
uplift at Churchill (Tyrrell, 1896; Johnson, 1939; Cooke, 1942; Williams, 
1949) as opposed to those who claim that up to c 2 meters (6.6 ft) of 
uplift per century is taking place at Churchill (Bell, 1897; Gutenberg, 
1941, 1942, 1954). 

It is possible that the ice which covered Churchill came from the 
Tcesdiviae, tosthe west of Hudson Bay; if this were the case, Churchill, 
being about 285 miles from the western ice divide, would have a modern 
land uplift of about one foot per century (again applying the modern 
Fennoscandian rate of uplift). 

The tide gage at Father Point, Quebec (being about 400 miles from 
the Laurentide ice divide), is located within the proposed area of up- 
lift and could be wndergoing an.uplift of c 0.5/ roay AN OTe: years (based on 
the modern Fennoscandia rate). The importance of uplift at Father Point, 
Quebec, lies in the fact that it is one of the tidal bench marks upon 
which the Canadian precise level net is based; furthermore it is the 
starting point for the determination of the new International Datum for 
the Great Lakes. 

The advisability of using the tide gage at Father Point may be 
Questioned for three reasons: 

(a) The gage is located within the area of postglacial uplift. 

If uplift is now occurring, and if corrections are not made for the up- 


lift, precise levelings taken a number of decades apart would show an 


Ql. 


apparent subsidence cr the Wand; 

(ob) the gage is situated near the small end of a "funnel" formed 
by the Gulf of St. Lawrence, the estuary of the St. Lawrence River and 
the St. Lawrence River. The convergence of the shorelines greatly in- 
ereases the influence of meteorological effects (see pp. 58-59, fierlisy 
which in turm makes difficult the determination Of an accurace meanmece 
level; and 

(c) the Father Point tide gage is on the eastern side of an ancient 
fault zone (Logan's Line), whereas the greater part of the Canadian level 
net and the whole of the Great Lakes region lies to the west of the fault 
zone. Although the fault zone is an ancient one, it would be prudenu mex 
to have the key tide gage of a leveling net and the greater part of the 
néetron Opposite sides ofia tauli zone: 

The modern zero isobase of postglacial land uplift cannot be in 
the known area of horizontality (i.e., south west of the Nipissing zero 
isobase). Moreover if the modern zero isobase were in the VICI Cy. com 
the southwestern shore of James Bay; then, as’ is indicated in Table 2-=B 
and Plate I, the postglacial uplift would take place northeast of a line 
which passes near the southwestern shore of James Bay and which curves 
eastward SO aS to Cross the St. Lawrence River just horth of Quebec, 
Quebec. Under these circumstances, ene maximum possible southern extent 
of postglacial crustal oon would include almost all of Lake Supe- 


rior,, less than the northern one-fourth of Lake Michigan, less than the 


oe, 


northern two-thirds of Lake Huron and less than the northern three-fourths 
of Lake Ontario. The minimum extent of the southern boundary of uplift 
on the other hand would be somewhere near James Bay—which would mean that 
modern land uplift in the Great Lakes region would be nonexistent. The 
actual southern limit of uplift should lie somewhere between these two 


\ 
extremes—probably closer to the Great Lakes than to James Bay. 


7 
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\ IV. PREVIOUS DETERMINATIONS OF RATES OF UPLIFT 


Rates of uplift in the Great Lakes region have been computed by var- 
ious investigators who used Gilbert's method of calculation based on the 
principle of "water leveling" (see pp. 7-11, 40, 68-69). Examples of the 
maghnuvude of vuplagrc in the Great Lakes are given in Table 4, which also 
indicates the direction of uplift between the gages of each pair of gagere 
The-Land is uplifted in a direction extending from the geographic toc_- 
tion of the first gage of the pair toward the second gage. 

The directions of uplift in Table 4 were found by projecting the 
isobases of former glacial lake shoreline features through the locations 
of the gages; the gage of a particular pair which was on the higher iso- 
base was considered to be the gage which was "upslope." Thus the direc- 
tions of land uplift in Table 4 are based on uplift which warped former 
shorelines, whereas the directions in Table 4 represent modern uplift 
calculated from lake-level gage records. 

If the directions of uplift between the pairs of gages in Table 4 
are compared with the directions of uplift between the same pairs of 
gages in Table 4, it will be found that the direction of uplift will be 
reversed for certain pairs of gages, e.g., Marquette-Duluth, Conneaut- 
Cleveland, Oswego-Charlotte, etc. In cases of discrepancy the directions, 


and wherefore: the raves, in lable, > must be anecorresge. 
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TABLE 3 


EXAMPLES OF PREVIOUS DETERMINATIONS OF RATES 
OF UPLIFT IN FEET/100 MILES/100 YEARS 


Lake Lake 
Superior Mich.-Huron 


| Max | Min | Max | Min | Max 


0.43% 0.399 |o.46e = 


Lake 
Ontario 


After 


Gilbert (1896- 
wT p. 636) 


Moore (1922, 
Dowel LOL) 


Gutenberg | 
(1933, p. 457) | 0.8" 0.12 


Gutenberg | 
ohio ieya) 061+ 90.040), 0,485 0.09" 10.55" 6.0 }1.34P 0.06” 


Moore (1948, | | | 
p. 700-701) oped O50 7 On 55 OOF 10 26 OnOSUa On ay> O1dios ON 


Land is uplifted in the direction of the second gage site. 


Harbor Beach-Escanaba 


a. Milwaukee-Port Austin Of 
b. Milwaukee-Escanaba p. Kingston-Cape Vincent 

e. Cleveland-Port Colborne q. Port Dalhousie-Kingston 

d. Charlotte-Sacketts Harbor r. Harbor Beach-Collingwood 
e. Marquette-Sault Ste. Marie s. Calumet Harbor-Harbor Beach 
f. Marquette-Duluth tee ror ovanley<rort Colborne 
g. Milwaukee-Harbor Beach u. Cleveland-Port Stanley 

h. Cleveland-Amherstberg v. Kingston-Oswego 

i. Cleveland-Buffalo w. Marquette-Port Arthur 

Jj. Oswego-Charlotte x. Marquette-Houghton 

k. Oswego-Toronto y. Conneaut-Cleveland 

lL. Duluth=Port Arthur z. Port Dalhousie-Oswego 

m. Michipicoten-Port Arthur aa, Oswego-Cape Vincent 

n. Calumet Harbor-Milwaukee 
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Karlier Discussions Of Errors 


As has been stated previously, the fundamental assumption which 
underlies water-leveling and the determination of Great Lakes rates of 
uplift is that the summer season mean lake surface is level and may be 
used as a plane of reference. The importance of this concept must be 
stressed, for obviously if the mean lake surfaces are not level during 
the summer months, comparisons cannot be made between the gage differences 
of two gages over a period of time due to the lack of a standard reference 
point. If gages are compared by means of a lake surface which is not 
bevel, tien the rates of uplitt which result from wthe*comparison’ are not 


Wich trite 


GILBERT'S INITIAL STUDY (1896-97) 

The problem of determining a level lake surface, as well as the ef- 
fects of the various other factors which influence the taking of accurate 
water-level measurements, was completely understood by G. K. Gilbert who 
made the first computation of modern rates of crustal movement in the 
Great Lakes. His understanding of the types of errors involved, their 
importance, and the means of avoiding or eliminating them was not approached 
by subsequent investigators. The following quotation is taken from Gil- 
bert's (1896-97, pp. 641-645) chapter "Plans For Precise Measurement"— 
plans which, if they had been followed and developed, would have largely 
eliminated the various types of error which prevent the accurate measure- 


ment of water-level elevations. 
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Gilbert (pp. 643-644) explained that: 


In bays and estuaries there are local temporary variations 
occasioned by the floods of tributary streams. 

There are solar and lunar tides, small as compared to 
those of the ocean, but not so small that they may be neglected. 
The wind pushes the lake water before it, piling 16 up 
on lee shores and lowering the level on weather shores. During 
great storms these changes have a magnitude of several feet, 
and the effect of light wind is distinctly appreciable. Even 
the land and sea breezes, set up near the shore by contrasts 
of surface temperature, have been found to produce measurable 

erLhecus: Om whe water Level. 

There is also an influence from atmospheric pressure. 

When the air is in equilibrium, sar sine ever occure, wie 
pressure’ is the same on all’ parts*or the. lake surface, and 
the equilibrium of the lake is not disturbed; but when the 
air pressure varies from point to point this variation of 
pressure is a factor in the equilibrium of the water Surtace, 
the surface being comparatively depressed where the air pres- 
sure iS greater and elevated) where 1G is less. 

When a storm wind ceases, the water not merely flows back 
to its normal position but is carried by momentum beyond, and an 
oscillation is thus set up which continues for an indefinite 
period. A similar oscillation is started whenever the equilib- 
rium is disturbed by differences of atmospheric pressure; and 
these’ swaying motions, called seiches, ..., persist for long 
periods. In fdct, they bridge over the intervals from impulse 
to impulse, so that the water, of the Great Lakes never comes 
to rest. 

These various influences work independently but simulta- 
neously; yandetheir: effects are blendéd “in thevactual osei dias 
tions of the water surface at any point. In using the water 
surface: for the purpose of precise leveling, it iswnecessary 
to take account of all such factors and make provision for the 
avoidance .or correction of the errors they tend to preduce. 

The gage employed for the determination of water height 
should be of some automatic type, giving a continuous record. 
This-is necessary in order that the study of the record may 
furnish data for the complete elimingmmhen of errors from Gide, 
seiches, and\landvandwseasbreezes. <a ier shouldsber installed 
as to be secure from settling... The height of its zero should 
be! readilyevyeritrabiles 

Near each station there should be at least three benches, 
constructed with special reference to permanence and stability. 
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They should be independent of one another and independent of 
other structures. 

Pressure of the air should be continuously recorded by 
a baragraph, carefully standardized. A wind vane and anemometer 
should give automatic records. 


Although Gilbert realized the importance of the factors which affect 
the determination of rates of uplift, he was limited, as in any pioneer 
study, by the data available. Gilbert attempted to avoid the effect of 
wind slope by using only those records which were taken on days of very 
light wind; he was no doubt unaware of the fact that wind-drift currents 
(hence the piling up of water) are the result of the effects of wind 
stresses which were applied up to 10-12 days before the day of measure- 
ment (Millar, 1952, pp. 536-337; Ayers, et al., 1958, pp. 112-115; Ayers, 
1959, pp. 4-5). 

Gilbert (pp. 637-638) discussed the other sources of errors as 
follows: 


The probable errors of the individual measurements are 
rather high, ranging from 14 to 50 per cent, and this suggests 
the possibility that the closeness of their correspondence may 
be accidental. It should be remembered also that at two or 
three stations there was reason to believe that the gage zeros 
were settling during the period in which the observations were 
made, and the results involve the doubtful assumption that the 
rate of settling was uniform. There is room for doubt as to 
the precision of the instrumental leveling; in only a few in- 
stances is the fact of duplicate measurements recorded, and 
single measurements are notoriously insecure. Error was 
doubtless admitted by ignoring the effects of barometric 
gradient. River floods may have introduced errors. ... There 
may also be personal equations of observers, especially as 
the gages at pairs of stations were not in every case of the 
same type. For all these reasons I am disposed to ascribe 
only a low order of precision to the deduced rate of change, 
and regard'it as indicating the order of magnitude rather 
than the actual magnitude of the differential movement. 
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"LEVEL" LAKE SURFACES 

The validity of the assumption that the Great Lakes water surfaces 
are level during the summer was questioned in 1897 by William T. Blunt, 
engineer of the U. S. Deep Waterways Commission, in "Effects of Gales on 
Lake Erie" (Blunt, 1897, pp. 155-168). 
Blunt (p. 156) declared: 


In the survey of the Northern and Northwestern Lakes the 
assumption was made that the mean surface of each lake was 
level within the limit of possible instrumental errors in trav- 
ersing its length, and all heights and gauges west of Oswego 
have been based on this assumption. The actual period of ob- 
servations: used: to transfer the level by lake surface was from 
May 11 to August 31, 1875; and while it is certain that the 
eminent officer then in charge of the ,survey took every pre- 
cautions GO obtain accurate results, it will sti liebe: ameicer 
of scientific if not of practical interest to have a verifi- 
cation of these elevations. When we consider the marked ef- 
fects of even light winds on the surface of Lake Erie, the 
very decided effects of strong continuous winds, and the extra- 
ordinary effects of gales, in connection with the fact that the 
great preponderance of winds is from the westward, the propo- 
sition that even the mean surface is level appears somewhat 
clouded; it at least requires verification. 


John F. Hayford, geodesist and engineer, completed an investigation 
in 1922 which also refutes the assumption that summer mean lake surfaces 
are level; the results of the investigation were published by the Car- 
negie Institute of Washington in Publication 41/, "Effects of Wind and 
of Barometric Pressure on the Great Lakes." 

The principal problem of determining accurate lake-level elevations 
was stated by Hayford (p. 2), who declared: 

As the investigation progressed, it gradually became more 


clearly evident that the largest and most serious errors 
encountered were those which arise from the fact that the 
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surface of any one of the Great Lakes at any given instant 
is not level except by accident. The surface has a slope 
av every poinu due to the influence of winds and barometric 
pressures. 


Hayford, for the most part, dealt with the short-term departures 
from level; however his discussion of tide gage records indicates that 
he also realized the importance of the seasonal effect of wind and baro- 
metric pressure (the conclusions reached apply equally well to lake-level 
gages). Hayford (p. 132) pointed out that: 


It, should, mot be overlooked in this connection that the pre- 
vailing winds and the prevailing barometric gradients tend to 
Deseeaconal . tOnbe repeated cach year, and thatwtherefore, the 
taking of a mean for several years is of only moderate effec- 
tiveness in reducing the error in the mean. The monthly values 
of mean sea-level at various tide gages support the statement 
by showing a seasonal variation, as a rule, and thereby inci- 
dentally indicating that the wind effects and barometric ef- 
fects are certainly decidedly appreciable in the monthly means. 


The final chapter of Hayford's book dealt with the "Application to 


Determination of Tilting of the Great Lakes Region." There Hayford 


(p. 134) said: 


The rate of tilting as derived was .OO42 foot per mile per 
century—an exceedingly small rate of change. The conclusion 
was derived from apparent changes of relative elevation of the 
water surface as measured at different gages on Lakes Michigan- 
Huron, Erie, and Ontario in different years. The amounts of 
change involved are of the order of 0.1 to 0.2 foot in a per- 
iod of 20 to 40 years. Evidently, when such small changes are 
in Question there is more chance of securing the necessary 
accuracy if corrections as large as those shown in Tables Nos. 
19 to 23, pages 80-96 of this publication, for barometric ef- 
fect and wind effects, are taken into account. 

The deductions of Gilbert are probably correct in the 
main. But a new investigation based on observed elevations 
of water surface corrected for wind effect and barometric 
effects would have greater accuracy and is desirable. 
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Sherman Moore's first determination of rates of crustal movement 
in the Great Lakes region (examples of which are given in Table 3) 
was. also published in 1922. In his first«paper, ‘Wilt of the-Hartheim 
Great Lakes Region," Moore devoted several paragraphs to the errors which 
affect lake-level gage readings; this discussion of errors is in sharp 
contrast to his second paper published in 1948 which contained only two 
sentences referring to errors (one remark is on bench mark stability and 
the other suggests that scatter of points may be due to varying wind and 
barometric conditions). 

One statement is particularly pertinent to the discussion of "level" 
lake surfaces. Moore (1922, p. 154-154) declared that: 

The lakes themselves are not always level. Prevailing winds 

and continued differences in barometric pressure cause tilt- 

ing of the surface. which may last for considerable periods of 

time. These inequalities’ will usually Dealance ane Tong per= 

Lod ol. Gime, but evens yearly meanris neveinece from them, 

effects. 

Apparently Moore did not realize that the second sentence of the quota- 
tion nullified the assumption—that the summer mean lake surface is 
level—upon which his method of determining uplift rates is based. 

Bene Gutenberg's initial paper on "Tilting Due to Glacial Melting" 
(1933) contains only one reference (on bench mark stability) to errors 
in the treatment on Great Lakes mates of uplitihe "Hisnsecond= paper, pub= 
lished in 1941, discussed the influence of the following factors upon 


the recording of tide gage elevations: meteorological effects, eustatic 


changes, water and shore conditions, local movements, effects of the 
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method of observation, effects of errors, and land movement. His obser- 
vations: relating to the measurement of lake-level elevations on the 
other hand were scanty. Gutenberg (p. 740) said: 
Absolute values (corresponding to the column 6 in Table 4) were 
MOpecalculaveoyeaceuley Brerini lucncedibyethe. generainchange in 
tiewake eve SduecstOepréecipi tation, flow sin rivers }ivertical 
movements of the region of the outlet, wind and other causes. 
imesddi tions whesGp.—(45Jeseids 
Apparently, there are meteorological conditions which affect 
the mean lake level at the various stations on Lake Ontario 
morescnean Thosesinaevhe rotherslakes and therefore aivnis to-be 
expected that the results for Lake Ontario are less accurate. 
It may be seen from the foregoing paragraphs that the existence 
of factors which prevent a lake surface from being a level surface was 
known to those investigators who used the lake surface as a datum plane 
in the comparison of pairs of gages in. their attempts to measure rates of 
crustal movement. However these investigators seemingly regarded the in- 
fluence of meteorological effects, ete., as constituting only a temporary 


disturbance of level and therefore continued to use uncorrected lake sur- 


face elevations in their calculations. 


OTHER ERRORS 

The existence of other types of error which affect the taking of 
lake-level elevations was known to those investigators who measured 
rates of crustal movement. This fact is demonstrated by the papers of 
Sherman Moore, engineer of the U. S. Lake Survey, who discussed the var- 


idus types of errors in his first paper on»rates of tilting (1922) and 
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in his gage histories. (1939-44). The following quotations from Moore's 

1922 paper point out the poor quality of the early measurements —measure- 

ments incorporated in the calculation of rates of uplift which, according 

to Moore (1922, p. 182), averaged 0.43 + 0.07 foot/100 miles/100 years. 
Moore °(1922 {pp 153-155 Gl) “stated “that; 


Practically, the determination is. complicated by errors 
arising from several causes. Gage readings are frequently 
incorrect due Gomistakes in reading “he eage, put moresire= 
quently) to-failure to check -the .zéro of the gage back to 
stable benchmarks sufficiently often. Benchmarks are not 
always stable, but are subject to local settlement. Bench- 
marks are frequently destroyed, and unless there.is more 
than one benchmark at the point, later records are not 
comparable. 

.. From 1860 the records are complete at some points, 
but until about 1872 the gages, in the majority of cases, 
were; poorly cered Tor, and. were not property rel orenceda 
fixed benchmarks. In the seventies more care was given to 
the gages, ’and levels connecting their zeros with faxed 
benchmarks were run at frequent intervals. Unfortunately 
the records of the actual gage readings are not available. 
Tabulations of lake levels referred to some plane of reference 
are in the pubbished reports, but the elevations are Trequently. 
inconsistent, contain some gross errors, and as a whole’ they 
are not fully reliable. ... During the period 1880 to 1899, 
the+ sages ‘apparently received bul dittle (care, “and “au *peincs 
where there existed unstable conditions, the records have but 
little value. Since 1899 the gages have been well cared for, 
self-registering instruments have been installed, and frequent 
eheck Jevels have been run," The Locations of many of the pres— 
ent gages are not the same as those of the earlier ones, “and 
even where the location is the same it has been impossible in 
some cases to recover the old benchmarks. 
At Harbor Beach the records began in 1875. The records before 
1900 are staff gage readings, and there is evidence that the 
gage was not well cared for. There were frequent changes in 
gages and benchmarks, but the majority of the benchmarks have 
been recovered and their elevations determined. The records 
as a whole are probably good, although considerable variation 
between individual years is noted. 
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For Lake Erie, the gage at Cleveland is standard. The 
record goes back to 1860, and although the history of the gage 
dseaustLimes jobscure, the, éleyations are believed to be good. 

At Buffalo a self-registering gage has been maintained since 
1899. The records before this time were staff readings made 
On various gages at different points in the Harbor. No level 
connections between the earlier and the later work can be 
found. 

- At Amherstburg there are records of a self-registering 
gage since 1899, but the observations scatter rather widely due 
PLOUADIV LOnarVverinole reall Aan, twne «river. 

Oswego, at the present time, is accepted as the standard 
pape for “Lake Ontario.” Its history in the earlier years is 
somewhat obscure, but the elevations since 18600, have been 
accepted as correct. Until about 1880 Charlotte was consid- 
ered the standard gage. Several attempts have been made to 
obtain something consistent out of the-records at this point, 
but never with any success. ... The gages and reference 
points have always, settled,: although it is believed that the 
original benchmark on the lighthouse has been stable. 

At Port Dalhousie, actual gage readings and direct levels 
to the benchmark are available for 1075, at the time of the 
transier of precise levels. The records for earlier years are 
referred to depths over the lower sill of the lock, and to 
make them agree with other records on Lake Ontario it is nec- 
essary to assume changes in the elevation of the sill as great 
as two feet, which of course is impossible. 

At Toronto there are records of the water surface eleva- 
tion since 1300, made by the harbor master. It is believed 
that they are all referred to the same benchmark, and that 
the latter has been stable. However, various devices have been 
in use to indicate within the harbor-master's office, the stage 
of the lake, and these have usually been faulty. This has re- 
sulted in errors, which, due to failure to check the accuracy 
of the indicating mechanism have at times been carried through 
several years. As a result the observations scatter very badly, 
but a mean line is probably very near to the truth. 

The gage records at Kingston fall into two groups, 1895- 
1901, and 1908-1919. The earlier records are depths over the 
invert of the dry dock, the elevation of which was determined 
by comparison with Tibbetts Point. The later group is referred 
to benchmarks the elevations of which were determined by 
another water level transfer. So far as is known there is 
no instrumental level connection between the two groups. The 
records of the self-registering gage at Tibbetts Point, main- 
tained since 1901, scatter badly due to local conditions, and 
the determinations of a line through them is unsatisfactory. 
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ERRORS DUE TO GAGE LOCATION 

Two types of gage locations are important in affecting the accuracy 
of lake-level gage readings, which in turn influences the accuracy of 
precise leveling and rates of uplift. The first type of location is the 
relation of the gage site to the local geography and underwater condi- 
tions, and the second type is the relation of the gage site’ to the tarea 
of horizontality of past uplift. 

The effect of local topography and depth of water upon the height 
of the water surface has long been known (e.g., Humboldt, 1849, pp. 309, 
410} Whittlesey, 1059) ‘pps 6, O-10;918, 2; Henry, 10s pe peu 
and the importance of gage site location with respect to the area of 
horizontality of former glacial lake shorelines was clearly explained. im 
Frank B. Taylor (1927) in a paper to the Michigan Academy of Science 
Paper entitled "The Present and Recent Rate of Land-Tilting in the 
Region of the Great Lakes." 

An exampilesor the effect of ‘gage, locatvon=near a rivert@ic 0 Venmes 
the monthly mean elevation sheets for Oswego, New York before 1932 (U. S. 
Lake Survey, n.d.) where a note at the bottom of the elevation sheet said: 

iherrecords have: been carciulilyeke purande anemic Liciike 5 

but do not always indicate the stage of Lake Ontario, as the 

gage readings are influenced by High stages of water in Oswego 

Rivere.'These high stages usually occur in spring: months, and 

give readings trom 0.27.5 tO-Os5 Lt. sOresevem more, -av0ve me 

level. of Lake @nvar Lo. 


The interconnécting ramifications of the erfects of lLake-level gage 


locations and their importance to precise leveling and accurate computation 
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of rates of uplift is very well brought out by the relationships of four 
gages on Lake Ontario. 

Sherman Moore (1941, pp. 2-3) in his gage history of Tibbett's 
PoLnceoec LAred: 


10. The gage record aft Tibbetts Point is of but little 
value, due to local conditions which resulted in failure of 
the gage to give Lake Level. From the first, the gage was in 
the well from which water was drawn for the boilers in the 
fog signal station. The bottom to the southwest of the point 
is very flat and shallow, and to insure a supply of water in 
the well at low lake stages, a channel several hundred feet in 
length to deeper water was cut in the bottom. This was filled 
with large broken rock. With westerly winds, the seas running 
over the flat bottom and beach would hold the water in the 
well above lake level, and the gage would read too high. 

14. The elevations at Tibbett's Point here listed do 
not give true elevations of Lake Ontario within about 0.15 
ft. The discrepancy is a function of the direction and ve- 
locity of the wind. 


The significance of the above quotation may be seen when its con- 
clusions are compared to the following material from the gage histories 
of the Canadian Hydrographic Service (n.d.): 

Elevations at Kingston on 1904 Datum are based on a com- 
parison of float gage readings for 1909 and 1911 to 1915, with 
water surface elevations at Tibbett's Point. 

Elevations at Toronto on 190% Datum are based on compar- 
isons of float gauge readings from 1907/7 to 1909 with water 
surface elevations at Tibbett's Point and from 1917 to 1925 
with water surface elevations at Kingston. 

Elevations at Port Dalhousie on 19043 Datum are based on 
a comparison of float gauge readings from 1914 to 1917 with 
water surface elevations at Kingston. 


If the elevations at Tibbett's Point are "of but little value," and 


if the elevations at Kingston, Ontario are based on the Tibbett's Point 
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readings, then their accuracy must be doubtful. Furthermore, if the 
elevations at Port Dalhousie and Toronto are based on the Kingston float 
gage elevations, they must also be of questionable accuracy. 

This chain of dependent elevations helps to explain why precise lev- 
eling (at its present state of accuracy) cannot be used to determine rates 
of crustal movement in the Great Lakes region (see pp. B7e90)* 

In view of the long history of the effect of topography and water 
depth upon the water-level and upon its recording at water-level gages, 
and especially in view of F. B. Taylor's explanation of the importance 
of gage site location with respect to the area of horizontality, it is 
difficult to understand why the importance of gage location to; the care 
culation of rates of crustal movement was not realized by those inves- 
tigators. who determined rates of uplift for the Great Lakes region. 

It is particularly difficult to understand how B. Gutenberg (1941, 
pp. 740-741) could say: 

Prom the fact that, there is practi ea) ly snevrelative 

change between Calumet Harbor and Milwaukee on Lake Michigan, 

and bevween Port. Stanley and Cleveland one lake Brie jaaieds 

eoncTuded that the zero isobase runs orwicenOruny Oimbiese 

stations, and consequently that the zero assumed for these 

two lakes approximates the absolute zero. This agrees with 

the findings of Taylor (1926) that, at least since the time 

of the Nipissing Great Lakes, which was about 5000 years 

ago according to Antevs (1939), there has been no noticeable 

Uplit to cine vie snes von. 
and then use Calumet Harbor, Milwaukee, Cleveland or Port Stanley as one 


gage of 10 pairs of gages out of a total of 14 gage pairs on Lakes Migh- 


igan-Huron and Erie. It is obvious that’ uplift cannot be measured in an 


ie a 


area where uplift has stopped. Nor can accurate rates of uplift be found 
wiere One gage Of a pair is in an area where uplift has stopped, the 
other gage is in an area of uplift, and the point between the two gages 


at which uplift begins is unknown. 


= LAKES RATES OF CRUSTAL MOVEMENT BY PRECISE LEVELING 

One attempt has been made to determine absolute rates of crustal 
movement in the Great Lakes area by means of precise leveling. The re- 
sults of this study are found in Sherman Moore's (1948) paper "Crustal 
Movement In The Great Lakes Area," pages 702-706 and Plate I. 

As a result of the factors which are discussed in pages 87-90 
of this paper and to several erroneous assumptions by Moore, these rates 
Of crustal movement are not valid. The initial point for the determi- 
nation of elevations for the Great Lakes area is Rensselaer (Greenbush), 
Newercrk, end 10 1S lor Lhis point that the first. erroneous value for 
the rate of crustal movement was computed. If the initial rate of 
crustal movement is incorrect, the values for the other points in the 
leveling net. are probably also incorrect. 

Moore (1948, p. 702) stated: 

Between mean tide, or half tide, at New York, which at 

that point are practically identical, and Rensselaer, near 

Albany on the Hudson River, there are levels in 1857, LES Tale 

and 1934, all run by the U. S. Coast and Geodetic Survey. 

By the levels of 1934, the land at Rensselaer is lower by 

0.849 foot than by the levels of 1877, corresponding to a 

subsidence at a rate of 1.49 feet per 100 years. The ele- 


vation by the levels of 1857 falls only 0.0 foot from a 
line through the other two points. This early determination 
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has now been used, as it was considered less accurabe than tae 

later Jdievels,- bub ats anclusion would have had only a -megiizioie 

erfecuwon pone rate. 

The specious reasoning underlying the determination of this rate of 
subsidence at Rensselaer is revealed by an examination of the U. S. Coast 
and Geodetic Survey precise leveling elevations of benchmark Gristmill at 
Greenbush (Rensselaer), New York, and by the leveling net adjustments 
which were made by the Coast and Geodetic Survey. 

The elevations (Comstock, 1876, p. 71; U. S. Deep Waterways Comm., 


LGO{s pp.) (O-(siHayrord, 1900, pps 449, 540; Hayford. 1905,, pp too. 


2095 (0s. Dos DOW eG, HL O14 pl O> JS arenas. Collows: 


Year Klevation in Heat 
ane Sse eo 
MeN We. 720 
18942 13.645 
1899 Vey (il 
1902 14.873 
19043 13.063 
1907 14...005 
1912 13.865 
1929 13.618 
19342 13.845 


“Level lines run Dy athe Usps Const 
and sGECdetilesourvicy. 


The most recent determinations of bench mark Gristmill based on the 


1929 adjustment (Gossett, 1961, letter) are as follows: 


13.619 feet (1902 leveling) 
132555 feet (1934 leveling) 
13.501 feet (1955 leveling) 


The changes of elevation in the above paragraphs must be examined 


with the following facts in mind (see p. 84): (a) First order leveling 
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‘by the U. S. Coast and Geodetic Survey began in 1878; therefore the 
levels of 1857 and 1877 were not of first-order precision; (b) the 
first adjustment to the level net occurred in 1899, partial adjustments 
were made in 1903, 1907 and 1912, and a complete adjustment was made in 
1929 (Rappleye, 1948a, p. 1). For these reasons the logical choice of 
the early elevation for Gristmill (bearing in mind the precision neces- 
sary to calculate rates of crustal movement) would be the elevation of 
13.577 feet resulting from the first adjustment in 1899. 

The explanation for the changes in elevations at Gristmill lies 
not only in the fact that better equipment and techniques were used in 
the later levelings (1894, 1934, 1955), but also to the fact that new 
leveling lines and better determinations of sea level were introduced 
into the precise leveling net with each adjustment or readjustment to 
the net. As was stated in the U. S. Deep Waterways Commission Report 
or oop (py 7 

In 1894 the Coast and Geodetic Survey ran a line of pre- 

cise levels along the Hudson River, starting from the bench 

mark at Dobb's Ferry. The superintendent gives the elevation 

of the bench mark on the gristmill as 13.645 feet, with the 

following note: "The difference between the above (13.645) 

and any former results is probably due to the more perfect 

determination of tidal level than to any ‘other cause. 

In other words, the adjusted elevations are not necessarily the re- 
sult of a subsidence or uplift at the bench mark, but, instead, represent 
the introduction of new information into the precise leveling net. The 


bench mark shifting shown by the elevations is almost purely a "paper" 


movement. 
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If the elevations for the years 1877 and 1894 had been used, the 
subsidence would have been 1.083 feet L/ years, Or) 0.41 feet/100 years. 
or, if the years 1899 and 1902 had been used, the subsidence would have 
been 0.296 foot/3 years Ox OO f feet/100 years. 

In his section on "Reduction to Sea Level," S. Moore (pp. 704-705) 


States: 


For, & correlated picture of the movement, as: a-whole.y che 
observations must be reduced to a common datum. ‘The only 
practical datum for this purpose seems to be sea level. - if 
one admits! a changing sea level there is no means of deter- 
mining whether the movement is uplift or subsidence. Great 
Variations in the relative elevation of lland and ocean level 
during geologic time seems well established, but it seems im- 
probable that there has been any appreciable progressive change 
Lithe «volume, of, oceanic water vim the ast) LOO-veare-. 


As was discussed in pages 45-48, 56-58, of this paper, there is a 
world-wide change in sea level of ¢ 0.36 foot/100 years,. as. well acre 
change in sea level at New York City (the starting point for Great Lakes 
leveling lines) corresponding to 0.78 foot/100 years. ‘One Gan, admins o 
changing sea: level and by measuring the change and correcting for! it, 
decide, whether @ movement is subsidence or uplift.) The heignt of the 


ocean level in the geologic past has little or mo bearing on) the: problem: 


Significance of Previous Determinations 


An investigation of the gage histories of U. S. Lake Survey gages 
(Sherman Moore, 1949-1944) and of the gage histories of the Canadian 


Hydrographic Service emphasizes the many opportunities, particularly in 
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the earlier measurements, for error to enter into the determination of 
lake-level elevations due to equipment limitations, shifting bench marks, 
incorrect or lost records, observer mistakes, gage locations, etc. 

Rates of crustal movement in the Great Lakes region are determined 
from the gage differences of pairs of gages. The gage differences in- 
clude not only the changes due to land uplift but also the residual ef- 
fects of wind set-up, the barometric pressure effect, the changing dif- 
ferences sitieeciche wmplitudes,. tidal Giiterences, Observer error, in- 
strument error,-bench‘mark changes, etc. While it is very probable that 
many of these effects are small in magnitude, or that many of them are 
opposite in sign and thus compensate each other, nevertheless the errors 
caused by these effects must be recognized and removed, or corrected, so 
that the true magnitude of the uplift may be approached as closely as 
possible. Because the magnitude of land uplift is so small, the com- 
bination of the various errors which are preserved in the lake-level 
gage records tend to completely mask its effect and make its detection 
with present methods difficult, if not impossible. 

The change in the character of a plot of gage differences over a 
ate of 95-96 years is illustrated graphically by Fig. 5 which shows 
gage differences between gage pairs on Lake Ontario and Lake Erie. The 
great decrease in the fluctuations of the gage differences as operator 


and instrument error was reduced by means of improved equipment, gage 


locations and gaging techniques is readily seen. 


(A) 
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Fig. 5. Gage differences on Lake Ontario (A. Toronto minus Port 
Dalhousie) and Lake Erie (B. Cleveland minus Port Colborne) show- 


ing decrease in range of fluctuations with improved equipment, 
locations, and techniques’. 
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The proportions of various types of error which distort true lake- 
level values changed throughout the period of record taking. In the 
earlier days of lake gaging, say before 1920, bench mark stability and 
instrument and operator error predominated, whereas the systematic errors 
Poiecd by meteorological effects, tides, harbor effects, etc., were. rel- 
atively small. After improved equipment, techniques, etc., removed the 
grosser aspects of instrument and operator error, the importance of mete- 
orological effects, etc., grew proportionately. 

At the present stage of lake-level gaging the greatest deviations 
from true lake-level elevations (on each lake—not referred to sea level) 
are due to meteorological effects and the effect of gage location. Be- 
fore worth-while determinations of crustal movement can be made, correc- 
tions for these effects, as well as for the remaining instrument and 
Operator errors, must be made. \ 

‘the following section discussing modern rates of land uplift around 
Lake Erie points up the difficulties which arise when all of the aspects 
of a problem are not viewed in a coordinated manner—in this instance 


the failure to examine and analyze all aspects of the problem caused a 


case of "mistaken identity" to arise. 


LAKE ERIE 
An examination of Plate I reveals that Lake Erie and its shores have 
been completely within the area of horizontality since Nipissing time, a 


period of approximately 3000 years. Despite the fact that no uplift has 
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oceurred for 3000 years, all of the determinations of rates of crusvam 
movement (with the exception of Horton and Grunsky, 1927, p. 32) have 
found uplift around Lake Erie—the rates varying from 0.009 foot/100 
miles/100 years (Moore, 1948) between Cleveland and Port Stanley (Guten- 
berg, 1941, finds zero uplift for the same pair) to 1.04 feet/100 miles/ 
100 years (Moore, 1922) between Cleveland and Amherstberg (see Table 3). 

As postglacial uplift probably ceased about 4000 years ago, the 
rates of the quantity called "crustal movement" must in actuality be 
rates of some other quantity. Because Great Lakes water-level gage rec- 
ords’ are not corrected for the’ factors which were discussed in Parca 
of this paper (pp. 41-76), the previously determined "rates of crustal 
movement" measured the change in the net difference of the accumulated 
effects and errors (chiefly meteorological) which are incorporated in 
the gage records. The "rates of crustal movement" would be more accu- 
rately labeled "average net setup" for the period of record of a given 
pair of gages. 

The same method is used to determine rates of crustal movement in 
the entire Great Lakes region, and as this method has yielded rates en 
crustal movement up to 1/04 feet/100 miles/100 years in an area where 
uplift had ceased several thousand years ago, and as this rate exceeds 
all of the Nipissing rates (Table 1) and all but two of the computed 
modern rates (Table 4), it is reasonable to assume that the value of the 
rates’ of modern erusval movement , even the existence et modernvcris tar 


movement, around the Great Lakes is in doubt. ) 


V. COMPUTATION AND CORRELATION OF METHOROLOGICAL EFFECTS 
AND GAGE DIFFERENCES 

As has been pointed out. previously, gage differences are the basis 
for computing rates of crustal movement and are also the basis for com- 
puting set-up between pairs of gages. The justification for using the 
seme Quantity as the basis for two entirely different phenomena is the 
assumption that winds over the Great Lakes in the summer months blow in 
opposing directions so that "the mean surface of each lake is level" 
(Comstock, 1882, p. 592) and, therefore, that set-updoes not exist to 
any significant degree; consequently gage differences represent movement 
of the land between the gages of each pair. 

If the winds over the Great Lakes do not neutralize each other's 
effects On» the water surtace during the summer months, i.e., if vector 
resultant winds from one quadrant predominate; then the wind from the 
prevailing direction. will exert its infiuence and net set-up will occur 
inethe direction of the prevailing wind. If net seteupdoes otcur, a plot 
of the gage differences over a number of years will measure the average 
net setsupfor that time, not the rates of crustal movement. 

The crucial point is whether summer winds can produce a predominant 
summer vector resultant wind whose setup can be correlated with the mag- 
nitudes of the gage differences. If such a correlation can be shown, the 
rates of crustal movement which have been calculated heretofore are not 


valid. 


Alb, 
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Previous Statements of Westerly Prevailing Winds 
Mark W. Harrington's (1895) comprehensive study Surface Currents of 
the Great Lakes included two tables of wind directions over the Great 
Lakes. In writing of the most frequent winds on the Great Lakes (as 
shown by compilation of tri-daily readings from 1871-1888) Harrington 
(p. vi) stated: 
If we take the months May to September, inclusive, the 
numbers are Wa, sein terol, OF OM. 4. 1 O'S Ue lOsso heme 
W., 13; NW., 7. This is 88 for these months, of which 43 per 
cent are SW; 15 percent W; and 8 percent NW.; or 56 percent 
from a westerly direction. 
The prevailing character of the westerly winds at the 
lake station is shown still more clearly in the resultant 
wind directions (Table II): 


The wind directions for the following stations were taken from 


Harrington's Tabse Tip. ie): 


June July Aug. Sept. Annual 
Buffalo Soe Wee 6.50 We USO OW SUR e ya ene tans 
Duluth Weis Hee’ Nee We NN SSW Nee Se We lar hee ane 
Oswego S01 Wea, SafoaWs S755 Woe! 6.20 Wen es ome 
Rochester! “S.G0lW. 670 We) NGL (O We = ROW, 8 ons ene 
Toledo SHOW SOU We S67 We Seago oN oon aus 


William T. Blunt's "Effects of Gales on Lake Erie" (1897, pp. 157- 
158, diagram no. 2) also contained tables and a graphic representation of 
winds as recorded at Toledo, Ohio (1891-95). Blunt's tables showed that 
winds from the southwest and west exceeded those from the northeast and 
east during eleven months of the year (including all of the summer months). 
Blunt (p. 158) stated that, "... the resultant movement and direction in 


the average months are most decidedly from the westward. ..." and "this 
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naturally tends to keep the mean level of the east end of the lake higher 
than that of the west end." 

Wind directions of all stations, except Cleveland, around the Amer- 
ican shore of Lake Erie (based on 15 years of observations), during the 
months of June, July, August and September blew from the southwest (from 
S 13° W to S 73° W) (Henry, 1902, p. 11). Regarding the wind directions 
at Cleveland, Ohio, A. J. Henry (p. 11) declared, "The unusually large 
number of southeasterly winds at this station is not clearly understood, 
unless as suggested by local forecast official Kenealy, of Cleveland, 
they are due to land and lake breezes." 

A modern summary of Great Lakes wind velocities (USWB, 1959, pp. 
13-14) confirms the earlier observations that the prevailing summer winds 
are from the southwest on Lakes Erie and Ontario and are from the westerly 
Quadrant on Lakes Michigan, Huron and Superior. 

As the preceding excerpts illustrate, it has been known for 60 or 
more years that the predominant wind direction for the Great Lakes is 
westerly; in addition, the prevailing southwesterly direction of summer 
winds on Lake Erie has also been recognized since the 1890's. 

Because it has been known at least since the early nineteenth cen- 
tury (Dwight, 1822, p. 96; Hall, 1843, p. 200) that southwest winds on 
Lake Erie cause a rise of water level at the eastern end of the lake, 
and because it has also been known, since the 1890's (if not before), 
that southwest winds predominate during the summer months, it is sur- 


prising that the assumption that the summer season mean lake surface is 
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a level surface has been questioned only by William T. Blunt (see p. 


102). 


Present sLrudy 


A review of the material in the preceding sections prompted the de- 
cision to test the validity of the assumption which underlies water-ley-= 
eling, as well as. the determination of rates of crustal movement, te 
test to be by correlating the gage differences of a pair of gages Onjume 
Great Lakes with summer wind velocities on the same lake. 

For a nmumber of reasons, Lake frie was an ideal localrioneror sie 
correlation study; the more important reasons being: 

(a) Lake Erie is oriented so as to be almost parallel with the pre- 
vailing winds of the summer months. 

(b) Water-level gages (at Toledo and Buffalo) are located at the 
ends. of the lake.: and. the direction of 4a line from Woledo ito sur to lowes 
within 3° of the direction of the lake's axis. 

(c) First-order U. S. Weather Bureau stations record wind velocities 
relatively near (9-18 miles) the gage sites and at two other points where 
the lake is divided approximately into thirds (Cleveland and Erie). 

(d) Lake Erie is shallow (average depth c 58 feet) and the gage 
Sites are at the converging ends of the lake; both factors increase the 
magnitude of set-up, which in turn, amplifies the relationship which 


exists between wind velocities andset-up (gage differences ). 
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(e) Lake Erie is located wholly within the Nipissing area of hor- 
izontality; therefore postglacial uplift could not have occurred during 
modern times, and the rates of crustal movement which have been calcu- 
lated must represent some other quantity. 

Similarly the gage differences of gage pairs in the other Great 
Lakes could also be correlated with wind velocities. The gage pairs 
chosen should meet three requirements: 

(1) The direction between the gage sites of the pair of gages must 
parallel, at least roughly, the direction of postglacial uplift (i.e., 
one gage should be about N 70° E from the other). 

(2) Both gages of a gage pair must be north of the Nipissing zero 
isobase, and 

(4) There must be nearby weather stations which record wind veloc- 
ities and whose records are available. 

These requirements restrict testable additional gages to two pairs. 
One pair (Duluth-Port Arthur) is located on Lake Superior, and the other 
pair (Toronto-Kingston) is on Lake Ontario. Although these gage pairs 
fill the requirements better than any of the other available gage pairs, 
serious deficiences remain which greatly limit their usefulness; however 
the deficiencies serve to point out the need for more complete meteorolog- 
ical and lake-level data on the Great Lakes. 

Owing to gage site locations as well as to the fact that the average 


angle between the summer wind direction (1950-59) and the direction of the 
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Duluth-Port Arthur axis was 41°, it is expected that the correlation 
between gage differences and wind velocities would be low. An examina- 
tion of current charts of Lake Superior (Harrington, 1895, Lake Superior 
ChertseMiilear, Lose. Hie. 5) reveals circumstances which would tend to 
support this expectation. The current from the area of Thunder Bay 
(Port Arthur) Ontario moves southwestward along the northwest shore 

of Lake Superior to an area somewhere between Grand Marais, Minnesota, 
and Devil's Island,..Wisconsin,“where It curves. ‘to the east. The secur. 
rent flowing along the northwest shore from Duluth moves northeastward 
to an area south-southwest of Grand Marais where it also curves to the 
east, or continues to curve around until it is flowing to the south- 
southwest. Thus the situation exists where the current direction ay 
one gage is 180° from the current direction at the other gage; as setup 
is dependent upon the wind-drift current (which in this case is flowing 
in opposite directions) it is to be expected that a very low or non- 
existent correlation would be found. 

The Toronto-Kingston gage pair on Lake Ontario meets the first two 
requirements satisfactorily, but proves to be inadequate on the third— 
that of nearby wind and barometric pressure observations. Wind and bar- 
ometric pressures are recorded at the Class I weather station at Toronto 
at the western end of the lake, but records of wind velocities and bar- 
ometric pressures are not available for the eastern end of Lake Ontario. 


The nearest available wind and barometric pressure observations were 
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taken at Trenton, Ontario, which is located only two-thirds of the dis- 
tence from Toronto to Kingston. “A’gap in the Trenton records for the 
years 1954 and 1954 was filled in by wind observations taken at Main Duck 
Island. Main Duck Island has an excellent location in open water at the 
eastern end of the Take, but the observations were recorded by a radio 
beacon station which is a class c weather station (the Canadian Mete- 
erological Division stations are classed as I, II, III and ¢). In 
addition the record at Trenton is available only to 1955. 

Because wind and barometric pressure information is not represent- 
ative of the eastern end of Lake Ontario (observations were made about 
54 miles from Kingston), the correlation of gage differences between 


Toronto and Kingston and wind velocities should not be significant. 


Lake Erie Winds 


VECTOR WINDS 

At the beginning of the study it was expected that mean wind veloc- 
ities for the summer months could be obtained for the locations and years 
under investigation directly from the U. S. Weather Bureau Local Clima- 
tological Data and the Canadian Meteorological Branch Monthly Meteoro- 
logical Summaries. However it soon became apparent that prevailing wind 
directions and average speeds from these compilations were not suitable 
for calculating set-up or for correlating with gage differences. 


The wind speeds and directions reported in the Local Climatological 
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Data are determined as follows: 


The prevailing wind directions for the month is the direction 
which has the greatest total number of hourly occurrences arrived 

at by the summation of hourly observations in Table B, Wind Direc- 

tion and Speed Occurrences, published in the Local Climatological 

Data Supplement ... (Fox, 1960). 

Average Wind Speed... at the foot of the column, enter the 

sum of the daily average hourly speeds, and the average of these 

speeds as obtained by dividing the sum by the number of days in 

the month. ... (U. S. Weather Bureau, 1959). 

Wind speeds and directions arrived at in this manner are suitable for 
work in climatology, for evaporation and cooling studies, etce., but noe 
for work which is concerned with transport, in this study with the trans- 
port of water. The transport of water by a wind-drift current is de- 
pendent upon the wind stress, which in turn depends upon vector winds. 

Because vector wind speeds and directions are not published (or, ex- 
cept in certain studies, determined) for the Great Lakes, it was neces- 
sary to compute the vector winds for Lakes Erie, Ontario and Superior. 
The length of time required to plot the wind vectors restricted the 
time period of correlation to ten years (1950-59) for Lake Erie and Lake 
Superior: The lack of wind information imposed the eight year period 
(1948-55) for Lake Ontario. 

The procedure which was used to plot the vector winds is given in 
pages 138-141. Tables 7-66, (pp. 180-227 of Appendix II) contain the 
daily, monthly and summer season vector winds for: (a) Lake Erie— 


Toledo, Ohio, and Buffalo, New York, (b) Lake Ontario—Toronto, Ontario, 


and Trenton, Ontario, and (c) Lake Superior—Duluth, Minnesota, and Fort 
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William/Port Arthur, Ontario. 

Over-water vector resultant winds computed for the summer season 
(June, July, August and September) for Lake Erie (1950-59) and Lake 
Ontario (1948-55) are summarized in Table 5. ‘The summer season vector 
resultant winds for Lake Superior are land station winds because an 
over-water:land wind ratio has not been determined for Lake Superior. 

All studies have indicated that the prevailing summer wind direc- 
tions for Lake Erie and Lake Ontario are from the southwestern quadrant, 
and that the summer prevailing wind directions for the other Great Lakes 
are from the western quadrant. Winds blowing over the water surface 
from the southwest cause a net set-up to occur—the water surface being 
depressed on the windward shore and piled up on the lee shore (see Fig. 
4b). Theset-upis a tilt of the lake surface—the tilt being up toward 
the northeast. 

Postglacial crustal movement in the Great Lakes region has con- 
sisted of an upwarping to the northeast; therefore when investigators 
(seeking to measure postglacial uplift by means of uncorrected lake- 
level gage readings compared with the presumed "level" summer lake sur- 
face) plotted gage differences for pairs of gages and found a northeast 
tilting, they apparently assumed that the change in gage differences was 
produced by crustal movement. Because investigators thought that the 
lake surface was level, whereas it was actually tilted upward to the 
northeast, the northeastern gages of the pairs appeared to be uplifted 


by the amount of the net set-up. 
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TABLE 5 


A, LAKE ERIE OVER-WATER WINDS (JUNE-SEPTEMBER ) 
MEAN VECTOR RESULTANT WIND FOR TOLEDO-BUFFALO 
Summer Season Over-Water Vector Resultant Wind 


Year ey COM. C mph _ eevee iM ee = BOT 
1950> 21h° (ssw) 265 1955 BOT?) SW) 
Osi! Pele Sw) Bike: 1956 225° = Caw) 
1952 e1h° (sw) 576 1957° 238° (wSw) 
1953 223° (SW) 505 1958 eho°  (WSw) 
1954 a we) ae is Te 237° = (WSW 


Bea on three months ie. August, eee 
eBased on three months (June, July, August) 


B. LAKE ONTARIO OVER-WATER WINDS ( JUNE-SEPTEMBER ) 
_MEAN VECTOR RESULTANT WIND FOR TORONTO- TRENTON _ 


ear foe mph oie Year ; From mas 
1948 2080" (WW) 64 1952 ehe° (WSW) 
1949 253° (WSW) 460 eb. 251° (WSW) 
1950 260° (WwW) 598 19548 270° (Ww) 
19pL one oy Boe 1955 246° (WSW) 


2Observations from Main Duck Island 1953-1954. 


Year 


1950 
1951 
1952 
aes 
1954 


C, ' LAKE SUPERIOR LAND SPATION WINDS ( JUNE-SEPTEMBER ) 
MEAN VECTOR RESULTANT WIND FOR DULUTH=PORT ARTHUR 
Summer Season Land CRESS Vector Resultant Wind _ 


oe iG Onl 
2627 (WwW) 
278° (WwW) 
26° (WSW) 
262° (Ww) 


cou 


(NNE ) 


_mph 


290 
146 
225 
216 


80 


Year 


LED 
1956 
1 ( 
1958 


197 


_ From ie 
266°- (W) 
337° (NNW) 
OT Oe ke CMe) 
266° (WwW) 


anbee 


(SW) 


_mph 


304 
582 
ee 
(59 


NS 


mph 
638 
476 
410 


Bey, 


mph 


1LO9 

86 
LS 
420 


ihe 
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When gage differences were plotted for a number of years and a best- 
fit curve drawn, the slope of the curve was called the "rate of crustal 
movement;" in fact, the slope of the curve represented, for the most part, 
the average net set-up for the time period being considered. 

It must be emphasized that the quantities which were determined and 
used in the following sections, i.e., gage differences, barometric pres- 
sures, and effective wind velocities, are net quantities; that is, they 
are residua of winds, barometric pressures, and wind and barometric pres- 
sure effects which occurred over the four month summer period. As quan- 
tities which do not represent phenomena that were actually observed, but 
rather averages of actual conditions, they should represent a quasi-steady 
state condition, and the relationships deduced from their study may or may 


not apply in detail to momentary, hourly or daily conditions. 


EFFECTIVE WINDS 

Summer season over-water mean monthly vector winds (summer mean 
monthly winds) cannot be used directly in the computation of set-up, or 
in correlating gage differences with set-up—they must first be converted 
to effective winds, i.e., winds which represent that portion of the over- 

water wind which causes the observed setup at the gage sites. 

Winds vary from 0% to 100% in effectiveness—100% effective winds 
are those winds which cause the observed gage differences (corrected for 
barometric pressure effects). Zero per cent effective winds are winds 
which blow at right angles to the 100% effective winds. Ideally, effec- 


tive wind velocities when squared and plotted against corrected gage 


Ab) 


differences would yield a + 1.0 correlation coefficient. 

Effective winds are usually found by resolving vector resultant 
winds into their components; in this case, effective. winds are functions 
of the cosines of the angles between the resultant wind directions and the 
direction of the lake axis (deviation angles). However, effective winds 
are defined as winds which cause the observed setup, and not as the simple 
components of resultant winds. Because resultant winds must exert their 
influence across the air-water interface, and because they blow over a 
eunved surface and not-a flat plane, effective winds may now be ecm. 
funetions of the resultant winds, instead they may bear some other rela= 
tionship to resultant winds. With this possibility in mind) errecuiye 
winds were computed in two ways: (1) one computation was made using 
the cosines of the deviation angles from 0-90 degrees; and (2) the 
other computation was made using a linear relationship, i.e., winds were 
considered to decrease 1.11% in effectiveness for each degree that the 
resultant wind directions deviated from the direction of a line con- 
necting the two gages under investigation (see pp. 139-140). 

The test to determine the correct method is by comparing the degree 
of correlation (shown by the correlation coefficient r) between the ef- 
fective wind velocities squared and the observed gage differences cor- 
rected for barometric préssure effect. If winds from one particular direc= 
tion are most effective in creating the set-up, the correlation coefficient 
for this direction and the set-up should be at a maximum. Thus the values 


of the correlation coefficients should increase as the wind directions 
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approach the 100% effective wind direction, be at a maximum at the 100% 
effective wind direction, and decrease as this direction is exceeded. 
If winds blow from a direction opposite to those which cause the quasi- 
steady state setup, i.e., they blow "down slope," they are called neg- 
ative effective winds. For example, a wind hae ULowe In -a direécvion 
180° from the direction of the 100% effective wind would be a -100% 


effective wind. 


Direction of Effective Wind 

Previous investigators (Hellstrém, 1941, pp. 17-18; Keulegan, 1953, 
pp- 102-104; Harris, 1954, p-. 38; Gillies, 1960, p. 37) have considered 
the 100% effective wind as blowing parallel to the lake axis, i.e., on 
Lake Erie the 100% effective wind would be from 248°. The choice of the 
lake axis direction as the 100% effective wind direction presumably fol- 
lows the reasoning of V. W. Ekman's theory of currents which postulates 
that the rotation of the earth does not deflect surface currents in shal- 
low water, and that, consequently, wind-drift currents flow in the direc- 
tion of the wind. In addition, Ekman's theory proposes that the slope of 
the water surface will always be in the direction of the wind in water 
of any depth. Because these conclusions do not agree with the results 
of this investigation, Ekman's theory will be examined briefly in the 
light of the observations which prompted its formation. 

V. Walfrid Ekman originated his now classic theory of ocean currents 


in 1902 at the suggestion of Fridtjof Nansen who discovered during the 
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drifting of the ice pack which held his ship, the "Fram, |: that thesice 
drift of a given wind deviated to the right of the wind direction. fFol- 
lowing the original Norwegian publication of the theory in 1902, Ekman 
expanded his theory which was then published in 1905 as "On the Influence 
of the Earth's Rotation on Ocean Currents." 

The theory represents conclusions deduced from a mathematical model 
of ocean currents incorporating several simplifying assumptions; these 
assumptions cause the model to differ from conditions as they exist in 
feture (see pp. 54-56). Ekman (1905) states several conclusions which 
are pertinent to this discussion as follows: 


Equations (5) then show that in the northern hemisphere 


ee 


In the southern hemisphere it is directed 45° to the left... 
(Goma) 


The above-mentioned result, according to which the sur- 
face-current's deflection from the wind-direction, is invari- 
ably 45°, seems rather strange; one would indeed expect the 
earth's rotation, to have -less- influence on, the curren tsicetae 
smaller its vertical component wo sin d fi, LOI 


The angle q between the wind and the surface-current, is 
not exactly 45°, when the depth is finite. ... and the angle 
of deflection q@ consequently depends on the ratio between the 
depth of the «sea d. and the Depth of Wind-Curren’ Daa Ee a/D 
iso smal tracui on, Of 1S small and sone -currenvesocewmearty, 
in thesdimectionsof the ewind.. \,As the depth inereases aces 
alternately smaller and greater than 45°. Thus for instance 
Qi. 214, Su ford 5=40,25.505 10. = 45. showed) =O), 54D aca tla 
ford, = 705 (> Dende. 45° ford =D. When dis greater 
than D, the deviations from the mean value a@ = 45° are quite 
insignificant, and the motion takes place almost exactly as 
On the deepasedu. (peel sh). 


The arrows represented without shaft-feathers give the 
direction of the slope; it is remarkable how nearly this direc- 
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tion follows the wind's direction (common for the whole plate) 
whatever be the depth of water. This shows clearly that the 
earth's rotation has no considerable deflecting influence on 


te a, Ee -tC—. 
Oe oO - 


deviate from the wind's direction). Its influence on the 
absolute magnitude of the mounting up is also found to be 


rather moderate, its effect being to diminish the inclination 

of the water-surface in the ratio 0,98 if d = 0,5D, in the 

ratio 0,77 1f d= 1,25 D, 0,f/l1 if d = 2,5 D, and exactly 

By eet iao ts dinginite ipres |) < 

The above quotations represent the chief deductions relating to the 
angle between the wind direction and the surface wind-drift current, and 
between the wind direction and direction of water surface slope. 

Although Nansen's original observations of the angle of deviation 
between wind direction and ice-drift are given as 20°-40° to the right 
of the wind (e.g., Ekman, 1905, p. 2; Sverdrup, et al., 1942, p. 492), 
an examination of Nansen's (1902; 1904) summaries of the original data 
reveals that the range of the deviation angles was not that uniform. 
Nansen's observations taken aboard the "Fram" from 1893-1896 are sum- 
marized in the histograms of Figs. 6 and 7 (the angles of deviation were 
compiled from Nansen, 1902, pp. 366-67, Table 9, Column 13). 

The histogram of Fig. 6 shows the distribution of the frequencies 
of occurrence of the angles between the wind resultant and the direction 
of ice-drift (after correction for the permanent current). The majority 
of the angles are in the 20°-40° range, but wide variances exist, e.g., 
the maximum angle was 80° to the right of the wind, the minimum angle was 
63° to the left of the wind, and mean angle was 28° (only 62% of Ekman's 


theoretical value). 
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FREQUENCY 

20 

18 

16 

14 Maximum angle = +80° 
Minimum angle = -63° 

12 Mean angle = eeee 

10 - = Left of wind resultant 


+ = Right of wind resultant 


ine) 


= x @ 
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Angle between wind resultant and direction of ice drift 


Fig. 6. Frequency of wind drift—wind resultant angles measured 
during drift of the "Fram," 1893-1896. 


FREQUENCY A. SHALLOW WATER (< l00 m.) 


B. DEEP WATER ()100m.) 


OnNDDWO 


GSS< 


MO tn OG Oe WoGye Ore ten 
Angle between wind resultant and direction of ice drift 


Fig. (. Comparison of frequencies of drift deviation angles for 
soundings in: A. shallow (<100 m) and B. deep (>100 m) water. 
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Sounding information from Nansen's (1904, pp. 9-11) paper "The 
Bathymetrical Features of the North Polar Seas" was combined with wind- 
drift information from the 1902 paper to produce the histogram of Fig. 
({—which represents the distribution of the angles of deviation of the 
ice-drift from the wind direction on the days when soundings were taken 
(56 soundings and deviation angles can be compared). The deep water 
histogram is the more interesting of the two histograms in that it shows 
that 40 of the 47 angles were less than 35° to the right, four angles 
were 59° to the right, and none of the angles fell in the 40.5°-55.5° 
range, whereas according to Ekman (1905, p. 10) the angle "is invariably 
45°." The compilation also includes instances of deviation angles bear- 
ing tO, une lett, of the wind direction instead of to the right as required 
by the theory. 

Nansen (1902, p. 378) gave the angle between wind drift and wind 
resultant as: (a) 26° from November 23, 1893, to November 23, 1894, 
(b) 34° from November 24, 1894, to November 28, 1895, and (c) 23° from 
November 28, 1895, to November 27, 1896. Other observations of the 
Wind-drift direction angles confirm this lower range of values. For 
example, according to Sverdrup, et al. (1942, pp. 623, 666), Brennecke 
found the drift of the "Deutschland" in the Antarctic Weddell Sea to be 
34° on the average, and Sverdrup reported the ice-drift over the North 
Siberian Shelf to average 44° to the right of the wind. In addition, 


G. E. Hutchinson (1957, p. 268) stated that R. Witting in 1909 measured 
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(109 observations) the angle of current deviation from the wind direction 
in water only 9 meters deep and found the deviation to be 33° to the wea. 
In each of these cases an explanation was advanced as to why the ob- 
served angles differed from the theoretical angle rather than why the 
theoretical value failed to conform to conditions as they were measured 
in nature. When empirical observations reveal a consistent difference 
between their values and the values obtained from a mathematical model 
of the phenomenon, it would appear that the model fails to account for 
some relationship which exists in nature,;and for that reason it should 
be re-evaluated and adjusted so as to conform to natural conditions. 
Ekman's statement that the slope of the water surface is always in 
the direction of the wind, rather than at some angle to the right of the 
wind, runs counter to intuitive reasoning of the subject. Because the 
surface slope is caused by the wind-drift current which flows at some 
angle to the right of the wind (in the northern hemisphere), it would 
seem that the surface slope should also be at a maximum in a direction 
at some angle to the right of the wind. Ina study "The Effect of 
Steady Winds on Sea Level at Atlantic City," A. R. Miller (1957, p. 30) 
found that a nomogram relating the factors which cause departures from 
mean sea level indicated that a deviation angle of 20° gave the greatest 
slope if the assumption were made that the greatest departure takes place 
when net transport is normal to the coast-line. Miller (p. 30) stated 


that: 
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If the net transport represents surface drift and the 

gradient wind is replaced by geostrophic wind, the empirical 

angle of 20° is comparable to the computed angles (do - 2, = 

19° to 28°) between surface drift and gradient wind. 

If Ekman's theory is followed, the direction of the effective wind 
would parallel the direction of a line connecting the two gages which 
record the water-surface levels; however, if, as suggested by Miller's 
nomogram, the slope occurs to the right of the wind, the effective 
wind would be from some direction to the right of a line connecting the 
gage sites. 

Lake Erie, for reasons of geographical location, relatively uni- 
form shallowness, wind measurements (first order weather stations around 
the lake) and water-surface measurements (levels taken by modern instru- 
ments and procedures ), presents an excellent situation for testing Ekman's 
hypothesis of water-surface slope direction in shallow water. 

As shown in Fig. 8, the direction of a line between the gage sites 
at Toledo and Buffalo is N 71° EH. If Ekman were correct with regard to 
to slope direction, an effective wind blowing from S 71° W (251°) would 
give the best correlation with the observed gage differences corrected 
for barometric pressure effect. On the other hand, if the greatest slope 
occurs to the right of the wind direction, the wind showing the best cor- 
relation with the set-up (in this situation where the slope direction is 
fixed) should come from some direction to the rieht of the Toledo-Buffalo 


axis. 
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TOLEDO-BUFFALO AXIS \® 
NSZICE: ea) 
<j 


Ped 
we 
_ » 
TOLEDO —e wn / 
// ew 50 MILES 
i eee. a 
‘ee 


100% Effective wind considered as coming from: 
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Fig. &. Lake Erie effective wind directions. Rela- 

tionship of wind directions to 100% effectiveness in 

creating observed setup at Toledo and Buffalo. 
VECTOR WIND AND EFFECTIVE WIND PLOTTING PROCEDURE 

The procedure outlined below was used to convert land station wind 

directions (prevailing) and speeds (average) to vector quantities. Wind 
directions and speeds for American stations were taken from the U. 5S. 
Weather Bureau's Local Climatological Data Supplements and those for 
Canadian stations were obtained from Meteorological Division's Monthly 
Meteorological Summary. 


The plotting procedure for vector winds is as follows: 
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(a) The four synoptic hour observed winds (or winds recorded 
closest to the hours of 0100; 0700; 1300; 1900) were plotted as vectors 
(scale of 1 in. = 5 mph); their resultant, when divided by four, gave 
the daily mean vector wind at the land station concerned. 

(b) In order to obtain over-water vector winds for the lakes under 
investigation, the daily land station mean vector winds were converted 
to over-water vector winds by applying empirically derived over-water:; 
land ratios. These ratios which have been determined only for Lake Erie 
(Hunt, 1958, pp. 28, 30) and Lake Ontario (Bruce and Rodgers, 1959, pp. 
10-11) are as follows: 

(i) Toledo: —Vyater/Viang = 1-59 for June-September, 1950-59, with the 
exception of June, 1950, 1952-54, and July, 1955, when V,,/Vi = 1.13; 
(46) Buffalo: —V,,/Vi = 1.13 for June-September, 1950-59, with the ex- 
ception of June, 1950, when V/V = 0.90, and September, 1956, when 
Wi = labo: 

(iii) Lake Ontario: —(Bruce and Rodgers "spring") Vy/Vi1 = 1.60 for June 
and July. 

(iv) Lake Ontario:—(Bruce and Rodgers "fall") Vw/Vi = 1.90 for August 
and September. 

(c) The 30 (or 31) daily over-water mean vector winds were plotted 
on a scale of 1 in. = 10 mph to secure the monthly over-water vector re- 


sultant winds. Monthly mean daily vector winds were found by dividing 


the value of the monthly over-water vector resultant winds by the number 


of days in the month. 
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(d) The number of monthly over-water resultant winds corresponding 
to the number of stations used were plotted (scale 1 in.=50 mph) as a 
progressive vector plot in which the June wind velocity of one station 
was drawn to scale, being followed by the June wind velocity at the other 
station; the wind velocities for July, August, and September were plotu- 
ted similarly. The vector sum was divided by the number of stations (two) 
to get the summer season over-water vector resultant wind [to be called 
the summer resultant wind ] for’ the; lake Under’ diseussion. 

(e) The summer resultant wind divided by four gave the summer 
season over-water mean monthly vector wind [hereafter called the summer 
mean monthly wind J. 

(f) Summer mean monthly winds were converted to effective winds be- 
fore the relationships with set-up were determined. 

Eifective winds on Lake Hrie were determined by two methods; 

(a) the first method is based on the assumption that the effective wind 
is a linear function of the vector wind; (b) the second method considers 
the effective wind to be a cosine function of the yector wind. 9 Usain 

the first method, the effectiveness of a given wind was found by: (a) 
subtracting the wind direction from the 100% effective wind direction 

to get the deviation angle, (b) multiplying the deviation angle by 1.11% 
(90° = 100%) to obtain the percentage of ineffectiveness, (c) subtracting 
the’ percentage of inetifeetiveness: from 100% to get the effectiveness per- 


centage, and (d) multiplying the resultant wind speed by the effective- 


ness percentage to get the effective wind speed. For example, adi the 
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direction of the 100% effective wind were 250°, winds from 160° and 340° 
would have zero effectiveness; winds from 70° would be -100% effective; 
and winds from, say, 220° and 280° would be 66.7% effective. Following 
the second method, the effectiveness of a wind was computed by: (a) 
determining the cosine of the deviation angle (equal to 100% effective 
wind direction minus the given wind direction), (b) multiplying the 
cosine by 100 to get the percentage of effectiveness, and (c) multi- 
plying the vector wind speed by the effectiveness percentage to get the 
effective wind speed. 

The resulting summer season over-water mean monthly effective vector 
winds [called effective wind or "V"] were then squared because wind 
stress, and thus setup, is a function of the wind velocity squared 
(see pp. 51-53). The quantity V“ was used in wind-slope calculations, 
eorrelations, etc. 

The effective wind velocities for Lake Erie which were computed 
from the summer mean monthly winds (1950-59) for 14 directions from 251° 
(on axis) to 208° (43° to the right of axis) are given in Table 67 
(Appendix II, p. 228). ;. ‘These effective wind velocities were squared 


and correlated with summer season mean corrected gage differences for the 


same years. 


Correlation of Effective Winds and Wind Slopes 


Observed gage differences are composed of wind slopes and barometric 


pressure effects; therefore, the effect of barometric pressure differences 
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must be computed and subtracted from the gage differences before the 
gage differences can be correlated with the effective winds. Barometric 
pressure effects were computed following the method on p. 229 (Appendix 
II) for Lake Erie (1950-59) and joerg Ontario (1948-55); the results of 
which are given in Tables 68 and 69 (Appendix II, pp. 241-232). 

The wind slopes (set-up) which remained after the barometric pres- 
sure effects were removed from the gage differences were then correlated 
and regression lines determined for each of the 14 directions listed in 
Fig. 8. The regression lines and correlation coefficients were calcu- 
lated by standard formulas (Wallis and Roberts, 1956, pp. 534-545; 


Goedicke, 19553 ps 163). The regression line formula 


Pers a WKY Or yy [= pee (10) 
was used, 
where: 
a = intercept 
b = slope 
x = VW = effective velocity squared 
y = gage difference 
and where 
6 esse ees hi 
Baie (ex)? 
n 
Cer =n - b — dey 
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Correlation coefficients were computed by means of the formula 


where: 


The results of the computation of correlation eee etente and re- 
gression lines using a linear relationship for effective winds and vector 
winds (13 directions) are given in Table 70 (Appendix II, p. 244). Six 
selected correlation coefficients and their effective wind speeds com- 
puted by the cosine method are given in Table 71 (Appendix II, p. 24h). 

A summary of the effective wind directions considered as the 100 per 
cent effective wind, together with their correlation coefficients, is 
given below. After a direction was selected as the 100 per cent effec- 
tive direction, summer mean monthly vector winds for 1950-1959 were con- 
verted to effective winds using the selected direction as the direction 
of the 100 per cent effective wind; these effective wind velocities were 
squared and correlated with observed set-up. This procedure was fol- 
lowed for 13 directions (from a direction parallel to the Toledo-Buffalo 
axis to a direction 43° to the right of the Toledo-Buffalo axis);. the pur- 
pose was to determine the 100 per cent effective wind direction which 
would correlate most closely with the observed set-ups (corrected for bar- 


ometric pressure effect). 


14h 


The fact that a correlation coefficient of 0.76 for ten pairs of 
observations has less than a 1 per cent chance of arising accidentally 
(Herdan, 1960, p. 166) points out the significance of the following 
correlation figures where 11 of the 14 coefficients were considerably 


greater than 0.76. 


100% eff. wind from correl. correl. 
wind — degrees to e0ei i. coeii. 
from Reavoteaxies (cos. ) (linear ) 
2515 OF 0.90 0.84 
248° ae 0.91 0.84 
2ha° of 0.85 
20° ald 0.05 
250" Ae C20K 
245 16° 0.90 
233° 13° 0.92 
EXGN 0} 2a” 0.94 OL95 
228° 23° 0.9% 0.98 
BOGS 28° 0.93 0.94 
altel 342 0.84 
aloe 38° Osa 
208° ls 0.89 0.68 


The correlation coefficients computed from plots of effective winds 
considered as cosine functions of the vector winds and gage differences 
ranges from 0.89 to 0.94; there is little change in the value of the co- 
efficients whether the wind is blowing along the direction of the Toledo- 
Buffalo axis or from a direction 43° to the right of the axis. in this 
case, it would appear that there is no definite 100 per cent effective- 
ness direction within at least 43° of the Toledo-Buffalo axis. On the 
other hand, the correlation coefficients resulting from effective winds 


assumed to have a linear relationship with vector winds shows the corre- 
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lation coefficients increasing to a maximum (the 100 per cent effective 
direction), then decreasing in value as the 100 per cent effective di- 
rection is passed. 

The progressive increase in the correlation coefficient (r) from 
0.84 when the wind is blowing along the Toledo-Buffalo axis direction 
to 0.98 when the wind is from 23° to the right of the axis, followed by 
a decrease in the value of r as the deviation angle exceeds 24°, sug- 
gests very strongly that the wind slope in Lake Erie occurs about 23° 
tO the right of the wind, and, therefore, that the wind-drift surface 
current also flows about 23° to the right of the wind. 

The very high correlation coefficients (up to 0.98) were unexpected 
in light of the errors discussed in pp. 68-76, as well as to the fact 
that wind velocities were determined at only two weather stations (Toledo 
and Buffalo). The only obvious solution seems to be that compensating 
errors have occurred and that vector wind velocities at Toledo and 
Buffalo (when converted with Hunt's ratios) are representative samples 
of the winds which blow over Lake Erie. 

The very close correlation between effective wind velocities and 
corrected gage differences leaves little doubt that the observed gage 
differences are actually wind slopes (net set-up), and consequently are 
not measures of crustal movement. The correlation also brings out the 
fact that the greater part of the error in modern lake-level records 
for Lake Erie, and probably for the other Great Lakes as well, is due to 


meteorological effects. 
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Lake Erie "Hindcast" 


A “"hindeast" of calculated gage differences for 1950-59 was pre- 
pared in an effort to test the validity of the assumption that the ob- 
served gage differences are due to wind slope and barometric pressure 
effect. The wind slope was calculated from the regression line for 
winds from 23° to the right of the Toledo-Buffalo axis (wind slope = 
Pees OX 10.8 v~ - 0.2543) and barometric pressure effects were taken 
from Table 68 (Appendix II, p. 231). ‘The results of these calculations 
are given in Table 6 (Appendix II, p. 179) and are compared graphically 
with observed gage differences in Fig. 9. [Note: The calculations 
were carried to three or four decimal places before the final rounding- 


off to two places in order to reduce the rounding-off error. |] 


Gage Difference 
(Feet) | | | | 
+02 


202 


——— Observed gage differences 


——- "Hindcast" set-up 


ik fe Se | RE eS se eee | | 
1950 1951 1952 1953, 1954 i959 1956 1957 1958 1959 


Fig. 9. Lake Erie "hindcast" set-up compared with 
observed gage differences for the summer seasons 
of 1950-59. Hindcast set-up = calculated wind 
slope + barometric pressure effect. 
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The comparison of observed gage differences and predicted set-up in 
Fig. 9 aids in emphasizing the fact that the observed gage differences 
are due almost wholly to meteorological effects, that the summer season 
mean lake surface is not level, and that uncorrected gage differences 
cannot be used in the calculation of rates of crustal movement. 

+ ( Present Status of Great Lakes Rates of Uplift 

The correlation studies of wind velocities and gage differences to- 
gether with information derived from the isobases of former shoreline 
features permits a resumé to be made of the VelLdiuy, Ola racessote crus tal 
movement in the Great Lakes region which have been calculated up to the 
present time. 

Postglacial uplift around Lake Erie ceased several thousands of 
years ago,and the quantities called "rates of crustal movement" are 
chiefly meteorological effects. Therefore all post-Nipissing rates of 
uplift for this lake are erroneous. 

All of the gage pairs on Lakes Michigan and Huron are either in 
the area of horizontality, or have one gage of the pair in the area of 
horizontality. For this reason all previously determined rates of up- 
lift for these lakes are not valid. Future determinations of possible 
uplift on Lake Michigan and Lake Huron must be determined from gage pairs 
which meet the requirements listed on page 122,and corrected gage records 


must be used for the determinations of uplift. 
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Gage differences, therefore calculated rates of uplift, from gage 
pairs north of the area of horizontality on Lake Ontario are probably, 
for the most.part, the result of meteorological effects; however, lack 
of data prevents a proper correlation study of wind velocity and gage 
differences from being made. When the necessary wind information and 
lake-level information becomes available, it will be possible to deter- 
mine whether gage differences represent net set-up,uplift, or a combi- 
nation of the two factors. 

Most of Lake Superior lies to the north of the Nipissing zero iso- 
base; therefore it is possible that postglacial uplift is occurring 
around Lake Superior. Although uplift may be taking place around Lake 
Superior, the present gage sites are so located that it is impossible 
to determine whether gage differences are caused by uplift, wind slopes 
or-some, Other re hee In the future, if the necessary meteorological 
stations and lake-level gages are established,:.it should be possible 
to separate gage differences into their component parts and to determine 
if one component represents postglacial crustal movement. 

Before accurate rates of postglacial uplift can be measured in 
those areas where it may exist, it will be necessary to revise present 
procedures for obtaining lake-level gage differences so as to eliminate, 
or compensate for, the influences and errors which have been discussed 
in the preceding pages. If the essential meteorological and lake-level 


data become available, and if the necessary precautions are taken and 
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corrections made, then it will be possible to determine whether or not 
postglacial crustal movement exists in the Great Lakes region. In 
addition, if uplift is now taking place in this region, accurate measure- 


ments of its rate may be made. )-) 


VI. CONCLUSIONS 


1. Up to the present the only accurate rates of postglacialieruge a 
movement in the Great Lakes region are those based on the elevations of 
former shoreline features of late glacial and postglacial lakes. Suen 
rates are derived from the differences in elevation between the zero 
isobases and the isobases of maximum deformation—differences which are 
measured in tens or hundreds of feet and which, result from thousands 
years, oi ditierential warping. 

The magnitude of these quantities is great enough for elevations 
determined by ordinary methods of spirit leveling to give accurate rates 
of crustal movement. 

Rates of uplift calculated from isobases of former shorelines are 
to be contrasted with modern rates of uplift for the Great Lakes region 
determined from water-level gage records taken over periods of tens of 
years. In addition, the differences in elevation of ae peas are meas- 
ured in hundredths of a foot, tenths of a foot and, in some cases, im 
Tee. 

2. Modern land uplift due to postglacial isostatic rebound ean 
occur only to the northeast of the Nipissing zero isobase, 1.¢€., norene 
east of the Nipissing area of horizontality. Because postglacial crustal 
movement cannot occur in the Nipissing area of horizontality all rates of 
uplift on Lake Erie are invalid; furthermore, in the other Great Lakes 


all rates of uplift which are based on pairs of gages in which one gage 
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of the pair is south of the Nipissing zero isobase are also invalid. 

3. Owing to a number of reasons (including the failure to consider 
the rise of mean sea level at the tide gage in New York City; mistaking 
the changes of elevation of bench mark Gristmill at Rensselaer, New York, 
which were caused by better determinations of sea level and adjustments 
of the precise leveling net for actual uplift; the inclusion of errors 
which have been discussed in this paper, etc.) the only determination of 
absolute rates of uplift by precise leveling in the Great Lakes region 
is valueless. 

4, The choice of the tide gage at Father Point, Quebec, as the 
initial point for the precise leveling establishing the new International 
Great Lakes Datum was a dubious one. The Father Point tide gage is only 
about 400 miles from the former Laurentide ice divide; thus it is within 
the area of uplift. The gage is at the narrow end of the converging 
Shores of the Gulf of St. Lawrence and the estuary of the St. Lawrence 
River. This situation increases the influence of meteorological effects. 
These factors (if corrections are not made) would prevent precise lev- 
elings made a number of years apart from being compared to the same datum. 

5. Water-level gage records furnishing the basic data for water- 
leveling and modern rates of crustal movement contain errors which reduce 
the accuracy of elevations in the Great Lakes region, and which render 
valueless previously determined rates of crustal movement. The errors 
include those caused by meteorological effects, gage location effects, 


instrument error, and operator. error. Early gage records of the Great 
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Lakes (up to about 1920) contained a greater proportion of instrument 
and operator error than error due to meteorological effects; gage records 
since that time, however, have an increased proportion of error due to 
meteorological effect owing to the reduction of operator and instrument 
error. 

6. The influence of the disturbing factors must be removed or 
rendered insignificant if gage readings are to be an accurate repre- 
sentation of the actual elevation of the lake surface, or if the records 
are to be used in the determination of uplift. 

Reduction in the size of errors due to meteorological efiectstma. 
well as to the effects of gage location, can be brought about by removing 
the water-level gages from their present locations near population cen- 
ters and shallow water, and relocating them at sites in deep water away 
from harbors, bays, rivers, inlets, and areas of man-made changes in the 
configuration of the shoreline and underwater topography. 

7. The assumption that the summer mean lake surface is level (the 
assumption which underlies the practice of waterleveling and the cal- 
culation of modern rates of postglacial crustal movement) is shown to be 
incorrect by the calculation of vector winds from Lake Erie and Lake 
Ontario for the summer months (June, July, August and September) of the 
years 1950-59 (Lake Erie) and 1948-55 (Lake Ontario). These determinations 
reveal a net vector wind from the southwest quadrant—a net vector wind 
which causes a net set-eupbetween water level gages located at opposite 


ends of the lakes. The set-upis expressed as a tilt in the lake surface 
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upward to the northeast. 

8. Rates of crustal movement calculated from best-fit curves of 
gage differences versus time represent, in the large part, the average 
net set-up between the gages of each pair for the time period plotted. 
This is demonstrated by the correlation of effective wind velocities 
on Lake Erie for the summer months of 1950-59 with gage differences of 
water-level gages at Toledo, Ohio, and Buffalo, New York. The corre- 
lation coefficient, r, was 0.98 with a 100 per cent effective wind 
direction from 23° to the right of the Toledo-Buffalo axis. 

9. Definitive correlation studies of wind velocities and gage 
differences necessary for the determination of the existence, or of 
the rates, of uplift cannot be made for Lakes Ontario, Huron, and 
Superior until extensive wind, barometric pressure, and lake-level data 
become available for the eastern end of ee Ontario and the northern 
and northeastern shores of Lake Huron and Lake Superior. When the nec- 
essary data are available, it will be possible to determine whether gage 
Gitzecrences represent: net setup, uplift, or a combination» of these: factors. 

10. Ekman's. classic theory of ocean currents, which represents 
conclusions deduced from a mathematical model of ocean currents incor- 
porating many simplifying assumptions, should be re-examined and modified 
on the basis of empirical observations. This recommendation is supported 
by an analysis of Nansen's original observations (the basis for Ekman's 
theory) which revealed wide discrepancies between the observations and 


the results of Ekman's theory. 
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11. The correlation study of water-level gage differences and ef- 
fective winds on Lake Erie presents an excellent opportunity to test 
Ekman's assumption that the direction of the water-surface slope is al- 
ways in the direction of the wind, and provides a new technique for de- 
termining the angle between the direction of water-surface slope and 
wind direction. 

12. The Lake Erie correlation study (effective winds [computed as 
linear functions of vector winds] vs. gage differences corrected for bar- 
ometric pressure effect) indicated that the observed water-surface slope 
was caused by effective winds whose 100 per cent effective wind was from 
24° to the right of the Toledo-Buffalo axis. 

14. As the wind slope is caused by the wind-drift current, and as 
the wind slope is caused by wind blowing from about 24° to the right of 
the slope direction; then the wind-drift current is also caused by a 
wind blowing from about 23° to the right; or, expressed another way, the 
surface wind-drift current on Lake Erie will be directed about 24° to 
the right of the wind. 

14. The thermocline established during the summer months probably 
acts as a temporary bottom. Therefore, if the thermocline acts as 4 
qQuasi-bottom, and since the depth to the thermocline is approximately the 
same in all of the Great Lakes; then the angle of deviation of surface 
currents from wind direction for Lake Erie (24° to the right of the wind), 
should also be representative of the angle of deviation for the other 


Great Lakes. 
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Cuap. XXH, §§ 1,2.J SPIRIT-LEVELING. 595 


CHAP DER XxX1I. 
ELEVATIONS OF THE GREAT LAKES. 


§ I. The elevations of the Great Lakes above mean tide sea-level, as determined by the Lake 
Survey, depend upon two distinct processes, viz., that of spirit-level measurements and that of 
water-level measurements. 

By the first process, starting from a bench-mark of known height above sea-level at Greenbush, 
New York, the elevation of a bench-mark at Oswego, New York, near the east end of Lake Ontario, 
was found. In like manner the differences in elevation of bench-marks at the following pairs of 
points were determined: Port Dalhousie, Ontario, near the west end of Lake Ontario, and Port 
Colborne, Ontario, near the east end of Lake Erie; Rockwood, Michigan, near the west end of 
Lake Erie, and Lakeport, Michigan, near the south end of Lake Huron; Escanaba, Michigan, near 
the north end of Green Bay,* and Marquette, Michigan, on the south shore of Lake Superior. 

3y the second process, depending on the assumption that the mean surface of each lake is 
level, the relative heights of the pairs of bench-marks for the respective lakes were determined. 
lor this purpose water-gauges were fixed pear these bench-marks, and tri-daily observations of 
the height of the water-surface at each gauge were made during the months of May, June, July, 
and August, 1875, this series of observations being taken as of sufficient extent to give a reliable 
mean. Assuming, then, that the mean surface of each lake for this period of about four months 
was level, the differences of the gauge-readings gave the relative heights of the zero-points of the 
two gauges on each lake, and as these zero-points were carefully referred to the corresponding 
bench-marks, the relative heights of these bench-marks were also known. 

As the surfaces of the lakes vary considerably in elevation from year to year, their mean eleva- 
tions can only be found by observations extending over a series of years. Such observations, con- 
sisting of tri-daily gauge-readings, have been made on Lake Ontario at Charlotte and Sacket’s 
Harbor, NY. Y.; on Lake Erie at Cleveland, Ohio, and Erie, Pa.; on Lake Huron at Port Austin, 
Mich.; on Lake Michigan at Milwaukee, Wis.; and on Lake Superior at Marquette, Mich. By 
comparing the observations made at Oswego in 1875 with those made during the same time at 
Charlotte, the elevation of the bench-mark at the latter place, to which the surface of Lake Ontario 
has been referred, becomes known, and thus also the mean elevation of Lake Ontario for the period 
covered by the observations at Charlotte. Similarly the mean elevation of Lake Erie has been 
derived from the observations made at Cleveland, the mean elevation of Lakes Huron and Mich- 
igan from the observations made at Milwaukee, and the mean elevation.of Lake Superior from the 
observations made at Marquette. 

The methods used and the results derived thereby, of which the foregoing is a brief outline, will 
now be given somewhat in detail. 


LEVELING BY MEANS OF THE SPIRIT-LEVEL. 


§ 2. For this work two parties were detailed, Assistant F, W. Lehnartz having charge of the 
first, and Assistant L. L. Wheeler of the second. During the year 1875 the lines from Greenbush 
to Oswego, and from Port Dalhousie to Port Colborne were leveled in duplicate, and a single line 
of levels was run from Gibraltar, near Rockwood, Mich., to Lakeport. The instruments used 
during this year were Stackpole level No. 1496, 11 inches focal length, object-glass 1} inches in 


* For reasons given in the sequel it is assumed that Lakes Huron and Michigan and Green Bay have the same altitude. 


APPENDIX II 


Table of "hindcast" for Lake Erie. 

Tables of Lake Erie vector wind velocities - 1950-59. 

Tables of Lake Ontario vector wind velocities - 1948-55. 

Tables of Lake Superior vector wind velocities = 1950-59. 

Table of Lake Erie over-water effective wind velocities (linear). 
Procedure for computing barometric pressure effect. 

Table of Lake Erie barometric pressure effects - 1950-59. 

Table of Lake Ontario barometric pressure effects - 1948-55. 


Table of Lake Erie gage differences vs. effective wind velocity squared 
(linear function). 


Table of Lake Erie gage differences vs. effective wind velocity squared 
(cosine function). 
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LAKE ERIE - VECTOR WIND VELOCITIES AT TOLEDO, OHIO - 1959 
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Barometric Pressure Effect 


The magnitude and direction of the tilting of the water surface 
caused by the inequality of barometric pressure at opposite sides of the 
lake were calculated by means of J. F. Hayford's (1922, p. 11) formula 
for the barometric pressure effect. The barometric pressures used in 
the computations were obtained from the U. S. Weather Bureau's Clima- 
tological Data, National Summary (1950-59) for American stations and 
from the Canadian Department of Transport, Air Services (Thomas, 1961, 
letter) for the Canadian stations. 


Hayford's formula is as follows: 


Hh - He = - (Mh - Me)(13.6)(Z5) = - (Mh - Me)(1.13) 
where 
H, = H> = barometric pressure effect, and 
H, = elevation of water at point 1 (feet) 
H> = elevation of water at point 2 (feet) 
M, = barometric pressure at point 1 (inches of mercury) 
Ms = barometric pressure at point 2 (inches of mercury) 


1.13 = 215-6 (density of mercury) , 1. ioe sonar nentoereen) 


1 (density of water) 


Example: —The summer season of 1950 with Buffalo as point 1 and Toledo as 
point 2 gives the following barometric pressure effect: 
Hi - Hs = = (30.012 = 30.025)(1-13) = + 0.015 


In this case, the barometric pressure at Toledo was greater than that at 
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Buffalo by - 0.013 inches; the greater pressure at the western end of 
Lake Erie depressed the water at the western end and caused a rise at 


the eastern end (Buffalo) of 0.015 foot. 
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TABLE 70 


LAKE ERIE GAGE DIFFERENCES VS. EFFECTIVE WIND VELOCITY SQUARED 
(LINEAR FUNCTION) 


Direction of Effective Wind Relative to Toledo-Buffalo Axis 


Gages On Axis Right of Axis 
nivreccnse Effective Velocity Effective Velocity Squared (V' 
Squared (V2) From 11° 
-0.19 2900 4000 
-0.04 8300 11200 
-0.07 TOO 10800 
-0.10 7600 10400 
-0.03 14900 19900 
-0.21 3100 4200 
+0.01 11700 15900 
-0.06 10400 13500 
+0.05 ia 33900 
-0.14 11700 
Correlation 
Coefficient 0.85 0.85 
Regression Line 1.06 x 107? V® -0.185 0.981 x 107? Vv -0.184 0.863 x 107? v= -0.186 0.826 x 107? v= =0.190 
Direction of Effective Wind Relative to Toledo-Buffalo Axis 
Gage® Right of Axis 
ieresenee Effective Velocity Squared (V“) 
From 21° 
-0.19 4200 5300 
-0.04 11900 14600 
-0.07 11400 14100 
-0.10 11000 13500 
-0.03 20700 19300 
-0.21 4500 5500 
+0.01 H 16600 19900 
-0.06 13900 11500 
+0.05 32400 26600 
-0.14 12100 10400 
Correlation 
Coefficient 0.87 0.95 
Regression Line 0.890 x 107? v® -0.201 0.984 x 107? V° -0.217 1.08 x 107? y* -0.152 1.21 x 107? v> -0.248 
Direction of Effective Wind Relative to Toledo-Buffalo Axis 
Gage® Right of Axis 
Di ffenece Effective Velocity Squared (V<) 
From 23° From 38° 
-0.19 7600 
-0.04 13000 
-0.07 20100 
-0.10 12500 
-0.03 12300 
-0.21 4100 
+0.01 14900 
-0.06 7200 
+0.05 16600 
-0.14 6700 


Correlation 
Coefficient 0.98 0.72 
Regression Line 1.25 x 107? Vv? -0.2543 1.37 x 107? V® -0.193 1.32 x 107? V@ -0.166 1.20 x 1079 v® ~0.216 


Direction of Effective Wind Relative to 
Toledo-Buffalo Axis 
Effective Velocity 


Squared (v2) 


From 43° Right of Axis 


Gage® 
Difference 


Correlation 


Coefficient 0.68 


Regression Line 1.29 x 1079 v® -0.207 


8Gage differences are corrected for barometric pressure effect. S 
Note: Effective velocity is in mph/month; therefore v* is in (mph/month)*. 
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